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Abstract 

In  this  article  a  fractional  step  0-method  is  described  and  studied  for  the  approximation  of 
time  dependent  viscoelastic  fluid  flow  equations,  using  the  Johnson-Segalman  constitutive  model. 
The  0-method  implementation  allows  the  velocity  and  pressure  approximations  to  be  decoupled 
from  the  stress,  reducing  the  number  of  unknowns  resolved  at  each  step  of  the  method.  The  con¬ 
stitutive  equation  is  stabilized  using  a  Streamline  Upwinded  Petrov-Galerkin  (SUPG)-method. 
A  priori  error  estimates  are  given  for  the  approximation  scheme.  Numerical  computations  sup¬ 
porting  the  theoretical  results  and  demonstrating  the  0-method  are  also  presented. 
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1  Introduction 


Numerical  methods  for  modeling  viscoelastic  fluid  flow  are  difficult  for  a  variety  of  reasons.  The 
modeling  equations  (assuming  slow  flow)  represent  a  “Stokes  system”  for  the  conservation  of  mass 
and  momentum  equations,  coupled  with  a  non-linear  hyperbolic  equation  describing  the  constitu¬ 
tive  equation  for  the  stress.  The  numerical  approximation  requires  the  determination  of  the  fluid’s 
velocity,  pressure  and  stress  (a  symmetric  tensor).  For  an  accurate  approximation  a  direct  approx¬ 
imation  technique  requires  the  solution  of  a  very  large  non-linear  system  of  equations  at  each  time 
step. 

The  fractional  step  0-method  [25,  23,  24]  decouples  the  approximation  of  velocity  and  pressure  from 
the  approximation  of  the  stress,  thereby  reducing  the  size  of  the  algebraic  systems  which  have  to 
be  solved  at  each  sub-step.  An  added  benefit  of  the  0-method  [25]  is  that  the  algebraic  systems  to 
be  solved  at  each  sub-step  are  linear. 
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The  fractional  step  0-method  was  introduced,  and  its  temporal  approximation  accuracy  studied  for 
a  symmetric,  positive  definite  spatial  operator,  by  Glowinski  and  Pirreneau  in  [11].  The  method  is 
widely  used  for  the  accurate  approximation  of  the  Navier-Stokes  equations  (NSE)  [26,  27,  14].  In 
[15],  Kloucek  and  Rys  showed,  assuming  a  unique  solution  existed,  that  the  0-nrethod  approximation 
converged  to  the  solution  of  the  NSE  as  the  spatial  and  mesh  parameters  went  to  zero  ( h ,  At  — ►  0+). 
The  temporal  discretization  error  for  the  0-nrethod  for  the  NSE  was  studied  by  Miiller-Urbaniak  in 
[18]  and  shown  to  be  second  order.  In  [7]  the  0-nrethod  applied  to  convection-diffusion  equations 
was  shown  to  be  second  order  in  time. 

The  implementation  of  the  fractional  step  0-nrethod  in  [25]  for  viscoelasticity  differs  significantly 
from  that  for  the  NSE.  For  the  NSE  at  each  sub-step  the  discretization  contains  the  stabilizing 
operator  — Au.  For  the  viscoelasticity  problem  the  middle-substep  is  a  pure  convection  (transport) 
problem  that  requires  stabilization  in  order  to  control  the  creation  of  spurious  oscillations  in  the 
numerical  approximation.  Marchal  and  Crochet  [17]  were  the  first  to  use  streamline  upwinding  to 
stabilize  the  hyperbolic  constitutive  equation  in  viscoelastic  flow.  A  second  common  approach  to 
stabilizing  the  convective  transport  problem  is  to  use  a  discontinuous  Galerkin  (DG)  approximation 
for  the  stress  [2,  1], 

Error  analysis  of  finite  element  approximations  to  steady  state  viscoelastic  flow  was  first  done  by 
Baranger  and  Sandri  in  [2]  using  a  DG  formulation  of  the  constitutive  equation.  In  [22]  Sandri 
presented  analysis  of  the  steady  state  problem  using  a  streamline  upwind  Petrov-Galerkin  (SUPG) 
method  of  stabilization.  The  time-dependent  problem  was  first  analyzed  by  Baranger  and  Wardi 
in  [3],  using  an  inrplict  Euler  temporal  discretization,  and  DG  approximation  for  the  hyperbolic 
constitutive  equation.  Ervin  and  Miles  analyzed  the  problem  using  an  implicit  Euler  time  dis¬ 
cretization,  and  a  SUPG  discretization  for  the  stress  in  [9].  Analysis  of  a  modified  Euler-SUPG 
approximation  to  the  transient  viscoelastic  flow  problem  was  presented  by  Bensaada  and  Esselaoui 
in  [4].  The  temporal  accuracy  of  the  approximation  schemes  studied  in  [3,  4,  9]  are  all  O(At). 
Ervin  and  Heuer  proposed  a  Crank-Nicolson  time  discretization  method  [8]  which  they  showed  was 
O  ((At2)).  Their  method  uses  a  three  level  scheme  to  approximate  the  non-linear  terms  in  the 
equations.  Consequently  their  approximation  algorithm  only  requires  linear  systems  of  equations 
to  be  solved. 

In  this  article  we  analyze  a  fractional  step  0-nrethod  for  the  approximation  of  viscoelastic  fluid 
flows.  Advantages  of  the  method  are  three  fold:  (i)  second  order  accuracy  with  respect  to  the 
temporal  discretization,  (ii)  only  linear  systems  of  equations  need  to  be  solved,  (iii)  the  linear 
systems  to  be  solved  only  involve  the  velocity-pressure  or  the  stress  unknowns  (resulting  in  smaller 
linear  systems). 

This  paper  is  organized  as  follows.  In  the  next  section  we  specify  the  problem  and  describe  the 
fractional  step  0-nrethod  for  the  viscoelastic  modeling  equations.  In  Section  3  the  mathematical 
notation  used  is  given.  In  Section  4  we  show  computability  of  the  algorithm,  and  present  the  a  priori 
error  estimates  that  support  the  method.  Two  numerical  examples  demonstrating  the  method  are 
presented  in  Section  7. 


2  The  Mathematical  Model  and  ^-Method  Approximation 

In  this  section  we  present  the  modeling  equations  for  viscoelastic  fluid  flow  as  well  as  a  fractional 
step  0-nrethod  approximation  scheme.  Following  the  description  of  the  0-nrethod,  several  definitions 
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used  to  formulate  the  problem  in  an  appropriate  mathematical  setting  are  given. 


2.1  Modeling  Equations 


The  non-dimensional  modeling  equations  for  a  viscoelastic  fluid  in  a  given  domain  Q  C  (d  =  2, 3) 
using  a  Johnson- Segalman  constitutive  equation  are  written  as: 


Re 


a  +  A 


+  u  ■  Vm 


+  u  ■  Ver  +  ga(cr ,  Vm)  J  - 

2ad(M)  =  0  in  H, 

(2.1) 

1  +  Vp  —  2(1  —  a)  V  •  d(M) 

—  V  •  a  =  f  in  H, 

(2.2) 

V  •  m  =  0  in  H, 

(2.3) 

m  =  0  on  5H, 

(2.4) 

m(0,x)  =  Mo(x)  in  11, 

(2.5) 

cr( 0,x)  =  ctq(x)  in  H. 

(2.6) 

Here  (2.1)  is  the  constitutive  equation  relating  the  fluids  velocity  m  to  the  stress  <x,  and  (2.2)  and 
(2.3)  are  the  conservation  of  momentum  and  conservation  of  mass  equations.  The  fluid  pressure  is 
denoted  by  p.  The  Weissenberg  number  A  is  a  dimensionless  constant  defined  as  the  product  of  a 
characteristic  strain  rate  and  the  relaxation  time  of  the  fluid  [5].  Re  denotes  the  fluids  Reynolds 
number,  f  are  the  body  forces  acting  on  the  fluid,  and  a  G  (0, 1)  denotes  the  fraction  of  the  total 
viscosity  that  is  viscoelastic. 

The  term  ga  and  deformation  tensor  d(u)  are  defined  as: 

ga{cr,  Vu)  :=  (crVu  +  (Vit)T  erj  -  (Vm<t  +  a  (V«)T) 

and 

d(W)  =  ^(vu  +  (Vu)T). 

The  gradient  of  u  is  defined  such  that  (Vm)*j  =  dui/dxj.  For  the  remainder  of  the  paper  slow 
or  inertialess  flow  is  assumed,  allowing  the  term  u  ■  Vm  in  (2.2)  to  be  ignored.  Note  that  when 
the  parameter  a  =  1  in  ga(<x,  Vm)  an  Oldroyd  B  constitutive  model  is  obtained.  For  existence  and 
uniqueness  of  solutions  to  (2.1)-(2.6)  see  [21,  12,  10]. 


2.2  0-method 

In  order  to  implement  a  fractional  step  0-nrethod  for  the  viscoelastic  flow  equations,  an  additive 
decomposition  is  used  for  equations  (2.1)  and  (2.2).  Here  we  introduce  a  splitting  parameter 
u  G  (0, 1),  and  define: 

Constitutive  equation: 

l  Ga  :=  wer,  (2-7) 

2 Ga  :=  (1  —  u)<t  +  A  (u  ■  V<r  +  ga(<7,  Vm))  —  2ad(tt).  (2-8) 
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Conservation  of  Momentum: 


i Fu  '■=  —2(1  —  a)V  •  d('it)  —  V  •  cr  —  f ,  (2.9) 

2  Fu  ■■=  0.  (2.10) 

Let  At  denote  the  temporal  increment  between  times  tn  and  tn+1,  and  for  c  E  {9,oj,a,  a}  let 
c  :=  1  —  c.  Also,  let  /(n)  :=  f(-,nAt). 

The  0-method  approximation  for  viscoelasticity  may  then  be  described  as  follows.  (See  also  [25,  23]). 


0-Method  Algorithm  for  Viscoelasticity 

Step  la:  (Update  the  stress.) 

a(n+e)  _  a(n) 


X- 


6  At 


+  i  Ga{n+e)  =  -2  G<tH. 


Step  lb:  (Solve  for  velocity  and  pressure 

u{n+6)  _  u(n) 


Re- 


6  At 


+  Vp(n+0)  +  i  Fu{n+e)  =  -2  Fu{n\ 
V  •  ■u(n+0)  =  0. 


Step  2a:  (Update  the  velocity  and  pressure. 
u(n+§)  -  u{-n+e"> 


Re- 


(1  -  29)  At 


+  Vp(n+§)  +  2Fu(n+^  =  -i  Fu[n+e\ 


V  •  u 


(n+e)  _ 


0. 


Step  2b:  (Solve  for  the  stress.) 

a{n+~e)  -  a (n+e) 

A  - 


(1  -  29)  At 


+  2  Gcr(n+^  =  -i  Go>+(?). 


Step  3a  and  Step  3b:  In  order  to  complete  the  temporal  advancement  to  time  tn+i, 
Step  la,  and  Step  lb  are  repeated  with  (n)  and  (n  +  9)  replaced  by  [n  +  9^j  and  (n  +  1) 
respectively. 


Note:  Due  to  the  chosen  decomposition  of  the  constitutive  and  conservation  of  momentum  equa¬ 
tions,  (2.7)-(2.10),  the  approximation  of  the  nonlinear  system  (2.1)-(2.6)  using  the  0-nrethod  only 
requires  linear  systems  of  equations  to  be  solved  at  each  step  in  the  process. 


3  Mathematical  Notation 


The  L2(D)  inner  product  and  norm  are  denoted  by  (•,  •),  and  ||u||  respectively.  The  Sobolev  Wp(Q) 
norms  are  denoted  by  |j-||Wfc.  We  use  Hk  to  represent  the  Sobolev  space  Wk,  and  ||-||fc  to  denote 


4 


the  norm  in  Hk .  Function  spaces  used  in  the  analysis  are: 

X  :=  Hl(yi)  :=  {u  G  H\Q)  :  u  =  0  on  30}  , 

S  :=  |<t  =  ((Tij)  :  (Tij  =  a ji  \  <Jij  G  L2(H);  1  <i,j  <  tfj 

n  {cr  =  (cry)  :  u  •  Ver  6  L2(fl),  VuGl}, 

Q  :=  Lo(n)  =  j<7  G  L2(0)  :  J  qdx  =  0  j  , 


Z  :=  <|v  el:  J  q(V-v)dx  =  0,  V<7Gq|. 


Recall  that  the  spaces  X  and  Q  satisfy  the  inf-sup  condition 

(n  Y7  •  v) 

inf  sup  ’  — r—  >  P  >  0. 

<?sQ vex  IMI  Hvlli 


A  variational  formulation  of  (2.1)-(2.3),  found  by  multiplication  of  the  modeling  equations  by  test 
functions  and  integrating  over  11,  is:  Find  ( u,a,p )  :  (0 ,T]  —>  X  x  S  x  Q  such  that 


A  (  -^-,r  j  +  (<t,t)  -  2a(d(n),r)  +  A  (u  •  Vcr  +  ga(tr,  Vu),r)  =  0,  V  r  G  S',  (3.2) 


Re  ,vj  -  (p,  V  •  v)  +  2(1  -  a)  (d(u),d(v))  +  (<x,d(v))  =  (f,v),  Vv6l,  (3.3) 

(V  •  u,q)  =  0,  V  q  G  Q,  (3.4) 

u(0,x)  =  Uo(x),  (3.5) 

cr(0,x)  =  <tq(x).  (3.6) 

As  the  velocity  and  pressure  spaces  X  and  Q  satisfy  the  inf-sup  condition  (3.1),  an  equivalent 
variational  formulation  to  (3.2)-(3.4)  is  given  by:  Find  (it,  cr)  :  (0,  T]  Z  x  S  such  that 


A  (,  ~dt,T)  +  (<T,T)  -  2a(d(it),r)  +  A  (it  •  Vcr  +  ga(cr,  Vu),t)  =  0 ,  V  r  G  S,  (3.7) 


Re(  — ,v  1  +2(l-a)(d(it),d(v))  +  (cr,d(v))  =  (f,v),  V  v  £  Z,  (3.8) 


it(0,  x)  =  u0(x), 
cr(0,x)  =  cr0(x). 


(3.9) 

(3.10) 


To  describe  the  finite  element  framework  let  Tf l  be  a  triangulation  of  the  discretized  domain  11  C  Rrf. 
Then 

n  =  UK,  K  G  Th. 

We  assume  that  there  exist  constants  ci  and  C2  such  that 

c\h  <  hx  <  C2 Pk i 

where  h j.  is  the  diameter  of  triangle  K,  px  is  the  diameter  of  the  greatest  ball  (sphere)  included  in 
K,  and  h  =  max ^7^  hx-  Let  Pk{A)  denote  the  space  of  polynomials  on  A  of  degree  no  greater 
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than  k  and  C(£l)d  the  space  of  vector  valued  functions  with  d  components  which  are  continuous  on 
Cl.  Then  the  associated  finite  element  spaces  are  defined  by: 


Xh  :=  {vGXnC(0)d':v|xGPfc(K)VA'eTft}, 

Sh  :=  {r  €  SnC(n)dxd:r\KePm(K)\/K  eThj, 

Qh  :=  {q€QnC(Si):q\K€Pq(K)VKeTh}, 

Zh  :=  {v  <E  Xh  :  (q,  V  •  v)  =  0  \/ q  €  Qh}  . 


Analogically  to  the  continuous  spaces  assume  that  Xh  and  Qh  satisfy  the  discrete  inf-sup  condition: 

(3-11) 


inf  sup  >  (3  >  0. 

q£QhveXh  ||?||  1 1  v  1 1 ! 


Let  U  and  S  denote  the  L2  projections  of  u  and  a  onto  Zh  and  Sh,  respectively,  and  define: 

A(n)  =  u(n)  _  u{n) ,  E(n)  =  w(n)  _  u  W  ? 


p(")  _  a(n)  _  S(-n\ 
W  _  .An)  _  .An) 


eu  =  u 


h  i 


F(n)  =  S(n)  _ 

=  <T(n)  -  <j (n) . 


We  define  the  discrete  temporal  operator 


dtj(n+1)  ■—  -C(^»+i)  / ifn) 

When  v(x,t)  is  defined  on  the  entire  time  interval  (0,  T), 

IM|oo,fc  :=  sup  ||u(.,t)||fe  ,  |M|o,fc  :=  /  \\v(-,t)\\2kdt) 

o  <t<T  \J  o  / 

The  following  discrete  norms  are  used: 


1/2 


Ml (*)  :=  1HA)II  • 


\v  Lt  :=  max 

°°’A'  l<n<N 


,(n) 


N 


l0,fc  •" 


^n=l 


,(n) 


4  Analysis 

In  this  section  we  investigate  the  numerical  approximation  method  corresponding  to  (3.7)-(3.10) 
First  the  discrete  variational  formulation  of  the  0-method  is  given.  Then  computability  of  the 
algorithm  is  shown  and  a  priori  error  estimates  given. 


4.1  Discrete  Variational  Approximation 

To  stabilize  the  constitutive  equation  a  streamline  upwind  Petrov-Galerkin  (SUPG)  discretization 
is  used  to  control  spurious  oscillations  in  the  approximation.  This  is  implemented  by  testing  all 
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terms  in  the  constitutive  equation  (except  the  discretized  temporal  derivative)  against  modified 
test  elements  of  the  form  TS(n)  where 


T5in)  :=  t  +  5u^’ -Vt.  (4.1) 

Note  that  if  5  is  set  to  zero  the  standard  Galerkin  method  is  obtained.  The  variational  formulations 
for  the  steps  in  the  0-method  approximation  are  as  follows. 

Step  la:  Find  <r^n+^  G  Sh  such  that 

J^t  (CTi"+0)’T)+w  (<rh+e)’T^)  =  J^t  (CT!l'l)’T)-(1-w)  (°’/(r)>T5w)-A  (win)  • 

-  A  (ga{a^\Vu{£)),T5(n))  +2a  (d(M^n)), r5(n))  ,  V  r  G  Sh.  (4.2) 


Step  lb:  Find  uj^+^  G  Zh  such  that 


J^t  («in+e),v)  +2(i  -  a)  (d(«in+e)),d(v)) 


Vv£  Zfr.  (4.3) 


Step  2a:  Find  uh  +  E  Zh  such  that 


(n+0)  \ 

ul  Tv  = 


(1-20)A  t\"h  ,VJ  (1-20)A  t\Uh  ,Vy 

-  2(1  -  a)  (d(4n+0)),d(v))  +  (f("+0),  v) 


Vv£  Zfr.  (4.4) 


Step  2b:  Find  <j\r  +  ^  £  Sh  such  that 


(1  -  2d)  At 


”+">,r)+(l-u,)  (^ 
+  A  (ga(tr  1"+^,  VitN^), 


\  .  (  (n+0)  (n+0)  \ 

>  T5(»+«)  J  +  A  •  VCTh  >  T5(n+§)  J 


T5("+e'))  ~2a  (d(Uhn+*h,TS(n+*) 


(1  -20)A  t 


(n+0)  \  /  (n+0)  \ 

»T)  >r5(«+e)J’ 


V  t  e  Sh-  (4.5) 


Step  3a:  Find  o‘|”+1^  G  <5^  such  that 


/  (n+l)  \  /  (n+l)  \  A  /  (n+0)  \  ,  .  /  (n+0) 

;,r5(„+e-)J  =  —  f  ,r  -(1-w)  (< 


—  A  u 


,(n+^)  Y7„(n+^ 


~Ts(n+§)  A  9o(<  ,V  -uyh  ),  Ts(n+§) 


+  2a  f-d(i*["+^),rs(„+«)>)  , 


V  r-  G  Sh.  (4.6) 
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Step  3b:  Find  6  such  that 


U.++1) 

9 At  \h 


,v)  +2(1  -  a)  (d (i*j?+1)),d(v)^ 

+  (f(n+1),v) 


—  (Jn+<J)  V 

OAt  Uh  ' 


(n+l)  ,/ 

<*h  >d(V 


VvG4  (4.7) 


4.2  Existence  and  Uniqueness 

Before  error  estimates  are  presented  the  computability  of  the  0-nrethod  is  shown.  This  is  accom¬ 
plished  by  proving  that  the  associated  coefficient  matrices  used  in  each  step  of  the  algorithm  are 
invertible. 

The  following  induction  hypothesis  is  used:  There  exists  a  constant  K  such  that  for  n  =  1 , ,N 


u 


(n) 


U 


(n+0) 


<  K. 


(IH1) 


The  justification  of  (IH1)  is  established  below. 

Lemma  1  (Step  la)  Assume  (IH1)  is  true.  For  6  <  Ch  and  At  sufficiently  small  there  exists  a 
unique  solution  o4n+6^  £  Sh  satisfying  (4-2). 


Proof:  Equation  (4.2)  can  be  equivalently  written  as 


A\  ( cr 


(n+6) 
rh  > 


T  = 


A 

6A t 


(ahn)iT)  ~  ((! 


-  A  •  Vo{n)  +^K,Vw[n)))  ,Ts(n)J  +2a  (d(u(hn)), rs(n)^  ,  V r  £  Sh,  (4.8) 


where 


Ali  ( cr^+6\r)  ■=  I  (t 


OAt 


(n+6)  \  /  (n+6) 

h  »T  +wk  ;  >  r<5(") 


Here  (4.8)  represents  a  square  linear  system  of  equations  Hx  =  b.  With  the  choice  r  = 
examining  the  individual  terms  in  Ai  we  see  that 


A 

9At 

ui  I  cr 


.(n+6) 

(n+6) 

h  i 

ah 

(n+6) 
h  ’ 

(n+6) 

"h 

A 

lAt 


(n+6) 


=  LO 


(n+6) 


and 


x  (  (n+6)  (n)  „  (n+6) 

^{<Th  X  -V(Th 


1+5 


—  [u^  ,V  ' ,  cr'h 


(n+6)  (n+6) 


< 

< 


U 


(n) 


lo5 
~2 

codKCh 


h 


Ch 


-l 


(n+6) 


-1 


(n+6) 


Provided  <5  <  Ch ,  and  At  <  (2A )/(9ujKC)  then  A\  <r^n+0^  >  0,  and  thus  the  her  (A\ ) 

{0}.  It  then  follows  that  (4.2)  has  a  unique  solution. 

Lemma  2  (Step  lb)  There  exists  a  unique  solution  u^+0^  £  Z/i  satisfying  (f.3). 

Proof:  Equation  (4.3)  can  be  equivalently  written  as 


A2  (u(hn+9),v)  =  (u^v)  +  ff(n+0),  v)  -  fo-in+0),d(v))  ,  VveZft, 


where 


A2  (4"+0),v)  :=  ^  («in+e),v)  +2(1  -a)  (d(u(™+0)),  d(v))  . 

Note  that  choosing  v  =  u^"'+0'1 

m  (»!,“W),«4"+,>)  =  fft  K”+',,.4"+<’)) +2<!  - «)  (d(l.r,).d(4”+<’)))  >  o. 


and  existence  and  uniqueness  of  a  solution  to  (4.3)  has  been  shown. 


Lemma  3  (Step  2a)  There  exists  a  unique  solution  uh  £  Z \  satisfying  (4-4)- 

Proof:  First  write  equation  (4.4)  as 


^3  = 


u[h+6\  v)  ~  2(!  -  «)  (d(win+0)),d(v)) 


(1  -  29)At  Vh  ’V  ^ 

+  (f(^V)-(<^®\d(v)) 


where 


•43(“"  •V)  ~  (1  -  20)At  >V 


For  the  choice  of  v  —  -u(”+^),  the  system  eC  fit)  +  "f  n)  +  A  >  0,  and  the  proof  is  complete. 


Lemma  4  (Step  2b)  Assume  (IH1)  is  true.  For  6  <  Ch  and  At  sufficiently  small  there  exists  a 
unique  solution  er[  +~^  £  satisfying  (4-5). 


Proof:  Write  (4.5)  as 

A4  (^+9\t)  = 


(1  -  26)At 


(ah+d)’T)  ~u  (*h+e)’T5(n+e))  +2«  (d(u(f+6'>),Ts{n+S)')  , 


A  (4n+V)  ==  (1_2 e)At 


(n+ei)  \  i  t-\  \  (  (n+®) 

h  >T  +(!-w)  > Ts(n+e ) 


+  A  •  V<Tin+0'l’7>+(i))  +  A  [9a{^h+<>\  V«ln+^),Ttf(„+j) 
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Estimating  the  terms  in  A4  cr 


A 


(n+e)  (n+0)\ 


J 


yields 


(1  -  26)A t 


(n+e)  (n+e)\  A 


J  (1  —  26)  At 


(n+e) 


(1  -  w)  (  cr 


(n+e)  (n+0)\ 


'  h  ’  ah 


J  =  (1  -  W) 


(n+e) 


(i-w)  O-^  ,Su 


(n+d)  x„.(n+0~)  .  vff(n+®A 


<  eo5 


+ 


(1  —  te)2  6 

4e0 


(n+e) 


/  (n+e)  (n+e)  (n+e)  A 

A(^«h  -Vo-^  ;  1  <Jh  J 


< 


< 


u 


(n+e) 


Ch 


-1 


(n+e) 


A  KCh 


-1 


A  u 


(n+§)  .T7^(n+0')  X..(n+®)  ■  Vcr^+0^  - 


Vcrl  ,  ATc 


=  A5 


(n+e) 


4n+i>  ■  v4n+i> 


( n+ ~9)  Y7,.(^l+®)^ 


A  (  ga{(T^"'J,Vuh 


J 


Jn+^vJn+^ 

ah  v  uh 


(n+e) 


<  4A 

<  4d1/2A  Vu 

—  ; 

<  Ad1/2\Ch~l  jju 

<  4d1^2XCh~1K 


(n+e) 


(n+e") 


(n+e) 


(n+e") 


(n+e") 


in+§)  Y7„.(n+°)\  X„.(n+°)  .  Vo-(n+^^ 


tyai&h  ~  ’  1 


),<Au 


Thus, 


4A 


(T(n+«)vu(n+9') 


(n+^)  .  Vcr^n+^ 


<  4c/1/2  A 


Vm 


(n+e) 


(n+e) 


du 


X  (n+®)  V7  (n+^) 

<K  •  Vcre 


< 


AdXC2h~2K25 

ei 


(n+e") 


+  Aei(5 


.  v<r(n+») 


Al(<X 


(n+e)  (n+e) 


)  > 


A 


(1  -  26)  At 


+  (l-w)  - 


—XKChT1  (  ^  +  4ci1/2  + 


(1  —  te)2  A 
4e0 

AdKCh~l8^ 

ei 


(n+e") 


+5  (A  -  e0  -  Aei) 


u(n+^  •  Va{n+§) 
ah  v 
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Choosing  eo  =  ei  =  <5  <  Ch  1,  and  A t  <  Ch  1  establishes  A^cr^  +^,cr^  +  ~^)  >  0.  Hence,  a 

unique  solution  exists  for  (4.2). 


The  unique  solvability  of  (4.6)  and  (4.7)  of  the  algorithm  follows  exactly  as  (4.2)  and  (4.3). 

4.3  Error  Estimates 

In  this  section  supporting  analysis  for  the  0-method  described  in  Section  2  is  given.  The  0-nrethod 
decouples  the  stress  updates  from  the  velocity-pressure  updates.  Though  the  modeling  system  is 
nonlinear,  the  updates  for  the  stress  and  velocity-pressure  only  require  linear  systems  of  equations 
to  be  solved.  Motivated  by  this  decoupling  (as  a  first  step  toward  an  analysis  of  the  fully  discretized 
system)  we  investigate  separately  the  stress  and  velocity-pressure  approximations.  First  a  0-nrethod 
for  the  constitutive  equation  given  by  (4.2),  (4.5),  and  (4.6)  is  investigated  assuming  the  velocity 
and  pressure  are  known.  Then  a  0-nrethod  for  the  Stokes-like  problem  given  by  (4.3),  (4.4),  and 
(4.7)  is  analyzed  assuming  the  stress  is  known.  A  priori  error  estimates  for  each  scheme  are  given 
in  Theorems  4.1  and  4.2,  and  a  discussion  of  the  proofs  is  given  below.  (Detailed  proofs  are  given 
in  section  5,  and  6.) 

It  is  convenient  to  define  the  following  notation: 

iih  :=  discrete  approximation  using  true  er, 

&h  ■=  discrete  approximation  using  true  u,  and  p. 


Theorem  4.1  (Assuming  u  and  p  are  known)  For  sufficiently  smooth  solutions  a,  u,  p  such 
that 

Iklloo  JCTtlloo  J^lloo  JWilloo  JVkloo  J(Vkt|loo>  and  ||(V'li)tt||00  <  M,  Vf<E[0,r], 

At  <  Ch2,  the  fractional  step  6-method  approximation,  given  by  Step  la,  Step  2b,  and  Step  3a 
converges  to  a  on  the  interval  (0,T]  as  At,h  — »  0,  and  satisfies  the  error  estimates: 

Ilk  “  ^11  loo, o  ^  Fa(At,h,5),  (4.9) 

Ilk  “  o'hlllo.o  -  Fa{At,h,5),  (4.10) 

and 

/  n 


| tt  •  V(cr  —  cr 


h)\\\e 


u(n+0)  .  v  ( a(n+§)  -  &{n+6^ 


1 

2\  2 


<Fa(At,h,S),  (4.11) 


where 


Fa  {At,  h,  6)  :=  C{At)2  ^  Iktttllo.o  +  lk«ll0)1  +  Ik*  ll0)1  +  lkllo,i 


+  lktt|lo,o  +  Ik*  lino  +  Iklln.n  +  t 


o,o  ^  ii"  no,o 


+  C{At)8  (|kllo,i  +  Ik* llo,i  +  Ikllo.o  +  Ik* llo,o  +  ^t) 


+  C  (, hm+1  +  hm  +  8hm) 


l0,m+l 


+  Chm+1  lktllo.m+1  +  C'tfllktlllo.O  +  Chm+1  IlklllocO  • 


(4.12) 
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Theorem  4.2  (Assuming  a  is  known)  For  a  sufficiently  smooth  solutions  u,  a,  p  such  that 
IMloo  <  M,  Vf  £  [0,  T],  and  At  <  Ch2,  the  fractional  step  6-method  approximation ,  Uh  given  by 
Step  lb,  Step  2a,  and  Step  3b  converges  to  u  on  the  interval  (0,T]  as  At,h  — »  0,  and  satisfies  the 
error  estimates: 

111“  “  “/ill loo, o  ^  Fu{At,h,5),  (4.13) 

and 

\\\u- uh\\\Q1  <  Fu(At,h,5),  (4.14) 

where 

Fu(At,h,5)  :=  Chk+1  K||0)fc+1  +  Chk  |||u|||0)fc+1  +  Ch«+1  |||p|||0i,+1 
+  C(At)2  ||wttt|lo,o  4“  C(At)2  1 1  "*-***  1 1  o,i  4"  C(At)2  ||ft*  llo,o 
+  C(At,)2CT  +  Chk  | Il'tt-H loo.fc+i  •  (4-15) 


Outline  of  the  proof.  In  [6],  and  [7]  the  analysis  for  a  fractional  step  0-method  for  the  convection 
diffusion  equation  was  presented.  The  proof  of  Theorems  4.1  and  4.2  are  done  in  an  analogous 
manner  to  the  proof  presented  in  [7].  Here  we  present  an  outline  of  the  proof  for  Theorem  4.1,  and 
note  that  the  proof  of  Theorem  4.2  is  accomplished  in  a  similar  manner. 


Step  1.  When  obtaining  a  priori  error  estimates  for  the  time  dependent  approximation  schemes  it 
is  useful  to  examine  a  unit  stride,  i.e.  the  terms  analyzed  are  distance  At  apart.  For  the  0-nrethod 
this  is  accomplished  by  considering  linear  combinations  of  the  approximation  methods  steps.  In 


order  to  obtain  unit  strides  from 
following  linear  combinations  are  formed: 


~  ( n )  , 
ah  to 


(n+1) 


~  (n—0)  , 

<  to 


(n+e) 


and  cr^n+0  ^  to 


&^+e)  the 


0Ai(4.2)  +  (1  -  26)  Af(4.5)  +  0Af(4.6),  (4.16) 

0At(4.2)  +  (1  -  26)  Af(4.5)  +  6 At  ((4.6)  with  n  -»•  n  -  1) ,  (4.17) 

0Af(4.2)  +  (1  —  26)  At  ((4.5)  with  n  — ►  n  —  1)  +  6 At  ((4.6)  with  n  — ►  n  —  1)  .  (4.18) 


Step  2.  Evaluate  (3.7)  at  the  midpoint  of  each  unit  stride  and  subtract  equations  (4.16)-(4.18). 
Rearrange  terms  to  obtain  the  equations: 


(n+!) 


—  <T 


(n+1) 


-  a 


r(n) 


—  (T 


(n) 


+  a  tBu 


(n+1) 


—  <7 


(n+1) 


=  a  tBlr 


a.  ~  (n+1)  ~  (n+S)  ~  (n+9)  ~(n) 

At,  <j,  u,  crKh  fieri  >T 


(4.19) 


+  AtB-2pos 


=  AtB2 


rem 


^At,  a,  u,  a 


(n+§) 

h 


~  (n+9)  ~  (n—9) 

ah 


(4.20) 
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(n+e)  _  ~  {n+0)\  _  (  (n+e- 1)  _  ~(n+6- 1) 


,  T  )  +  A tBip 


Jn+6 )  _  ~  (.n+6) 


,T 


=  At%.m  (At,  u, tt'r"-11,  4”+#), 4”>,  r)  ,  (4.21) 
where  Fipos,  B-2pos,  and  Bspo3  denote  the  positive  part  of  the  operators. 

Step  3.  Use  cr =  f  ^  +  F^n\  and  choose  r  =  F^n\  in  (4.19)  to  obtain  an  expression  of 
the  form 


_pi(™+1) 


p(n) 


+  A  tBi  F 


Jn+l)  j_,(n+l) 


<  &tBirem  ( At,  cr,  u ,  T,  F 


(4.22) 


Assuming  that  =  0,  (4.22)  is  then  summed  from  n  =  0  to  n  =  (  -  1  so  that  the  expression 
telescopes  to 


/-i 


+  At£Si- 


n=0 


p(n+ 1)  p(n+ 1) 


=  AtTZ]  ( At,  <x,  it,  f ,  F )  . 


(4.23) 


A  similar  approach  is  taken  with  (4.20)  and  (4.21)  where  r  is  chosen  to  be  Fl'n+0'l  and  F^n+d\ 


respectively,  giving  equations  for  F1'1  ^ 
together  to  form  a  single  equation. 


,  and 


.  These  three  equations  are  then  added 


Step  4.  Suitable  inequalities/estimates  are  then  applied  to  the  terms  in  the  equation. 


Step  5.  Gronwall’s  lemma  is  applied  to  get  an  estimate  for 
using  the  triangle  inequality  we  obtain  the  error  estimate  for 
o-G- 1+0)  _  5-d-1+0) 


i- 1+0 


+ 


J-6 


rG)  -  *(0 


+ 


and  then 


+ 


+ 


In  Section  7  the  numerical  results  presented  use  continuous,  piecewise  linear  approximations  for  <r, 
and  p,  and  a  continuous,  piecewise  quadratic  approximation  for  it.  With  these  approximations  we 
have  the  following  estimates: 


Corollary  1  For  Sh  the  space  of  continuous,  piecewise  linear  functions,  At  <  Ch2,  and  cr,u,p 
sufficiently  smooth,  the  approximation  cry t  satisfies  the  error  estimate: 

1 1| if.  •  V  (cr  —  cf/j)  || \g  <  C  ((At)2  +  A t5  +  hS  +  h  +  5)  .  (4.24) 


Corollary  2  For  Xy,  the  space  of  continuous,  piecewise  quadratic  functions,  and  Qh  the  space  of 
continuous,  piecewise  linear  functions,  At  <  Ch 2 ,  and  cr,  u,p  sufficiently  smooth,  the  approximation 
Uh  satisfies  the  error  estimate: 


I u  -  uh\ 


xl<C((At)2  +  h2) 


(4.25) 
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5  Proof  of  Theorem  4.1 


In  this  Section  the  proof  of  Theorem  4.1  is  accomplished  using  the  steps  outlined  in  Section  4.  We 
define  the  upwinded  test  element  using  true  solution  ul'n'1  as  : 

Tj(„)  :=r  +  «(n*-VT.  (5.1) 

For  notational  simplicity  the  ‘tilde’  is  dropped  from  the  discrete  approximation  using  the  true 
u,  and  p. 


5.1  Set  up  linear  combinations: 


Three  linear  combinations  are  set  up  as  stated  in  (4.16),  (4.17),  and  (4.18). 

Linear  combination  1:  The  first  linear  combination  is  given  by 

0(4.2)  +  (1  -  20)  (4.5) +  0(4.6), 

and  yields  the  espression: 

#i  (<Th,  t)  :=  A  (dt  (o-jT+1))  ,t)  +  0a; 

+  (1  -  20)  u  (a^+0),Ts{n+9))  +  0a;  (cr^+1),  Tg(n+S)'j 

+  Qu  (^T^+e\Ts{n+s^j  +  6u  ,Tsin))  +  0A  •  V(r[n+Vi(n+6-)j 

+0A  (u(n)  •  Vff[n),T{W)  +  0A  (ga{<T^h+e\ 

+0A  (ga«,VMW),Tj(„))  -  (1  -  0)2o  ^d(u(re+0)),r(5(7l+e-)) 

-02a  ^d(u(r,)),T,5(„)^  =0,  V  r  G  Sh. 

The  true  solution  for  the  constitutive  equation  evaluated  at  tn+ 1  satisfies: 


(5.2) 


(5.3) 


-A(t*(n+3)  •  VfT(n+^  +  9a(<T(n+^,Vu(n+^)),r)  ,  Vre5.  (5.4) 
An  appropriatly  upwinded  version  of  the  true  soulution  is  obtained  by  considering 


0  *  (5.4) 


+  0  *  (5.4) 


rs(n) 


s(n+e) 


Next  we  add  and  subtract  terms  to  expression  (5.5)  yielding: 

^(fr,r)=Xf(r),  Vr£S, 


(5.5) 


(5.6) 
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where 


X”(r)  :=  A  [dt^n+l)  -  <rjn+^,r)  +  \Q  •  Vrj 

+  A0  •  Vrj  +  u9  (V n+e ')  -  <r(n+i) , 

+  w0  (o-(n+1)  -  er ("+^),r(5(n+e-)^  +  w(l  -  20)  (Vn+(?)  -  (r(n+5),T5(„+e-)^ 

+  Cj0  (cr(n)  -  cr(n+2),r5(n)^  +  Co9  (cr^n+^  -  a^n+ 2) ,  T^n+S^j 
+  OX  (u^  •  Vo»  -  w(n+D  • 

+  0A  (u(n+e~)  •  Vcr(n+^)  -  u(n+  5)  •  V<r(n+^),r<5(„+e-)) 

+  0A  (ffa(<r(n),VnW)  -  9a(a(w+^),  V>+i)),rJ(„)) 

+  0~A  (ffa(>+0~),  VU("+e')) -ga(a(n+2),Vu(n+12)),Ts{n+S)) 

+  92a  (d(w("+2))  -  d(ti(n)),r5(n)^  +  92a  (d(-u("+2))  -  d(it(n+^)),  T^n+S^  .  (5.7) 

Subtraction  of  (5.3)  from  (5.6)  gives: 

-Bi(ah,r)  =Xf(r),  V  r  <E  Sh,  (5.8) 

the  first  expression  that  will  be  bounded  to  complete  the  analysis. 

Linear  combination  2  :  The  second  linear  combination  gives  a  full  discrete  time  step  of  size  At 
between  t,n_Q  and  tn+i-6»  once  we  multiply  the  resulting  expression  by  At  in  a  later  step.  Using 
(4.17)  divided  by  At  yields, 

B%{(Th,  t)  :=  A  [dt  ,r^j  +9uj  (o-[n+0),  t5(„))  +  (1  -  20)  u  (a^+e) ,  T5{n+g)^j 

+  0uj  (°{n)> T, 5(«-»))  +  C1  -  2d)  &  ^in+^ ,  T5(„+e-)'j  +  0c5  (o{n),  t5(„)^ 

+ ^  ^>r5(«-e)) + (!  - 26))  a  •  Vffln+Vj(n+9-)j 

+  0A  («W  •  V<r ^ ,  r5(n) )  +  0 A  •  Vo-J^ ,  rg(n-8) ) 

+  (1  -  20)  A  ^3a(o,in+^,  V«^),T{(n+<-)j  +  0A  ((/«(4n),V«(nVi(")) 

+  0A  (ffa(ff[n^),VM(,i^)),rJ(„-C))  -  (1  -20)  2a  (d(-u(n+0)), 

-  02a  ^d(w(n^),  T5(n)^j  -  02a  ^d(ri(re_e,)),  r5(n-e)^  =0,  \/  r  £  Sh.  (5.9) 

The  true  solution  to  the  constitutive  equation  evaluated  at  the  midpoint  between  tn-o  and  tn+i-e» 
is 

°  =  "A  [^dt  e),T)  ”  (^n+h~e)^)+2a  (d (U("+M)),T) 

-  A  (u(n+2~e)  .  Vff(n+3-fl)  +  ga((r(n+2-e) ,Vu(n+2-9)),r^  ,  VreS.  (5.10) 
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In  a  manner  similar  to  (5.5)  we  form  the  linear  combination 


0(5.10) 

+  (1  -  20)  (5.10) 

+  0(5.10) 

Ts(n) 

"^(n+e) 

and  write  the  expression  as 


%r)=^(r),VTeS, 


(5.11) 


(5.12) 


where 

2?(r)  :=  A  (dtcr(n+d)  -  a^~e\r^  +  A  0  ,  Su™  •  Vrj 

+A(1  -2 0)  (^-(r{n+^~e\5u^n+^-e)  -Vi-j  +  A0  ■  Vr) 

+9u  (<x(n+0)  -  >+^),T5(n))  +  (1  -  20)u;  (V"+0)  -  >+^0) ,  T^)) 

+0w  (<x(n)  -  o-(n+5_0),T5(„-e))  +  (1  -  2 0)  w  (<r(n+0)  -  a-(n+5-0),r5(n+e-)^ 

+0c5  (crW  -  <>+M,t5(„))  +0u  (a^-V  -  *(n+*-e) , 

+  (1  -  20)  A  (u(n+fl)  -Vcr(n+0)  --u(n+i-0)  •  V<r(n+5-e),T(5(„+e-)) 

+0A  (u(")  •  Vf rW  -  w(n+3-0)  •  Vfr(n+3-0),r5(„)) 

+0A  (u(n-0)  •  Vcr(n-0)  -  w(n+|-e)  .  Vo-(n+3-0),T5(„-e)) 

+  (1  -  20)  A  (^(crM),  VuM)  -5a(^(n+^0),VW("+^-0)),r5{n+e-)) 

+0A 

+0A  [ga{cT^-e\yu^)  -^(^("+3-®),  Vu^+s-®)),^-,)) 

+02a  (d(W(n+3-e))  -  d(ti^),r5(n))  +  (1  -  20)  2a  (d(>+3-0))  -  d(tx(n+0)), rtf(n+S)) 

+02a  (d(u(n+3-0))  -  d(M(n-e)),r5(n-e))  .  (5.13) 

An  error  expression  is  now  formed  by  subtracting  (5.9)  from  (5.12): 

#2  (<r,  t)  -  B%(crh,  r)  =  2^(r),  V  r  G  SV  (5.14) 

Linear  combination  3:  We  form  the  third  linear  combination  (such  that  when  multiplied  by  At 
a  full  discrete  time  step  between  tn+o~  l  and  Ui+0  nray  be  analyzed)  by  considering 

0  *  (4.2)  +  (1  —  20)  *  ((4.5)  with  n  — >  n  —  1)  +  0  *  ((4.6)  with  n  — >  n  —  1) .  (5.15) 

Using  (5.15)  we  obtain 


B^(ah,r)  :=  A  (o{n+e))  ,r)  +  0w  (o|n+<?),  r^))  +  0w  (o{n),  Ts{n-e)^ 

+  (1  -  20)  w  (crd1+0_1),r(5(„-e)^  +  0w  (cr^,  t5(„)^  +6lj  (cr^~d) ,t5(„-o)) 

+0 A  (rx^-0)  •  Vcr^-0) ,  r6(„-e) )  +  0 A  («W  •  V«r  ,  r5 (n) ) 

+0A  (ga(a-^~6\  VM(n“0)),r^-9)j  +  0A  (#a(cr£,  Vu(n)), 

-02a  ^d(M(n));T(5(n)^  -  02a  ^d(«(n-0)),  rS(„-e)^j  =0,  V  r  £  S;,.  (5.16) 
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Evaluating  the  true  solution  at  the  midpoint  between  tn+e- 1  and  tn+g: 


0  =  -A  2),rj  -  (>+0-D,r)  +2a  (d(t>+<?-§)),  r) 

-  A  (u(n+0-^  -vAn+e-^)  +  ga{a(n+e~^ ,Vu(n+9~12)),T^  ,  V  r  G  S.  (5.17) 
An  appropriately  upwinded  true  solution  is  found  using  the  linear  combination 

0(5.17)  +0(5.17)  .  (5.18) 

Ts(n)  Ts(n-8) 

Adding  B^(a,r)  to  both  sides  of  expression  (5.18)  gives 

^(<t,t)=X3"(t),  Vt  G  S,  (5.19) 

where 

2£(t)  :=  A  +  A0  5«(n)  ■  Vrj 

+A0  ^a[n+e^\su^  •  Vrj  +  6u  -  <t(™+0-§),  T(5(n)) 

+0cu  (o»  -  a(n+e~12),rs(n-e))  +  (1  -  20)  w  (V^1)  -  a(n+e~^ ,Ts(n-e)) 

+6u  (cr^  -  <r(n+e~2))T(5(ll)^  +0+  -  a(n+e~^ ,  T(5(„-fi)) 

+0A  («(""*)  •  V«r(n“®>  -  w(n+0-i)  •  V(rH4);T(5(„_,)J 
+0A  (u W  •  Vtr(n)  -  w(n+0-i)  • 

+0A  (ga(*(n-d\Vu(n-^)  -  <7„(<r(n+^),  VUH-^)),T{M) 

+0A  (5a(<r",VMW)  -  5a(>+^),  V>+^)),r5(„)) 

+02a  (d(u(n+9~^)  -  d(w^),r5(„)^ 

+02a  ( 'd(u(n+d~ s))  -  d(M(n-^),T5(„-e))  .  (5.20) 

Subtraction  of  (5.16)  from  (5.19)  yields  the  final  error  expression  needed  for  analysis  of  the  consti¬ 
tutive  equation 

#3(cr>T)  -#3 (°+>r)  =  J3(t)5  Vt  G  Sfc.  (5.21) 

Bounds  are  now  found  for  (5. 8), (5. 14),  and  (5.21). 

5.1.1  Bounding  B™(<t,t)  —  i3”(crfe,r) 

In  this  section  bounds  on  the  term  B™(cr,  _p(n+1))  _ (0+,  i/’(n+1))  in  expression  (5.8)  are  examined. 
Recall  for  this  proof  we  have  dropped  the  ‘tilde’  from  <x h-  Thus, 

rM  =<r(n)_5(n))  and  F(n)  =5(n)  _a.(n)_ 
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It  is  also  convenient  to  use  the  notation: 


:=  FM  +  5u^  -VF(u\ 


Here, 


^(o-,F("+1))  -  B^(ah,F^)  =  A  ( dt(e^+1)),F^+ x))  +  9u  [e^+6\  F^) 

,  ft.  ,  (Jn+1'>  p(n+1)  A  I  a  na\  ,  ,  (An+e)  jp{n+l)  \ 

+  0UJ  ^eCT  ,  Fg{n+S)  J  +  (1-26 i)u  yefj-  ,  F &{n+-e) ) 


+  OiJ  [  6 


,(n+0)  pin+i) 


+  dCo  leZ’F 


(n)  p,(n+l) 


a  >  -1  $(n) 


+  ex  (u(n+»)  ■  Ve^+6\F^H^  +  ex 


?>)  •  Ve(n)  F(n+1) 
a  v  (  a  ,  r  £(n) 


+  6>A  ga(e 


+  ex(ga(ei\Vu^),F^ 


Each  of  the  terms  in  (5.23)  is  now  bounded. 


where 


>  —  f  _p(n+1)  “  _  p(n+1)  p(n) 
~  At  \ 

>  A  /"  F'(n+1)  2  _  r(«)  A 

-  2At  v  r 


(5.22) 


(5.23) 


A  (dt(e^+1)),  F(n+1))  =  A  (dt(F(”+1)),F{"+1))  +  A  (dt(r(n+1)),  F(n+1))  ,  (5.24) 

A  (dt(i'1(n+1)),  F(n+1)^  =  —  ^F^n+1\  F^n+1^  — 


(5.25) 


A  (dt(r(n+1)),F(n+1))  <  ^  di(r(n+1))  2  +  a  F(n+1)  2. 

Each  of  the  next  eight  terms  in  (5.23)  is  rewritten  as  the  sum  of  four  terms.  Starting  with 

du  (e^+0),  F^+1})  =  Oco  ^n+0\F^n+1^  +  Oco  ^n+0\6u{n^  •  VE(n+1)) 

+6uj  (r^n+e\  F{n+1^  +  dio  (r^n+e\du^  •  VF(n+1))  . 


(5.26) 


we  have  the  bounds 


9uj  ^F('n+e'1 ,  F^n+1^  <  Ouj  F^n+6>)  _p(«+1) 


<  —  pin+e)  2  ,  ^  p(n+ 1) 

“2  2 


0cj  (V(re+e),<5w(n)  •  VF(n+1)) 


<  Quo  F{n+0)  5u{n)  •  VF(n+1) 


<  Sdu  F(n+0)  u{n)  $  VF(n+1) 

OO 

<  60uM$  F{n+0)  VF(n+1) 

<  6ujM2^  \\F(n+e) II2  +  52Qi0  II  vir,(n+1) II 2  j 


(5.27) 


(5.28) 


(5.29) 
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tfw(r(n+0\F(Tl+1))  <  ^  ||r(n+0)||2  + ^  ||F(n+1)||2,  (5.30) 

and 

9uj  (r{n+e\5u(n)  •  VF(n+1))  <  r(n+,9)  2  +  520uj  VF(n+1)  2.  (5.31) 

Splitting  the  next  term  into  its  pieces  gives: 

0lu  (e£+1),  =  9uJ  (F(n+1),  F<n+1))  +  9uj  ^n+1\duA+e)  .  VF(n+1>) 

+0u;  (r(n+1),  F(n+1))  +  6»w  (T(n+1),  <5w(n+e)  •  VF(n+1))  .  (5.32) 

Bounding  each  term  we  have 

9lo  ^n+1\  F{n+1^  =  9lu  F(n+1)  2,  (5.33) 

9u  (F{n+1\5u(n+°)  -VF(n+1))  <  6uJ^2d  F(n+1)  2  +  520cj  VF(,i+1)  2,  (5.34) 

(r(n+1),F(n+1))  <  ^  r(n+1)  2  +  y  F(n+1)  2,  (5.35) 

and 

flea  (r(n+1\^(n+e~)  -VF(W+1))  <  9uJ^d  r(n+1)  2  +  529uj  VF("+1)  2.  (5.36) 

Next 

{l-29)u(e^+0\F^H))  =  (1-2 9)  u  (f^+6\  F^) 

+  (1  -  29)  lo  (F(n+0),S uA+'d)  •  VF(n+1)) 

+  {1  -  29)  lo  (r^n+e\  F^n+1^ 

+  (1  -  20)  lo  (r^n+e\  6u(n+e)  •  VF(n+1))  .  (5.37) 

Bounding  each  term 

(1  -  29)  u  (V n+9\  F(w+1))  <  ^  ~ -  F^n+9^  2  +  ^  ~  U  F("+1)  2,  (5.38) 

(1  -  29)  u  (yF[-n+e\5uA+~e)  •  VF(n+1))  <  (x  ~  2(9j  ujM29-  F(n+o)  2 

+  (1-2 9)uj52  VF(n+1)  2  (5.39) 

(1  -  20)  lu  (r(n+g),  F(n+1))  <  ^  ~  2^  ^  r(n+6,)  2  +  zAAbi  F(n+i)  2  (5.40) 

(1  -  20)  co  (r(n+9\  5u(n+e)  •  VF(n+1)j  <  (x  ~  2^| u M  9  r(n+e )  2 

+  52{l-29)uj\\vF{n+1)\\2  (5.41) 
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Splitting  the  following  gives: 


04)  (e^+0\  =  k o  (^jp(r!+e'),F(n+1^  +  04)  ,5u(n+°)  ■  VF{n+l^ 

+  fa  (r(n+0) ,  F(n+1))  +  eco  (r(n+0) ,  5u(n+e )  •  vf("+1))  , 

and  these  terms  are  bounded  as 

04)  (V(n+e~),F(n+1))  <  ^  f(”+0')  2  +  |  F(n+1)  2,  (5.42) 

0£>  (F(w+e~),^(n+^)  •  VF(n+1))  <  du}^  d  F(n+§)  2  +  S2ecj  VF("+1)  2,  (5.43) 

eco  (r(n+e~),F(n+1))  <  ^  r(n+e')  2  +  ^  F(n+1)  ^  (5.44) 

and 

0£  (r("+*),<fa(n+*)  -VF(W+1))  <  r(n+0')  2  +529Cj  VF(n+1)  2.  (5.45) 

Continuing 

04)  (eff,  F^+1})  =  04)  (FW,F(n+1))  +  04;  (fW,  £u(n)  •  VF(n+1)) 

+04)  (r(n),  F(n+1))  +  04)  (V(n),  <5u(n)  •  VF(n+1))  .  (5.46) 

These  terms  are  bounded  as: 

9u[F(n\F^n+1^  <d-^~  F™  2  +  y  F(n+1)  2,  (5.47) 

04)  (f(w),  Ww)  •  VF(n+1))  <  0ujA^^'  pW  2  +  <5204)  VF(n+1)  2,  (5.48) 

04;  (r(n\F(n+1))  <  ^  ||r(n)||2  +  y  ||-^(n+1)||2,  (5.49) 

and 

04)  (r(n),Wn)  •  VF(ri+1))  <  0ujA*2(*  r(n)  2  +  529Cj  VF(n+1)  2.  (5.50) 

Next  the  convective  terms  are  handled.  Thus, 

0A  ^w(n+^)  •  =  0A  («(n+^)  •  VF^n+§) ,  F(n+1)) 

+  0A  (u(n+0)  •  V F(n+d) ,  5u(n+e)  •  VF(n+1))  +  0A  (u(n+e)  •  vr(n+0),  F(re+1)) 

+  0A  (u(n+e">  •  Vr(n+e) ,  •  VF(n+1))  .  (5.51) 

We  bound  each  part  in  the  following  manner: 

0A  (u(n+0)  •  VF(n+0),F(n+1))  <  9Xd^L  i?(«+ !)  2  +  0A  VF(n+0)  2,  (5.52) 
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e\  (u(n+^)  •  VF(n+^,k(”+e')  •  V.F(n+1))  <  50X(^M  VF(”+1)  2 

r  /i  \  r  7i  r9 


+  %!||V^)f,  (5.53) 

<9A  (n(n+g)  •  vr(n+e),F(n+1^  <  F(n+1)  2  +  0A  vr(n+0)  2,  (5.54) 

ex  (u(n+^  ■  vr(n+°) ,  8u(n+^  •  vf("+1))  <  ^  exdM  vr(n+0)  2 

7  A  ^  t  rO 


2  I 
OXdM2 


Similarly  the  second  convective  term  is  split  as: 
OX  (u (">  •  Ve^}, F^+1})  =  0A  (w(n)  •  VF(n),. 


VF(n+1)  2  .  (5.55) 


+  0X  ( u W 


.  (u(n)  •  Ve^}, F^+1})  =  6X  (u (n)  •  VFW,F("+1) 

+  ex  (u(n)-vrW,F(n+1))  s 

and  bounded  term  by  term  with 

0A  («(")  •  VFW,F("+1))  <  0X^2 

ex  (u^  •  vF(n),<jM(n)  •  vF(n+1))  <  S0X^M2  1^ 
0A  (u(n>  •  vrW,F("+1))  <  gA(^M2 


1})  +  ex  (uin)  •  VF(n\5u{n)  •  VF(n+1)' 
+  0A  (w(n)  •  Vr(n), 5u(n)  •  VF(n+1))  , 


(5.56) 


F(n+1)  2  +  Q\  VF(n) 


(5.57) 


56XdM2  n  2  SOXdM2  „,nx  2 

VF(n+1)J  <  - - -  VF(n+1)  + - - -  VF(,i)  ,  (5.58) 


ex  ^(n)  •  vr(n),F(w+1))  <  exd^1  F{n+1)  ~  +  ex  vr(n)  2 

and 

ex  (u (n)  •  vr^,  W")  •  vF(n+1))  <  ex ^m2  vf("+1)  2  +  6~ex^M~  VI 

Lastly  the  ga(-,-)  terms  are  considered.  Splitting, 

ex  (gaiety, Vu(.n+°)),F^H^j  =  ex  (ga(F(n+e~),Vu(n+iP>),F(n+1^ 

+  ex  (ga(F(n+°),  Vw(n+0')),5ti(n+^  •  vF(n+1))  +  0A  (so(r(n+®),  Vt*(n+I 


„  ,„mA  eXdM2  „  ,n  2  52eXdM2 

■  vf(  +1)  i  < _  v F'  +1)  h _ vr*-  -1 

)  -  2  2 


(5.59) 


(5.60) 


/ 

/■Vf  y  +0A  ^a(r(n+0'),vJn+e')),F(n+1)y 
+  0A  ^a(r(n+0'),  V«(n+®)))(y«(n+®)  •  VF(n+1))  ,  (5.61) 


and  bounding  each  of  these  terms  gives: 


wr  unroc  icimo  ^ivco. 

0A^a(F(n+e'),Vu(n+0)),F(n+1))  <  0a|  ^(fM.VmH')) 

<  M0x\\vu(n+^\\  II  F(n+ 


<  4d0X  Vu(n+®)  F(n+e') 

oo 

<  AdeXM  F^n+e)  F(n+1) 


_p(n+!) 

'  ||.Fln+1)| 
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eXlqJF 


■  VF(n+1A  <  0X5 


Ad2  0X5  Vm 


4d20X8M 2  F(n+0)  VF("+1) 


w(n+e)  .  y_p(ra+l) 


i?1(Tl+0)  V 


<  4d40XM4  F(n+0')  2  +  520A  VF(n+1)  2 ,  (5.63) 


(5.64) 


,  f("+1)N)  <  40AdAf  r(n+e')  2  +  6XdM  F(n+1)  2, 


Similarly, 


n+^)  •  VF(n+1))  <  Ad49XM4  r(n+^)  2  +  J20A  VF(n+1)  2.  (5.65) 


0A  ^(e^VtiW)^'*;1))  =  6X  (ga(F(n\VuW),Ftn+1^ 

+  9X  [ga(F(n\  Vu(n)),<5u(n)  •  VF(n+1))  +  0A  ^0(r(n),  V«(n)), 

+  0A  ^o(r(n),  V«W),  <y«(n)  •  VF(n+1))  ,  (5.66) 


and  has  the  following  bounds: 


9x(ga{F(n\VuW),F(n+1^  <40XdM  F{n)  2  +  OXdM  F(n+1)  2,  (5.67) 

6x(ga{F{n\Vu(n)),5u(n)  -VF(n+1))  <4d40XM4  F(n)  2  +  6>A<52  VF(n+1)  2,  (5.68) 

6x(ga(T{n\Vu^),F(n+1^  <40XdM  T(n)  2  +  OXdM  F^n+1)  2,  (5.69) 

0A(sa(r(n),Vu(n)),£u(n)  -VF(n+1))  <40Xd4M4  T(n)  2  +  <S20A  VF(n+1)  2.  (5.70) 


5.1.2  Bounding  £>2(ct,t)  —  B^^cry^T) 

We  continue  in  this  section  by  determining  bounds  on  the  terms  B^icr,  r)  —  £>2(07^  t)  in  expression 
(5.14).  Note  that  the  value  of  r  is  chosen  to  be  _p(n+^).  Thus, 

f(”+«))  -  BS(ah,  f("«))  =  A  (dte!"+‘,),F("+>>))  +  fc  ^+"). 

+  (l-2e)Je^>,d;::}))  +eJeff,F^  +(1  -  We^F^) 


(n)  (n+0)A  /  (n+fii)  (n+f?) 

cr  >  -F  ,5(„_9)  )  +  (1  it/Jwleg-  ,  -*5(„+e) 


(4n).-r)w#))  +  »  +  (1  -  2«)  A  (>+»)  ■  Ve^.F^ 


<T  >  '  fi(n-9) 


+0X  .  Ve^,  F^})  +  OX  (u^  • 


VT  Jn-e)  w(n+0) 
VeCT  1*  g(n-9) 


+  (1  -20)  X  ^(e^.V^^),^)  +0A  ^(e^.Vi*^),^^) 

+^A  L(e&-0),Vu^),F^+_^  • 


(5.71) 
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Each  term  in  (5.71)  is  now  bounded.  Here 


A  (A  ,  F(n+^  =  A  (dt  (f("+®))  ,  F(n+^))  +  A  (dt  ( T 


with 


and 


A  (dt  F 


,F 


> 


A  dt.  r 


2A  t 

F(n+e)^  <  ^  ^(n+e) 


^(n+0)  _  p(n-6) 


2  A 
+  4 


d*  r 


Then 

0u;(e^,F-  7  I 


(n+0)  =  0W  [_p(n+0)i  F(n+e)\  +  0(J  /'_p(n+e)>  5u(n)  .  V_p 


+  fr(n+0),  F(n+e')N)  +  duj  (r(n+e\  5u (**)  • 


and  the  terms  are  bounded  as: 

9u  (F^n+e\F 


< 


9id 


pin+e) 


did  (F^n+e\6u{n)  -VF(n+fl))  < 


2 

9uM2d 


2  Odd 

+  T 


p(n+e) 


pi(n+e) 


Odd  fr  (n+e\F 


< 


Odd 


p(n+0) 


+ 


+  520ud 

II  F(n+0) 


V_p(n+^) 


and 

9u  (V  (n+e),Wn)  •  VF 
Splitting  the  next  term: 

An+e)  „(n+e)\ 


< 


9udM2d 


p(n+(9) 


+  52did 


VF(fl+fl) 


(1  -  29)  u  [e(Z+e\F[*J^  =  (1  -  29)  u,  [f^+9\  F 

+  (1  -  29)  u;  (Vn+0),  Su(n+9)  ■  VF 
+  (1  -  29)  u  (r(n+(9),  F(n+^))  +  (1  -  29)  u  (r^n+e\du(n+9)  •  VF 

Bounding  these  terms: 


(1-2  9)dd(F^n+e\F 


< 


(1  -  29)  u: 


p(n+9) 


2  (1-2 9)  id 

+  -T±- 


p(n+e) 


(5.72) 

(5.73) 

(5.74) 

(5.75) 

(5.76) 

(5.77) 

(5.78) 

(5.79) 


(5.80) 

(5.81) 


(1  -  29)  id  ^n+d\5u(n+§)  •  VF(n+e')) 

(1  -  2 9)  id 


< 


4ei 


(1  -  29)  id  (r(,l+9\  F(n+0^  < 


j?(n+0) 


(1  -  2 9)  id 


+  <52  (1  -  29)  W£! 

p(n+0) 


U 


(n+^)  .  v f(,1+®) 


2+  (l-20)w 


p(n+e) 


,  (5.82) 
(5.83) 
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and 


(1  -  29)  uj  (Y(n+0\  8u(n+°)  •  VF(n+^)) 

<  ^  ~2e^u  r(n+0)  2  +  <52  (1  -  29)  wei  «("+*)  •  Vf(,1+®)  2  .  (5.84) 

4ei 

Next 

9uj  (e{^\  =  9uj  (V n\  f(h+9^  +  9uj  (f™,  Su^n~^  •  Vf('1+^ 

+  9u  (r(n),  F<"+0))  +  9uo  (^r(n),  5w(n_e)  •  VF(n+0))  ,  (5.85) 

and  each  term  is  bounded  as 

9u  (FW,F(n+0'))  <  ^  ||F(n)||2  +  Y  |F^||2’  (5.86) 

9uj  (Vn),<5-u(n~e)  •  VF(n+e))  <  6ujA^  - -  F(n)  2  +  529lu  Vf(,,+9)  2,  (5.87) 

0w(r(n\iKn+*))  <  ^  ||r(n)||2  +  Y  |[F(n+^||2,  (5.88) 

and 

9co  (r(n),<5-u(n~e)  •  VF(n+e~))  <  0ujA^  r(n)  2  +  629lu  VfM')  2.  (5.89) 

The  term 

(1  -  2 9)u^+6\f^+6^  =  (1  -  2^(fH1,fM) 

+  (l-29)u(F(n+e\5u(n+e)  ■VFf’1^))  +  (1  -  26»)c5  (r(n+0),F^+^) 

+  (1  -  29)  Cj  (r(n+e) ,  5u(n+e)  •  VF(n+(f))  ,  (5.90) 

is  bounded  as 

(1  —  29)  u  ^T’(n+61) ,  _F(n+e))  =  (l  —  29)  u  F^n+e)  2,  (5.91) 

(1  -  29)  Q  [F^n+e),  5u(n+9)  ■  Vf('1+9)) 

<  (1  ~  26)UJ  f(u+9)  2  +  (1  -  29)  wei  «("+*)  •  VF(,1+®)  2  ,  (5.92) 

4ei 

(l-29)Cu(r(n+9\F(n+9)^  <  ||r(n+e~)||2+  ||jF("+e)||2 ,  (5.93) 

and 

(1  -  29)  u  (r(n+e),5u(n+e)  •  VF(,1+9)) 

<  (1  ~  r(n+^)  2  +  52  (1  -  29)  u>6!  u(n+9)  •  VF(n+9)  2  .  (5.94) 

4ei 
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The  term 


6Cj  (e{£\  F^J^  =  OCo  ( F(n\  f(”+0))  +  9Cj  (f< n\8u W  •  VF(n+0)) 

+  ecu  (r w,  +  eco  (r(n\  <Wn)  •  vF(n+0')) 

with  bounds 

(F(n\F("+0))  <  ^  F(n)  2  +  y  F(n+0)  \ 

9u  (f< n\5u (n)  •  VfH'))  <  0U}M  d  F(n)  2  +  62 9~  VF(n+e) 

eco  (r(n),  f("+0))  <  ^  r(n)  2  +  ^  F(n+0)  2 , 

and 

eco  (r(n),fo/n)  -vfM)  <  9ujA^  ^  r(n)  2  +  629Cj  VF(n+0)  2. 

Considering, 

9Co  (e^~0),  F^t^  =  0£o  (Vn“0),  F(n+e'))  +  0w  (f^”0),  5u^n~^  •  VF(n+0)) 

+  (r(n_0>,F(n+0))  +  du>  (r(n-e\ 5u(n~e')  •  VF(n+0))  . 

We  have  the  bounds 

6£j  ^F(n_0),  f("+0)^  <  —  F(n_0)  2  +  —  _p(n+^)  2; 

9u  (Vra~0),  5u(n~e'>  .VfM)  <  6ujA^  d  F (n“0)  2  +  526i(h  VF(n+0)  2, 

6>o>  (r(n_0),F(n+0))  <  ^  r(n_0)  2  +  y  F(n+^  2, 

and 

go;  (r(n~0),<5tt(n~0)  •  vf(n+^)  <  6ul^  ^  r(n_0)  2  +  52ocb  vf<t+0)  2. 

Examining  the  convective  terms: 


(5.95) 

(5.96) 

(5.97) 

(5.98) 

(5.99) 

(5.100) 

(5.101) 

(5.102) 

(5.103) 

(5.104) 


(1  -  29)  A  ^w(n+0)  •  Ve^+^F^H  =  (1  -  29)  A  («("+*)  •  VF(n+0), iKn+0)) 

+  (1  -  29)  A  (u(n+0)  •  VF(n+0),5u(n+0)  •  Vf("+9))  +  (1  -  29)  A  (w(n+0)  •  vr(n+0),  F(n+0)) 

+  (1  -  29)  A  (u(n+0)  •  Vr(n+0),d'w(n+0)  •  VF<T+0))  .  (5.105) 

Bounding  each  piece,  using  u | X}  =  0,  and  V  •  u^n+^  =  0 
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Thus, 


(1  -  29)  A  (w(n+e)  •  =0,  (5.106) 

(1  -  29)  A  (w(n+fr>  •  VF(n+0~) ,  8u(n+ ®)  •  VF(,1+^  =  <5  (1  -  29)  A  u(n+^)  •  Vf(n+e')  2  ,  (5.107) 


(1  -  29)  A  (u(n+0~)  •  vr(n+0'),  F 


(1  —  29)  XM2d  II  jp(n+e)  II2  +  (1 


-  20)  A 


(5.108) 


(1  -  29)  A  (u(n+0~)  •  Vr(n+0),dM(n+0)  •  VF(,1+9)) 

<  <52  (1  -  20)  Aei  L(n+0~)  •  VF(n+0')  f  +  (1  ~  2f )  AM^ 


(5.109) 


The  next  convective  term  is  split  as 

9\  fu™  •  Ve£},  F^J"^  =  9X  (u^  •  \/F(n\F^n+§A  +  9X  ( u (n>  •  VF (n>,  <Wn)  •  VF 


and  bounded  by 


+  0A  f M(n)  •  vr(n),  F(re+e')")  +9x( u{n)  •  vr(n),  5u(n)  •  VF 


2  +  0A  VF(n)  2, 


9X  (u^  •  VFW,Wn)  •  VF(n+e~))  <  SeXM2d  VF(n+e)  2  +  89XM2d  VF(n)  2 


0A  •  VF^.fM)  <  dXA^2(*  F 


9X  (u (n)  •  vr(n\  F(n+0'))  <  F 


2  +  0A  vr^  2, 


(5.110) 


(5.111) 


(5.112) 


(5.113) 


9X  (u{n)  -VF  {n),Wn)  •  VF 


9XM2d 


2  529XM2d  2 

h - vr(n) 

2 


(5.114) 


0A  fw(n-0)  •  Ve^“0),  F^J)  =  0A  fu^n~^  •  VF^n~e\  F 


+  9X  (u^-v  •  vf^-^^uM)  .  vf(,1+0) )  +  0A  (yn~0)  •  vr(n-0),  f( 

+  0A  fu(n“0)  •  Vr(n-0),<5u(n-0)  •  VF 


with  the  bounds 


9X  (u{n~e)  •  VF(n“0),F 


9XM2d  „ 

<  - - -  F 

4 


2  +  0A  X7F(n~e)  2 


(5.115) 


(5.116) 
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0A  [ySn~e)  •  VF^n~e\5u^e)  •  Vf(w+i?))  <  SdX^  ^  VF(n+e)  2  +  6dXI^  d  yF{n-e)  2, 

,  2  (5-117) 

e\  [u{n~e)  •  vr(n~e),F(n+g~))  <  <>xxl^lt  F(n+d)  2  +  ex  vr(n-0)  2,  (5.118) 

and 

ex  fu(n~9)  •  vr(n"0),5ti(n“0)  -vfH'))  <  exx^  vf(”+9')  2  +  52exxi2^  vr(n_0)  2 . 

(5.119) 

For  the  first  ga{ •,  •)  term 

(1  -  26)  X  [ga{e^ ,  V^^)),  =  {l- 20)  X  (ga{F^n+§) , 

+(1  -  26)  X  (ga{F^n+9),  Vu(n+§)),8u(n+9)  •  VfH'))+(1  -  26)  X  (ga(T^n+§) , 

+  (1  -  20)  A  (ga(r(.n+e\\7u(n+e)),  5u(n+e)  •  VF(n+0))  ,  (5.120) 

with  each  part  bounded  by 

(1  -20)  A  (ga(F(n+9),Vu(n+9)),F(n+9)'j  <  (1  -20)  A4Md  F(n+0)  2,  (5.121) 

(1  -  20)  A  (ga{F(n+§),  Vu(n+§)),6u(n+§)  •  VF(n+0~)) 

<  4  (1  -  20)  Ad4M4  F(n+e)  2  +  52  (1  -  20)  A  VF(n+e)  2,  (5.122) 


(1  -  20)  A  ^a(r(n+f?),  VmH)),  f(”+0)) 

<4(1-  20)  XM2d2  r(n+e)  2  +  (1  -  20)  A  Ftn+d)  2  ,  (5.123) 

and 


(1  -  20)  A  (ga(T(n+§),  Vu(n+9)),8u(n+°)  •  VF(n+0~)) 

<4(1  —  20)  Ac/4M4  r(n+e)  2  +  52  (1  —  20)  A  VF(n+e~)  2.  (5.124) 

Next, 

(ga(e^\  Vt*(n>),F^)  =  0A  (Sa(FW,V«W),FH')) 

+  0A  (5a(F(n),Vu(n)),Wn)  •  VF(n+e))  +0A  ^a(r(n\Vu(n)),F(n+e)) 

+  0A^a(r(n),Vw(n)),<5u(n)  •  VF(n+0))  ,  (5.125) 

and  has  the  bounds 

0A(#a(F(n),Vu(n)),iKn+^)  <46XM2d2  F(n+e~)  2  +  0A  F(n)  2,  (5.126) 
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0A^o(F(n),Vt*(n)),<5w(n)-VF(n+®))  <40XM4d4  F(n)  2  +  529\  VF^n+§)  2,  (5.127) 

6X  (gfl(r(n),Vit(n)),F(n+e'))  <  4d\M2d2  F(n+e~)  2  +  9X  T(n)  2,  (5.128) 

and 

ex  (3fl(r(n),vuw),  5u(n)  •  vF(n+e~))  <  AexM4d4  r(n)  2  +  52ex  vf(,1+0_)  2 .  (5.129) 

The  last  term  in  this  grouping  is 

ex  {ga{e^\Vu^),F^}^  =  ex  (ga(F^~e\Vu^),F(n+9)) 

+  ex  (ga(F(n-e\Vu{n-e)),8u(n-0'>  •  VF(n+^))  +  ex  (ga(T{n-e\Vu(n-d)),  F(n+§)^) 

+  ex  (ga(r(n-0\  Vu(n~0)),  5u{n-e)  •  VF(n+0))  ,  (5.130) 

and  we  use  the  bounds 

ex(ga{Fin-d\Vuin-9)),F(n+9^  <4eXM2d2  F^n+e )  2  +  0A  F(n-0)  2,  (5.131) 

0A  ^a(F("-e),VM(n_0)),(5M(n-0)  •  VF(n+0)J  <46XM4d4  F^n~e)  2  +  826X  \7F(n+9)  2,  (5.132) 

#a  (c/a(r(n-e),vw(n“0)),F(n+e'))  <  AexM2d2  f(™+®)  2  +  0a  r(n~0)  2,  (5.133) 

and 

9A(fl0(rM)1ViiM)),iuM)  •  VF(n+®))  <40AAf4ci4  r(n"0)  2  +  52eX  VF(n+e~)  2.  (5.134) 


5.1.3  Bounding  B^(a,r)  —  B^(ah,r) 

In  this  sub  section  we  determine  bounds  on  the  terms  in  expression  (5.21),  B^(<t,t)  —  B^(ah,r), 
for  r  =  F(n+0).  Thus, 


;  (cr,  F('n+d'>)  -  B$(crh,  F(n+e ))  :=  A  (dte%+0) ,  F^)  +  (e, 


(n+6)  p(n+0) 
(7  J  *  5(n) 


I  a,  .  (An)  rp(n+ 0) 
+  ecu  \ea  ,  E  <j(n_9) 


I  /i  _  9m  ,  ( An+Q— 1)  p(n+0) 
+  (1  2fc'jo;lecr  >  *  50-0) 


1  2-  ( An~e)  j?{n+e) 
+  leCT  )  <5(n-0) 


+0A  («("“*)  •  Ve^},  F^-l)  +  («W  •  Ve£},  F^}) 

+0  A  (ga (e% ~6) ,  VU("-fl) ) ,  F^ )  +  eX  (ga (e£} 1  ,Vu<n)),F 

Each  of  the  terms  in  expression  (5.135)  are  now  bounded.  Here 


(5.135) 


A  (dt  (e^+0))  ,  F^n+e^  =  X  (dt  (V n+e ,  F(n+0))  +  X  (dt  (r(n+0>)  ,  F^n+9^  ,  (5.136) 


A  (dt  (F^n+9^  ,  F^n+0^  >  —  ^  p(n+e)  2  —  p(n+e)~l 


(5.137) 
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and 

A  (dt  (r(n+0))  ,F(n+0))  <  j  dt(Tin+e)^)  \\  F(n+0)  2.  (5.138) 

Next, 

9lo  (e^+0),  =  9co  ^F(n+0\  i^n+0))  +  9co  (f(-11+0\  6u^  •  VF(n+^) 

+  9co  (r(n+0),  F(n+e))  +  9co  (r(n+0),  5u(n)  •  VF(n+0))  .  (5.139) 

Bounding  each  term  we  have 

9co  (Vn+0),F(n+0))  =  9uj  F(n+0)  2,  (5.140) 

9co  (Vw+e),<5-u(n)  •  VF(n+g))  <  du}A*~d  F{n+e)  2  +  529uj  VF(n+(,)  2,  (5.141) 

9co  (r{n+9\  F(n+e^  <  ^  r(n+0)  2  +  y  F(n+0)  2,  (5.142) 

and 

9co  (r (n+e),  Ww)  •  VF(w+e))  <  r(n+0)  2  +  529lo  VF(n+0)  2.  (5.143) 

Again,  expanding  the  term 


9lo  (e(£\  Fg(n-6))  =  [F(n\F^n+e^  +  9co  {F(n\5u^n~^  •  \7F(n+e^ 

+  9co  (r(n),  F(n+e))  +  9co  (v{n\5u^n~e)  •  VF(n+0))  .  (5.144) 

Bounding  each  of  these 

9co  ^n\F(n+e^  <  ^  ||-F(n)f  +  Y  ||F(n+0)||2>  (5.145) 

9co  d  F^n)  2  +  629co  VF(n+^  2,  (5.146) 

0w  ^r(n),F(n+0))  <  ^  rW  2  +  y  ^(n+<?)  2>  (5.147) 

and 

0u;  (r(w),<fa(n~g)  •  VF(w+e))  <  9ujK^d  r(n)  2  +  629lu  VF(n+fl)  2.  (5.148) 

Next 


(1  -  29)  co  (e£+fl-1),  =  (1  -  20)  w  F(n+0)^ 

+  (1  -  20)  co  ^F(n+0“1),  Wn“0)  •  VF(n+0))  +  (1  -  20)  w  (r(n+0-1},  F(n+e)) 

+  (1  -  20)  w  (r^n+e~1\  5u{n~e)  •  VF(n+0))  ,  (5.149) 


and  has  bounds 


(1  -  20)  CO  ^(n+e-1),  F(n+e)^  <  C1  2°)  u  F(n+e- l)2+0_?p^  F(n+0)  2,  (5.159) 
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(1  -  29)  uj  ^(n+0-1),  5u{n-e)  •  \7F(n+e) ) 

^  (1  -29)uM2d 

~  A 


F{n+e-i)  2  +  s2  {1-29)  LO  VF(n+0) 


and 


(1  —  29)  u  ^r(n+61-1),  F(n+e^  <  — _ : p(n+e-i)  2  |  ^  29)  m 


j?(n+6) 


(1  -  29)  UJ  (r(n+®_1),  Su(n~0)  ■  VF(n+0)) 

^  (1  -  29)  uM2d 
~  4 

Splitting 

9Cj  (e^\F^e))  =  OCo  (F(n),FM) 


r(n+e-i)  2  +  s2{l-29)u  VF(n+0) 


+  9uj  (Vn),  5u{n)  •  VF(n+fl)) 

+  to  (r^, F(n+0))  +  to  (r^,  W")  •  VF(n+s^ 


we  have 


to  f F(n)  F(n+61^  <  —  F(n)  2  +  —  F(n+e) 
\  J  ~  2  2 


\  /  Z  II  II  Z  II  II 

to  ( F(n\5u W  •  VF(n+9))  <  euA*  ^  FW  2  +  52to  VF(n+0) 


and 


to  (r(n),F(re+,9))  <  ^ 
to  (r(n),<Wn)  •  VF(n+()))  <  6uJM2d 


|rH||2  _|_  ^  ||_p(«+0) 


The  term 

to  (egr^F^)  =  to 


r(n)  2  +  529uj  \7F(n+e) 


+  to  (F(n_0),<hi(n_0)  •  VF(n+0)) 

+  to  (r(n_0),  F(n+0))  +  to  (r(n_0),<Su(n_0)  •  VF(n+^ 


Here  each  term  is  treated  in  the  following  manner. 

to  f F^n~e^>  F(n+eA  <  —  F(n~d )  2  +  ^  _p(«+0) 

V  ’  /  -  2  2 

to  (Vn~e),fru(n~e)  •  VF(n+e))  <  ^  F(n_e)  2  +  e)2to  VF(n+0) 


and 


to  (r(n“0),F(n+0))  <  ^ 


|r'“-«»l|2  +  ? 


J?(n+d) 


Qu  {v{n-0\8u{n~e)  •  VF(n+e))  <  ^  r(n“0)  2  +  52to  VF(n+0) 


(5.151) 

(5.152) 


(5.153) 


(5.154) 

(5.155) 

(5.156) 

(5.157) 

(5.158) 


(5.159) 

(5.160) 

(5.161) 

(5.162) 

(5.163) 
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Examining  the  convective  terms, 


9X  ( u •  Ve^“0),  =  OX  [u(n~^  ■  yF{n-e\F^n+e^ 

+  ex  ( u (n_0)  •  vF(n_0),  5u(n~e)  •  vf(n+<?))  +  ex  (u{n~e)  •  vr(n_0),  F(n+e)) 

+  ex  (u^-0')  ■  yT^n-0\Su^n-^  •  VF(n+0))  .  (5.164) 

Here  each  part  is  bounded 

ex  F(n+0)  2 +  ex  2,  (5.165) 

ex  [u^n-G)  •  VF(n_0),  5u^~0)  •  VF(n+0))  <  S^X(^m2  vF(n“e)  2  +  5^m2  VF(n+e)  2  , 

(5.166) 

0A  (Vn~e)  •  vr(n~e),  F(n+e))  <  F[n+0)  2 +  ex  vr(n“e)  2,  (5.167) 

and 

ex  (u^-v  •  yT^n-e\5u^n-^  •  vF^w+e))  <  ^a^m2  vF(n+^  2+<y2^M"  vr(n-0)  2.  (5.168) 

The  second  convective  term  is: 

0A  (u(n)  •  Veff,  F$f)y)  =  ex  (u(n)  •  VF{n\  F(n+0))  +  0X  (u(n)  •  VF(n),  <5u(n)  •  VF(n+0)) 

+  0A  •  Vr^,  F(n+(9^  +  0A  (u^  •  VTW.JuW  •  VF(n+fl))  ,  (5.169) 

with  bounds  given  by 

0A  (Vn)  •  VF(n),F(n+e))  <  6X^2  F(n+e)  2  +  0A  VF(n)  2,  (5.170) 

0A(u<n>  ■  VF^\5u^  ■  VF^n+0^  <  5eX(^2  ||vF^  |2+  SO  MM2  |VF(n+e)  2,  (5171) 

ex  (rt(n)  •  vr(n),F(n+e))  <  eX(*xl2  i?’(»i+e)  2  +  0a  vr(n)  2,  (5.172) 

and 

ex(u^ .  yr(n\du^  •  vF(n+0))  <  ex ||vj,(n+e)||2+  s2ex^M~  ||vr^||2.  (5.173) 

Lastly  the  ga(-,-)  terms  are  considered.  Specifically 


ex  [ga(e^-9\Vu^),F{;(^  =  ex  ^ga(F^-0\Vu^),F^+e^ 

+  OX  [ga(F(n-e\yu^n-d)),5u^n-e)  •  VF(n+<?))  +  0A  ^a(r(n“0),  Vw(n"e)),  F(n+e)) 

+  ex  (ga(T{n-0\  Vw(n"0)),  5u^n~e)  •  VF(n+fl))  .  (5.174) 
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Handling  each  part 

0A(5a(F(n^,Vw(n_0)),F(n+0))  <A6XM2d2  F(n~e)  2  +  6\  F(n+e)  2,  (5.175) 

d\(ga{F(n-e\VvSn-e)),5u{n-6)  ■  VF^n+e^  <AOXM4d4  F(n~e)  2  +  526\  VF(n+0)  2,  (5.176) 

ex  (j0(rM),Vu(n“tf)),F(n+J))  <  46\M2d2  r(n-0)  2 +  ex  F(n+0)  2,  (5.177) 

and 

d\(ga{T^n-e\yu^n-e)),5u{n-e)  ■  VF^n+6^  <AOXM4d4  r(n-0)  2  +  526\  VF(n+0)  2.  (5.178) 
For  the  last  term: 

9X  (ga(eff,Vu^),F^d))  =  6X  ^(fW.VwW)^^) 

+  9X  (ga(F(n\  Vu(n)),<$i*(n)  •  VF(n+fl)) 

+  0A  ^a(r(n),Vw(n)),F(n+0))  +6X  ^ga(T^n\VuM),6u^  •  VF(n+0))  ,  (5.179) 
and  each  term  is  bounded  by 

6X  ^a(F(n),Vw(n)),F(re+e))  <  A9XM2d2  F(n)  2  +  #A  F(n+e)  2,  (5.180) 

ex  (ga{Fi-n\Vu^),5u^  •  VF(n+^  <  40XM4d4  F^n+e)  2  +  <520A  VF(n+0)  2,  (5.181) 

ex  ^a(r(n),vw(n)),F(re+e))  <  AexM2d2  F{n+e)  2 +  ex  r(n)  2,  (5.182) 

and 

ex  (ga{V^n\yu^),5u^  •  VF(n+e))  <  A6XM4d4  T(n)  2  +  629X  VF(,i+0)  2  .  (5.183) 

5.1.4  Bound  the  I\l(Fl'n+l))  Terms  in  (5.7) 

The  bounds  on  the  Z”(r)  term  are  now  found  with  r  =  F^n+1\  This  term  is  broken  into  two 
pieces  in  the  following  manner: 

2^(F(n+i))  =  Intaig(F^)  +  Intaiup(F (n+1)), 

with 

Inter lg(Ftn+^)  =  A  ( '  dt <r(n+1)  -  F(n+1)^  +  w  (Wn+e)  +  e^n+l^  -  a inH) ,  F^n+1^ 

+  cD  (Wn)  +  ea(n+§)  -  o-Ks),  F^n+1^ 

+  A  (W(n)  •  Vo-(n)  +  0u(n+e)  •  Vcr(n+0)  -  ■u(n+^)  .  V<r(n+^ ,  F(ri+1)) 

+  A  (%(ff(B),VuW)  +^(<tH“),V«H"))  -^((r^^V^+^.F^1)) 

+  2a  (d(u(n+5))  -  0d(u(n))  -  M(«M),F(,,+1>)  , 
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and 


Into- lnp(-F(n+1))  =  A 6  -VF(n+1)j 

+\e  ^-o-fn+^,<y«(n+®)  •  VF(n+1)j 

+u6  (at n+ 9)  -  fj(n+i), 5w(n)  •  VF(n+1)) 

+ujQ  (Vn+1)  -  cr(n+l),5w(n+0)  •  VF(n+1)) 

+w  (1  -  20)  ( '(T^+e ')  -  cr(n+5),  •  VF(n+1)) 

+u>0  (Vn)  -  <r(n+^),(Wn)  •  VF(n+1)) 

+u>0  (An+d)  -  a(n+5),Su(n+^  •  VF(n+1)) 

+0A  •  Vcr(n)  -  ii(n+  s)  •  Vcr(n+5),,5w(n)  •  VF(n+1)) 

+0A  (u(n+0')  •  Vcr(n+0)  -  w(n+5)  •  Vo(n+^,5u(n+e)  •  VF(n+1)) 

+0A  (ffa(<r(n),  VuW)  -  Vu(n+5)),  JuW  •  VF^1)) 

+0A  (Sa(<r(^,V«(^) •  VF^+1)) 

+02a  (d(u(n+5))  -  d (u(n)),<Wn)  •  VF(n+1)) 

+02a  (d(t>+^))  -  d(-a(n+^)),5w(n+0')  •  VF(n+1))  . 

Bounding  each  of  the  Intolg(F^n+1^)  terms: 

A  (dtcr(n+1^  -  <r|”+^,F(n+1)^  <  ^  dtcr{n+l)  -  o-["+^  +A  F^n+1)  2  ,  (5.184) 

u  (t 3a(n+e )  +  0cr(n+1)  -  o-(n+^) ,  F(n+1))  <  |  0cr(r!+e)  +  0o>+1)  -  <x(n+i)  2  (5.185) 

+  cj  F(n+1)  2,  (5.186) 

w  (Wn)  +  0cr(n+e')  -cr(n+^),F(n+1))  <  |  j|0cr(n)  +  0<r(n+e~)  -  <r(n+^  ^ +oj  ||.F(n+1) ||2  ,  (5.187) 

A  (W{n)  ■  Va(n)  +  6u(n+e )  •  Vcr(n+e)  -  «("+§)  •  V<r(n+^ ,  F(n+1)) 

<  ^  0w(n)  •  Vcr(n)  +  6u(n+e)  •  Vo-(n+0)  -  w(n+i)  •  Vcr(n+5)  +  A  "  ,  (5.188) 

A  (09a(^),VwW)  +  05a(o-(^),V>^))-5a(fr(n+i),Vw(”+^)),F^+1)) 

<  ^  ||05a(^),V«W)  +  05a(>+^,VW(n+0'))  -5a(<r(n+3))Vl*(n+3))|2 

+  A  ||.F(n+1)||2  ,  (5.189) 
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and 


2a  (d (w(n+5))  -  0d («("))  -  0d (w(n+0')),F(n+1) 

|d(w(n+3))  -  0d(u^)  -  0d(w(n+e')) 


<  - 

“  2 


+  2a 


_p(«+i) 


Bounding  each  of  the  Inter iup(F^n+1^)  terms 


A 0  ^~-a\n+^  ,5u{n)  -VF(n+1)j 


< 


A 982M2d 


X9  -a 


t  ? 


.  VF(n+1)  1  < 


4 

A  9S2M2d 


—a 


(n+i) 


—a 


(n+5) 


+  xe 


+  a  e 


^jp(n+l) 


V  F^n+1^ 


lu6  (Vn+0)  -  <r(n+!) ,  <$«(">  •  VF(-n+1^  < 


u682M2d  II  (n+e)  _cr(n+i)||2 


+  u;0 


\7  _p(ri+1) 


cu0  (tr(n+1)  -  a(n+^,5u(n+e)  •  VF(r!+1))  < 


u982M2d 


r(n+l)  _  rr(™+|) 


+  u>9 


^-p(n+l) 


cr  '  ’  —  <jy 

2 


5 


w  (1  -  20)  (Vn+0)  -  cr(n+l) ,  <^(n+e~)  •  VF(n+1))  <  ^ - 20)  5-M~d 


g.in+8)  _  c 


L o9  (Vn)  -  <r(n+3),Wn)  •  VF(n+1))  < 


uj052M2(1 


+  u(  1-  20) 

o»  -  cr(n+5) 


Y_p(«+1) 


uj6  ^a(n+e)  -  o-(n+i) ,  5u(n+0)  •  VF(n+1^  < 


u982M2d 

4 


2 

+  UJ0 

Y^(n+l) 

cr(n+0’) 

-  cr(n+5) 

+  u)0 


0A  (w(n)  •  V<x(n)  -  w(n+i)  •  \7a(n+^,5u(n)  •  VF(n+1)) 


< 


9X52M2d 


u 


(")  •  Vcr(n)  -  w(re+l)  •  Vcr(n+^)  2  +  0A  VF(n+1) 


0A  («("+*)  •  V<r(n+f?)  -  «("+§)  •  Vcr(n+I),<5w(n+0)  •  VF(n+1)) 


< 


9X52M2d 


u 


(n+0)  •  Vcr(n+^)  -  «("+§)  •  Vcr("+i)  2  +  0A  VF(n+1) 


2.  (5.190) 

(5.191) 

(5.192) 

(5.193) 

2 

(5.194) 

r(*4)  2 

(5.195) 

2,  (5.196) 
2 

(5.197) 

2,  (5.198) 

|2,  (5.199) 
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ex  (3a(<rW,V«W)  -  ga(An+1*) ,Vu(n+1*)),5uW  ■  VF("+1)) 

<  0X52^2d  \\ga(^n\  VuW)  -  5a(^(n+^ ,  V«(n+3))  ||2  +  0A  ||vf(”+1)  ||2  ,  (5.200) 

0~A  (9a(>+*),  Vt>+*))  -  <7a(>+^),  •  VF^+1)) 

<  0X52 M20  ||ga(<r(w+»),  Ww+e~))  -  ga(a(w+i),  Vn(w+^))|2  +  6\  |VJ,(W+1)|2  ,  (5.201) 

02a  (d(u(n+^)  -  d (u(n)),<kt(n)  •  VF(n+1)) 

<  0a82M2d  |d(M("+5))  -  d(M^)|2  +  92a  ||VF(W+1)||2  ,  (5.202) 

02a  (d («(n+5))  -  d(w(n+0')),<5«(n+0')  •  VF(n+1)) 

<  0a6^d  ||d(u(n+D)  -  d(-u(n+e'))||2  +  02a  ||VF(W+1)||2  .  (5.203) 

5.1.5  Bounding  the  Zr?(F(n+®))  Terms  in  (5.13) 

Here  we  bound  the  terms  in  Z^F^”-1-®))  in  expression  (5.13)  in  two  phases,  first  handling  the 

standard  Galerkin  terms,  and  then  the  upwinded  terms.  Consider 

In^F(n+8))  =  Z”9(z(n+0'))  +  T'lwp(F(n+d)),  (5.204) 

where 

l2g{F(n+°))  =  X  (dt*(n^  -  <T^n+*~6\  F(”+S)'\ 

+  lu  (ecr^+e)  +  (1  _  20)  +  da(n)  _  a(n+±-e) ^  F(n+e)^ 

+  u  ((1  -  20)  o-(n+*)  +  6a^  +  -  <r(n+M ,  f("+®)) 

+  A^  (1  —  20)  u(n+e')  •  Vcr(n+0')  +  9u^n)  •  Vo-(n)  +  0u(n"e)  •  Va^n~9) 

-  vin+^~d)  •  Vtr(n+5~e),F(n+^ 

+  a((1  -  20)  ga(<r(n+§) ,  Vu(n+9))  +  9ga(^n\Vu^) 

+  dga{(T^-d\  Vu(n-0))  -5a(o-(n+^0),Vu(n+^0)),F(n+0')^ 

+  2a^d(u(n+5-0))  -  0d(u(n))  -  (1  -  20)  d{u(n+9)) 
-ed(u^~^),F(n+S)\ 


35 


and 


Z2  uP(F^n+^)  =  X9  (-0tn+®-*\<$u(n)  •  VF(n+*)) 

+  A(1  -  20)  (^~-cr\n+^~e\5u^n+^-e)  •  Vf(n^j 

+  A 6  ^_4n+^Wn~0)  •  VF(n+0')j  +  6io  (a^n+e)  -  <x(n+f-0) ,  8u^  •  VF(,i+^ 

+  (1  -  26)u(a{n+e)  -  cr(n+5-e),8u(n+°)  •  VF(,1+9)) 

+  0u ;  (Vn)  -  cr ("+3-*)  ,<*«("-*)  • 

+  (1  -  26)  Co  (cr(n+,f>  -  cr (n+^e) ,8u(n+°)  •  VF(n^)) 

+  6Cb  (Vn)  -  <r(n+3_e),<y«W  •  VF(n+0)) 

+  6oj  (Vn“e)  -  cr(n+^“e),5u(n“0)  •  VF(n+fl)) 

+  (1  -  20)  A^«(n+0~)  •  Vcr(n+®)  -  u(n+  s-0)  •  Vo-(n+^-e),5w(n+0')  •  VF(n+0')j 
+  0A  (Vn)  •  Vcr(n)  -  w(n+5-0)  •  Vo-(n+5“0),  5w(n)  •  VF(,1+9)) 

+  0A  (Vn“e)  •  V<r(n“e)  -  u(n+^~e)  •  Vcr(n+i-0),5w(n“e)  •  VF(n+e)) 

+  (1  -20)A^a(er("+^),V«(n+*))  -  ga(a(n+5-e) ,Vu(n+t-°)),8u(n+°)  •  VF(n+^ 

+  0A  ^a(cr(n),VM(n))  -^a(fr(n+5-0),Vw(n+5-e)),<5M(n)  •  ViKn+^)) 

+  6X  (ga{a^n~d\  Vu{n~6))  -5a(cr(n+5-0),Vu(n+^e)),dM(n"0)  •  VF(n+f^) 

+  02a  (d(w(n+i-0))  -  d(u^),(5w(n)  •  Vf(n+9')) 

+  (1  -  20)  2a  (d(?>+5-0))  -  d(«(n+e~)),<5u(n+0~)  •  Vf(n+®)) 

+  02a  (d(w(n+5-0))  -  d (u^-e'>),8u{n-e'>  •  VF(n+0’))  . 

Bounding  the  2^ g(F(n+^)  terms: 

A  (^cr(n+0~)  -  tr^~0\  F(n+°)'\  <  ^  dta(n+§)  -  a[n+^  +  A  || F^n+°)  ||2  ,  (5.205) 

u  (da^n+e)  +  (1  -  20)  cr(n+0)  +  6a (n>  -  a(n+l2~e) ,  F(n+0')) 

<  j  ||0cr(n+f?)  +  (1  -  20)  o>+0)  +  0<r(">  -  cr(n+M  ||2  +  a;  ||F(n+0~)  ||2  ,  (5.206) 

w  ((1  -  20)  cr(n+^)  +  dtrW  +  6a{n~V  _  cr(»+§-e),  F(n+^)) 

<  ^  |(1  -  20)  fr(n+0')  +  0o-(n)  +  0o>-0)  -  cr(n+5-e)  ||2  +  £  ||F(n+<9~)  ||2  ,  (5.207) 
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A  (  (1  -  20)  u(n+d )  •  Va(n+e)  +  0u<n)  •  V<r(n)  +  0-u(n-0)  ■  -  u(n+|-0) .  Ver(n+3-0) ,  F(n+0)^ 

<  ^  |(1  -  20)  u(n+0')  •  V<r(n+0~)  +  0i/n)  •  +  Qu^-^  •  Vcr(n-0)  -  u (n+3“e)  •  Vcr(n+3-e)  ||2 

+  A  F(n+e)  2,  (5.208) 


A (  (1  -26)ga(An+9),Vu(n+d))  +  6ga(^n\Vu^) 


6ga{a^-e\  V«(n-fl))  -  ga(*(n+*-e) ,  Wn+^)),  F(n+0~) 


-4  (!  -  l^daicr 


+  6ga(a^,Vu^) 


+  6ga{(r(n-9)yu(n-e))  -  ga{a(n+^-e\\7u(n+^-e))  +A  f( 


(5.209) 


2a^d(-u(n+3-0))  -  0d(-i/n))  -  (1  -20)d(Jri+0'))  -  6d{u(n~V),F(n+°)'\ 


<  |  j|d(w(n+5-0))  -  0d(t*W)  -  (1  -  20)  d(w(n+e'))  -  0d(w^ 


+  2a  F(n+e)  .  (5.210) 


Bounding  the  terms  in  Z^  (F 


A  0  - 


!^(n+i  g)^M(n)  .  Vip(n+g)^  <  X6S2M2d  _<T(«+I  e)  +A0 


(5.211) 


A(1  -20)  (^-crJn+2  0^w(n+i-0)  •  VF(n+^ 


au-m)^  2  +  a(1-20) ||vf(? 


^_^(n+i  e)^M(w-g)  ,VF(n+e)"j  <  A eS-M-d  _<j.(n+  {  «)  +  A0  ||vF(n+e~)||\  (5.213) 


A0  - 


(5.212) 


0w  (er(n+0)  -  cr(n+5-0),5u(n)  •  VF 


6io52M2d  11  aBI2 


r(n+<9)  -  o-(n+2-°)  +0w 


(5.214) 


(1  -  20)  u  (<7(n+0)  -  a(n+^-e\du(n+§)  •  VF 


(1  Wu52  ||CT(n+0)  _  <T("+|-e)||2  +  (1  _  20)  cuei  |M(n+0')  •  VF(n+0')||2  ,  (5.215) 
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0lu  -  a^-0),  6u(n~V  •  Vf(n+®)) 


<  °^2M2d  a(n)  a(n+i-e) 

~  4 

2  +  0a;  VF(n+0')  2, 

(5.216) 

(1  -  20)  u>  (cr(n+0)  -  o-(n+^~0),Su(n+0)  •  VF(n+0')) 

<(l-20)w  ^(n+e)  ^(n+i-0)  \S2(1  20)  c5ei 

4ei 

u(n+0)  .  y p(n+e) 

(5.217) 

0a;  •  VfM) 

^  0&S2M2d  ^  CT(n+i_0) 

“  4 

2  +  0a;  VF(n+0~)  2, 

(5.218) 

0a;  (<r(n_0)  -  er(n+5-0) ,  Su{n~0)  •  Vf("+e)) 

^  9Co52M2d  ^n_e)  a(n+i-e) 
4 

2  +  0a;  VF(n+0~)  2, 

(5.219) 

(1  -  29)  \(u(n+0)  •  Vcr(n+0)  -  u(n+ M)  •  Va(n+^-0),5u(n+0)  •  VF(’1+0)j 

<  ^  ~  26^  A  w(n+0~)  •  Vcr(n+^)  -  u(n+5“0)  •  Vcr(n+5-0)  2 
4ei 

+  52  (1  -  29)  Aei  w(n+0)  •  VF(n+0)  2  ,  (5.220) 

9X  (u{n)  •  VcrW  _  u{n+h~e)  .  Vo-(n+5_0),d'u(n)  •  VF(n+<9)) 

<  U(V)  .  Vcr(n)  _  u{n+\-e)  .  v<T(n+|-e)  2  +  Qx  VF{n+0)  2  ?  (5-221) 

0A  •  V<T(n_0)  -  ■u(n+5-0)  .  Vcr(n+5-0),  <$u("-*)  •  VF(,,+9)) 

^9X6HPd  u(n-«).  Vo.(n-e)_w(n+i-e)  .  Vo.(n+i-e)  2  +  VF(n+0')  2,  (5.222) 

(1  -20)A^a(cr(n+0'),Vw(n+^))  -ffa(cr(n+5-e),Vw(n+^-0)),<5ti(n+0')  •  VF(n+^ 

<  ^~2g)A  5a(crH,VtiH')) -5a(fr(n+^-0),Vw(n+^-0))  2 

4ei 

+  52  (1  -  20)  Aei  ■u(n+e)  •  VF(n+0)  2  ,  (5.223) 
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0A  (ga{^n\Vu^)  -5a(cr(n+i-e),Vw(n+3-0)),Wn)  •  Vf(n+0')) 

<  dX62^2<*  jga(^n\\/u^)  -  ffa(cr(n+|-e),Vti(n+i-0))||2  +  0A||vF(n+0')||2,  (5.224) 

0A  (54<7<1,-e),Vu(n-0))  -  £/a(o-(n+^"e),  Vw(n+^“0)),<5ii(n"0)  •  VF(n+0')) 

<  9X62^2d  |ffa(cr(n-0),Vu(n-0))  -  ffa(cr(n+|-e),Vti(n+i-0))||2  +  0A||vF(n+0')||2,  (5.225) 

02a  (d(u(n+5_e))  -  d (tiW),Wn)  •  VF(n+^)) 

<  9a52^2(*  ||d(u(w+M)  _  d(n(n))|2  +  02 a  ||vF(n+e~)  | 2  ,  (5.226) 

(1  -  20)  2a  (d(w(n+5-0))  _  d(w(n+e~)),  Su(n+d)  •  Vf(n+9')) 

<  (1  ~2g)a  d(w(n+5-0))  _  d(w(n+e'))  2  +  52  (  1  -  20)  2aei  u(n+0~)  •  Vf(n+fl')  2  ,  (5.227) 

2ei 

and 

02a  (d(w("+5-0))  -  d(ii(n-0)),5u(n-0)  •  VF(n+^)) 

<  92a52M2d  ||d(n("+^-g))  -  d(«(w-*))|2  +  02a  ||vF(n+e~)  ||2  .  (5.228) 

5.1.6  Bounding  the  T^(F<'n+e'))  Terms  in  (5.20) 

Bounding  the  T^{Fl'n+^)  terms  in  (5.20)  is  done  in  two  phases,  first  the  standard  Galerkin  terms 

and  then  the  upwinded  terms. 

I. J(F(n+0))  =  T^g{F{n+e))  +  I!}up{F{n+e)),  (5.229) 

where 

2£fl(F(n+*>)  =  A  (dttr^+V  -  a-[n+e~^\F{-n+e^j 

+  u  (ocr^  +  0<r(n-e)  -  a(n+e~^\  F^n+e^ 

+  u  (e<T(n+V  +  0o-(n)  +  (1  -  20)  a(n+e~1)  -  ,  F^n+e^ 

+  A  (eu{n~e)  •  Vcr(n“0)  +  0u(n)  •  Vo-(n)  -  u(n+e_^)  •  Vo-(n+0“^),F(n+0^ 

+  A  (dga(tT(n-e\Vu^-^)  +0ga(an,Vu^)  -  ga{An+e-^  ,Vu(n+e~^),  F^n+e^ 

+  2a  (d(u(n+e~s))  -  0d(u^)  -  0d(u(n-^),F(n+0))  , 
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and 


l%up{F(n+e))  =  A0  ^-c t{U+ 6  ^\du(n)  •  VF(n+(,)j  +  A0  ^,5u^n-e)  •  VF(n+s)j 

+  Oco  (a(n+e)  -  <r(n+e-^),<5-u(n)  •  VF(n+0)) 

+  Olo  (Vn)  -  a(n+e-^,5u^-d)  •  VF(n+0)) 

+  (1  -  20)  u  ^(T(n+0“1)  -  cr(n+0-i) ,  5u(n~e)  •  VF(n+fl)) 

+  02  (Vn>  _  .  VF(n+^) 

+  02  (Vn-0)  -  fr(n+0“^),<5u(n-^  •  VF(n+0)) 

+  0A  (^M(n“0)  •  V(T(n_0)  -  u(n+e~^  •  Vo-(n+e-5),,5ii(n-0)  •  VF(n+0)) 

+  0X  («(")  •  Vtr(n)  -  u(n+e~ 2)  •  Vo-(n+e“i),(5it(ri)  •  VF(n+9)) 

+  OX  [ga{<T^n~e\yu^n-e))  -  5a(o-(n+0“^),  Wn+^)),<Wn-0)  •  VF(n+0)) 

+  0A  ^a(<rn,  Vit(n))  -  5a(o-("+0-^),  Vt*(n+®_3)),<Su(n)  .  VF(n+9)) 

+  02a  (d(u(n+0-^))  -  d(w(n)),<5u(n)  •  VF(n+0)) 

+  02a  (d(u(n+,H))  _  d («("“*>),  <Su(n-fl)  •  VF(n+0))  . 

Bounding  the  terms  in  2 3 g(F^n+^): 

X  (dta(n+e)  -  a[n+e~^\F^+e^j  <  j  dta{n+e)  -  cr\n+e~^  +  A  ||F(n+0)||2  ,  (5.230) 

2  ( 0(T ^  +  0(7^-^  -  a(n+e~^ ,  F^n+e^  <  j  9a (n)  +  9a(n~^  -  a(n+e~^  2  (5.231) 

+  2  F^n+^  2,  (5.232) 

u  (ea^n+d)  +  0a^  +  (1  -  20)  cr^-1)  -  (rK^Mj 

<  |  ||0cr(n+0)  +  0O-W  +  (1  -  20)  o-^-1)  -  cr(n+e4)||2  +  w  ||^^+0)||2  ,  (5.233) 

A  (Wn_e)  •  Ver(n-0)  +  0w(n)  •  Vcr(n)  -  u(n+6~^  •  Vo-(n+0_^),  F(n+f?)) 

<  ^  0w(n“0)  •  Vcr(n-0)  +  6u(n)  •  Vcr(n)  -  •u(n+,9-|)  •  Vcr(n+0-i)  2  +  A  F(n+e)  2  ,  (5.234) 

A(0<7a(o-(n_e),Vw(n-0))  +  05a(cr,l,Vw(n))  ~  ga(cr(n+e-^,Vu(n+d-^),F(-n+e^ 

<  j  \\ega(^n-e\Vu^)  +  9ga(vn,VuW)  -ffa(<T(n+e-^),V«(n+0-^)|2 

+  A||F(n+0)||2,  (5.235) 
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2 


2a  (d (u(n+0~^)  -  9d («("))  -  9d(u^-^),F^n+e^ 

|d(u(n+^))  -  9d{u ("))  -  0d(w(n“0)) 


+  2a 


j7>('«+0) 


Bounding  the  l^up(F<'n+e'>)  terms: 

xe 

A0  ^-4n+0”^,^(n_0)  •  VF(n+0A 


< 


A  952M2d 


—a 


{n+e~\ ) 


< 


xes2M2d 


—a 


(n+0-i) 


+  xe 


+  A  9 


\/jr(n+e) 


^p(n+e) 


9co  (Vn+0)  -  cr(n+e-^,6u^  •  VF(n+0)) 


< 


9uj52M2d 


r(n+e)  _  {n+e-l) 


+  9lo 


■^jp(n+e) 


9uj  (a(n)  -  cr("+0“5))(Jw(tl-0)  .  VF(n+0)) 


< 


9uj82M2(1 


o»  _  a(n+e-h) 


+  9uj 


\/jr(n+e) 


(1  -  2 9)u(a(n+0-V  -  An+e-^,8u^~d)  •  VF(n+e >) 


< 


(1  -  29)  u82M2d  ||^(n+(?_i)  __  (n+e_i)"2 


+  (1-2  9)u 


yj?(n+e) 


9Co  (<r(n)  -  (rH“s),5uW  •  VF(n+f?)) 


< 


9uj52M2(1 


o-(n)  -  cr(n+0“§) 


+ 


y^i(n+6>) 


(er(n-0)  -  cr(n+,?-i),  8u^n-e)  •  VF(n+0)) 


< 


9C852M2d 


Xn~9)  _  ^(n+e-l) 


+ 


y_p(n+6>) 


<9A  •  Vcr(n-0)  -  u(n+e~ s)  •  •  VF(n+fl)) 


< 


9X52M2d 


u 


(n-0)  •  Vcr(n_e)  -  w(n+0_D  •  Ver(n+0-5)  +  #A  VF(n+0) 


9X  (u(n)  •  V<r(n)  -  -u(n+e-5)  .  •  VF(n+0)) 


< 


9X52M2d 


u 


W  •  Vtr(n)  -  w(n+0_5)  .  V<r(n+0_5)  2  +  9\  VF(n+0) 


.  (5.236) 

(5.237) 

(5.238) 

,  (5.239) 

,  (5.240) 

,  (5-241) 

,  (5-242) 

,  (5.243) 

,  (5-244) 

,  (5.245) 
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9X  (5a(<7<B-e),Vu(n-0))  -  ga{(r(n+e-^ ,Vu(n+9-^),6u(n~e)  •  VF(n+0)) 


<  jga(^n-e\Vu^n-^)  -  ga{a(n+d-^,Vu(n+e-^))^  +  eX^VF^n+9)j2  ,  (5.246) 


ex  (ga(anyu(n))  -  ga(An+e-^ ,Vu(n+e-^),8u{n)  •  VF(n+0)) 


<  6X5^^  j9a(an,Vu^)- ga{a(n+e-^,Vu(n+e-^)j2  +  9X^VF^n+e^2  ,  (5.247) 


02a  (d(u(n+e_^))  -  d (tiW),Wn)  •  VF(n+0)) 


<  eaS~M~d  ||d(u(n+0-i))  _  d(MM)||2  +  92a  ||vF(n+0)||2  ,  (5.248) 


92a  (d(u(n+0-t))  -  d(u(n-^),5u(n-^  •  VF(n+<?)) 

<  |d(u(n+«-§))  _  d(u(n-^))||2  +  Q2a  ||vF(n+0)||2  .  (5.249) 

5.1.7  Bring  all  bounds  together: 

The  goal  is  to  now  bring  all  the  bounded  expressions  together  into  a  single  expression  so  that  the 
discrete  Gronwall  lemma  may  be  applied,  and  the  error  estimate  established.  Thus,  bringing  all  the 
bounds  on  the  three  linear  combinations  (4.16),  (4.17),  and  (4.18)  together  into  a  single  expression 
we  obtain 
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A 

2A t 
+9oo 


p(n+ 1)  _  jp(ri) 

]?(n+l)"2 


+ 


F(n+e) 


jr(n+d) 


< 


+  (1  -  29)  Co 
+  (1  -  29)  (8X  -  e\52  (2  +  3A  +  2a)) 

2( 


+  9oo 


pi(n-9)  "  _|_  jp (n+0)  _  jp(n+9)-l 

jp(n+9)"2 


U 


(n+9)  _  yjp(n+9) 


fW  ir  (  ex  +  9i0  +  ^  +  AOXMd  +  40  AM2  d2  +  80AM4d4  +  +  39uJ^2d 


+ 


+ 


+ 


11  V 

jp(n+l) 
jp(n+9-l ) 
jp(n+9) 


^l  +  4A  +  2a  +  0  +  0w  +  (l -20)w+ —  +  ^ 


9oo  9toM2d  XdM 2 


+ 


+  2  XdM 


) 


2  /  (1  —  2 9)  to 


1  + 


M2d 


4i4  ,  ,,mu2j2  ,  i1  20)W 


+ 

p(n-0) 

+ 

p(n+9) 

’( 

/  r.  .  _  ~  9uM2d  „  0/1,  0w  (1  -  20)  wM2d 

[  1  +  5A  +  0A  +  2a  (1  —  2$)  cu  - 

V  2  2  4 

9u>M2d  oo  XdM2  ,  ia 

+ - b  -  + - b  40AM  g?  +  49XM  d  + 

4  2  2  4ei  y 

-  +  0A  +  +  40AM2  J2  +  4AM4riA 

24  y 

1  +  4A  +  w  +  0  +  0(1  -  2oo)  +  (1  -  20)  A  +  2a  +  ^ 


+  — — 2g)  AM^  +  40 x(jM  +  9XM2d  +  4(1-  20)  A Md 

2  12  ,  (1  —  2^)  ^ 


+ 


+ 


+ 


+ 


yp{n-9) 

*( 

yF(n+9) 

2( 

\7jp(n+s) 

1 

\/p(n+l) 

*( 

+4  (1  -  20)  A M4J4  +  40AM4d4  +  160AM2  d2  + 
8XM2d 


X  + 


+  j20A^J  +  vf(b)  2  ^ 


3A0  1  + 


+)) 


(  2A  +  20  +  30A  +  4a0  +  +  9XA^  +  Q8XM2d  +  2d2  (A  +  20) 

.  ,  „  XdM2  8 XdM2 


+ 


1  +  3A  +  2a  + 


2  2 
8  XdM2  XdM2 

- r - b  — r— 


+  82  (uo  +  29  +  2X  +  (1  -  20)  oo  +  2 6u> 
+  82oo  (2  -  0)  +  28 2do  +  24'2A^ 


+  Hp  +  Hg  +  82Hupa  +  82Hupb- 

Here  the  terms  Hp,  Hg,  Hupa,  and  Hupb  in  (5.250)  are  defined  by 


(5.250) 
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^  dtv^-cr^  dtAn+§)  -  a^-9)  2 

+  ^  dt<j(n+V  -  *(n+d~^  +J  9a(n+e^  +  9a(n+l^  -  a(n+12)  2 
+  -  0o-(n+^  +  (1  -  29)  ain+^  +  0<jW  -  a(n+12~9)  2 
+  ^  |0<r("+«)  +  OffW  +  (i  _  29)  a(n+e~^  -  o-(n+e“5)  ||2 

+  -  Qa(n)  +  da(n+§)  -  <r(n+5)  2 
4 

+  j  || (i  _  29)  <r(n+e~)  +  0<rW  +  ^cr^-0)  -  <r(n+3-®)  ||2 

+  -  9a W  +  -  cr(n+0-i)  2 

4 

+  -  0MW  .  VcrW  +  0«(n+0)  •  Vcr(n+0)  -  W(n+§)  •  Vcr(n+§)  2 

4 

+  ^  ||  (1  -  20)  t*(n+®)  •  Vo-(n+e~)  +  0«W  •  VcrW 

+0u(n-9)  •  Vf7(n-e)  -  u(n+^-0)  .  Vcr(n+^-0)  2 

+  A  0tt(n-fl)  •  V<r(n-e)  +  0«W  .  VcrW  _  u(™+°-h)  .  Vcr(n+0_5)  2 
4 

+  ^  |^a(o-W,  V«W)  +  ^((rM.VuM)  -  9a(cr("+^),  V«(n+3))||2 
+  ^  (1  -  20)5a(ff(n+fl‘), V«(n+fl'))  l%(.W,  VuW) 

+  ^a(cr(n-0),Vw(n-0))  -ffa(cr(n+5-0),V-u(n+5-0)) 

+  ^  ||^a(o-(n-0),  Vu^)  +  0ga(an,  V«W)  -  ffa(CT(n+0-^) ,  Vu^s))  ||2 
+  |  ||d(Jn+5))  -  0d(t*W)  _  0d(«(n+®))|2 
+  |  |d(w(n+i-0))  -  0d(t*W)  -  (l  _  29)  d(u(n+®))  -  6»d(^W-e))||2 
+  |  ||d(w(n+0-3))  -  0d(t*W)  -  0d(^W-e))||2  ;  (5.251) 
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2 


H 


XM2d 


upA 


—a 


(n+§) 


2  XM2d 

+  — - — 


—a 


(n+i-0) 


2  XM2d 

+  — - — 


—a 


(n+e-\) 


+- 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


u9M2d 


r(n+0)  _  (n+|) 


Lo9M2d 


r(«+!)  _  o-(«+§) 


w(l  -29)M2d 


9u>M2d 

4 

u9M2d 


M+°)  -  (r{n+h) 


(n+e)  _  (n+i-0) 

2  w0Af2d 

U  U  '  ^  7 

+  4 

r(n+§)  _  fr(n+  5) 


2  (1  -2  6)u 


(1  -  29)  Cj 
4ei 
9ujM2d 


r(n+°)  _  rrC'^+l-6 


4ei 


2  9ojM2(1 

+  - : - 


M+e)  _  (n+i-e) 


rW  _ 


4 

9ujM2d 


O-W  _  o-O^+iH) 

cr(n+0)  -  cr(n+e_5) 


9uM2d 


+ 


4 

9loM2(1 


Xn~e)  _  cr(»+5-0) 
r(")  -  rr(n+0"5) 


(1  -29)ujM2d 


r(n+9- 1)  _  (n+e-|) 


+ 


9uM2d 


>)  _  rT(n+0-3; 


9uM2d  W^n-9)  _  cr(ri+e_i)|'2 
4 

9XM2d 


4 

9XM2d 


u 


M  .  v<r(n)  -  u(n+§)  •  Vcr(n+5) 


u 


in+'e)  .  Vcr(n+0)  -  u(n+!)  •  Vcr("+i) 


(1  -20)  A 
4ei 
9XM2d 


u 


in+°)  .  Vcr(n+^)  -  w(n+3-0)  .  Vcr(n+^“e) 


4 

9XM2d 

4 

9XM2d 

4 

9XM2d 


u 


M  .  Vcr(n)  -  w(n+3-0)  .  V<x(n+iH 


(n_0)  •  V<T(n_0)  -  M(n+iH)  •  Vcr("+i-0) 
(n_0)  •  V<T(n_0)  -  u(n+ e~^  •  Vcr(n+e_5) 


u 


W  .  Vcr(n)  _  u{n+0~l)  .  V<x(n+e_^) 


(5.252) 


45 


H 


e\M2d 


upB 


(1  -20)  A 
4ei 
0XM2d 
4 

e\M2d 
4 

ex  M2d 
4 

exM2d 


4 

0aM2d 

2 

0aAf2d 


,  V>+0~) )  -  Sa((7(n+s) ,  V«(n+*)) 
5a(o-(n+®),  VuM)  - 

ga(<r(n\VuW)  -  ga(a(.n+12-e)  ,Vu(n+12~d))  2 
ga{o-(n~e) ,  Vn^-"))  -  ga{cr(n+^-e) ,  Vu(n+2~e) ) 

ga{a(n~e\Vu(n-0))  -  ga{a(n+e-2) ,Vu(n+e-2)) 

ga(*n,VuW) -ga(a(n+e-2),X7u(n+e-2))  " 


d(u(n^))-d(wW) 


2  0aM2d 

+  — r - 


d(«(n+3))-d(Jn+0')) 


d(u(n+M)  -d(u^) 


2  (1  -20)  a 


26! 


d(u(n+s-®))  -d(u(n+e~)) 


e2aM2i  ||d(u(”+l-o))-d(u<”-'») 

d{U(n+9-2))  -  d{u^n~^) 


4 

0aM2d 


2  0aM2d 
H - ^ - 


d(ti(n+0-^))  -d(i*W) 


(5.253) 


and 

*  -  * 

+ 


+ 


dt  ( r(n+0) 


+  4 


+ 


+ 


p(n-0) 

p(n+0) 

p(n+l) 

p(n+d) 


dt  (r(n+1) 
u>  uM2d 


+  4 


dt  r 


0A  +  ^  +  40XM2d 2  +  4AM4J4 


2  /'w  0cu  0(jjM2d  (1  -  20)  w  (1  -  20)  w  (1  -  20)  wA^d^ 

1  2  +  T+3^^  +  ^^  +  2  + 


2  ^0u;^  ^ 


p(n+0—  1) 


2/(l-20)w  (1-2  0)uM2d\ 

2  +  4 


2  ^0c5  (1  —  20)  a>  BuM2d  (1  -  20)  w 

2~  +  2  +  4  +  4ei 


+40AM0  +  4  (1  -  20)  XM2d 2  +  40A M4d4  +  4(1-  20)  AM4d4 

MA  +  fc  +  3*4  +  iSUtd  +  2+llM  +  3  A+l  +  wxmW 

2  4  4 


+ 


+ 


+ 


p(n) 


yp(n-0) 


vr(n+0) 


'  S2XM2d\ 

°x+—) 


+ 


vr(n) 


30A  +  3 


520XdM2\ 


0A  + (1-20)  A +(1-30) 


S2XM2d  t  (1-20)  A M2d  20S2XM2d\ 

4+  +  2  J 


+ 


Multiplying  (5.250)  by  2At/X  and  summing  from  an  initial  time  t  =  0  to  time  t  =  l  a  telescoping 
of  terms  is  achieved.  An  inverse  inequality  is  also  applied  to  the  ||V.F||  terms  on  the  RHS  of  the 
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expression  (5.250).  Bounds  on  ||.Fe||J  and  [[i7,1  e||2  (the  terms  that  do  not  telescope)  are  stated 
and  proved  in  the  appendix.  Thus, 


pZ  — 1 


_ if 


TlZ  +  l  —  6 


:il  —  6 


?l+0 


+At  ^  ^  — - Ouj 


n—0 


p(n+ 1) 


n— 0 


n(n+6>) 


At  —0a; 


+AfE  T  (1  -  20)  (6X  -  e±52  (2  +  3A  +  2a))  L(n+®)  •  VF(n+0 


n—0 

l 


n—0 
2 


r»  (n+0) 


<  At^  F(n)  2  y  (  20A  +  20a;  +  30a;  +  89XMd  +  80AM2  d2  +  160AM4  d4  +  9u>M2d  + 


n—0 


{  3  0a)M2d 


+ch~ 2  (exe  (2  + mm2)) 


A*E  F 


(n+l) 


n—0 


2  \  (  2  +  8A  +  4a  +  29  +  26Cb  +  2  (1  -  29)  u  +  0w  + 
A 


6u>M2d 


XdM2  +  4A  dM 


+  Ch~ 2  2  +  6A  +  4a  +  SXdM2  +  XdM 2  +  2 S2io  (2  -  9)  +  \d2u  +  4<52A 


+  A  tJ2  F 


{n+6- 1) 


n—0 

l 


2 1 
A 


(1  -  29)  u>2  1  + 


M2d 


+  At^  F 


( n+8 ) 


n—0 


2 1 


A  ^2  +  10A  +  2 9u  +  2u  +  29X  +  4a  +  9uM2d  +  2  (1  -  2 9)lo  +  9uj 

(1  -  29)  u>M2d  9u>M2d  A  -  ,2  ;2  (1  -  20)  w 

+  2 - A - +  — - —  +  u  +  XdM2  +  80AM4  d4  +  80AM2  d2  +  - — - — — 

2  2  26i 

+Ch~ 2  ^2  +  6A  +  4a  +  XdM2  +  6XdM2  +  2 S2  (u  +  29  +  2X  +  (1  -  20)  oj  +  2 0w)  j  j 


A*E  F 


(n-S) 


n—0 


/ 


Y  (  ^  +  20A  + 
A 


£>M2d 


80AM2  d2  +  8AM4  d4 


+  C/i"2  2A  +  5XM2d  +  2S29X 


+  AiE  |F(n+^ 

n—0 

+ 


2i/ 

-  (  2  +  8A  +  2w  +  20  +  20(1  -  2u)  +  2  (1  -  20)  A  +  4a  +  9u>  +  (1  -  20)  Q 
9u>M2d  (1  -  20)  A M2d 


89XdM  +  29XM2d  +8(1-  20)  A Md 


2  2 
r4  J4 


4 44  ,  ooa\n.r2j2  ,  (1  2^)  ^ 


+8  (1  -  20)  AM4d4  +  80AM4  d4  +  320AMM 


2ei 


+  Cd"2  ^4A  +  40  +  60A  +  8a0  +  A M2d  +  0AM2d  +  29SXM2d  +  AS2  (A  +  20) 
*  2  1  2  /  2  z  2 

At  E  IHr  +  Ai  E  +  Ai  E  +  At  ^  j62HupB. 


T  2j 


f2  J 


(5.254) 


n—0  n=0  n—0  n—0 

To  allow  for  the  application  of  the  discrete  Gronwall’s  Lemma  the  inequalities  stated  in  (5.255), 
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(5.256),  and  (5.257)  must  be  satisfied. 


1  >  At- 1  29X  +  29uj  +  39ti  +  89XMd  +  89XM2d2 +  169XM'ld* +  9uMzd  + 

A  y  Z 

o  ^ 

2  +  8A  +  4a  +  29  +  29Cj  +  2  (1  -  29)  oj  +  9lo  +  — +  XdM2  +  4A dM 
+  ChT 2  (  6A0  (  2  +  )  +  2  +  6A  +  4a  +  <5AdM 


r4  j4 


3  9u;M2d 


+XdM 2  +  2 52lu  (2  -9)  +  452oj  +  4<52A 


(5.255) 


1  >  At  A  ^2  +  2(1  —  2d)  uj  (l  +  +  10A  +  29uj  +  2u)  +  20A  +  4a  +  9uM2d  +  2(1  —  26)  iv 


~  (1  -  2d)  uM2d  9u)M2d 

+9uj  +  - - i - + - - - 


T2j2  ,  (1  26)  LO 


+  u  +  XdM 1  +  89X APcP  +  89XM  d  + 

2  2  2ei 

+ChT 2  (  2  +  6A  +  4a  +  XdM 2  +  8  XdM2  +  2  82  (lo +  29 +  2X+ (1  -  29)  uj  +  2  9lu)  )  )  (5.256) 


u>M2d 


1  >  At f  (  2  +  u  +  29X  +  +  89XM2d2  +  8X M4d4 


+8A  +  2lo  +  29  +  29(1  -  2w)  +  2  (1  -  29)  X  +  4a  +  9Cb  +  (1  -  29)  to 
+  fMM2d  +  (1  -  29)  XM2d  +  8yX(jM  +  20XM2d  +  8(1-  29)  XMd 


+8  (1  -  29)  XM4d4  +  89XM4d4  +  320AM2  d2  + 

+  Ch~ 2  (2X  +  8XM2d  +  2529X  +  4A  +  49  +  69X  +  8a9 


4  44 


f  2  42  1  (1  29)  Cj 


2ei 


+XM2d  +  9XM2d  +  29SXM2d  +  4<S2  (A  +  29) 


r2j 


2,1 


(5.257) 


Summarizing  expressions  (5.255),  (5.256),  and  (5.257)  gives  a  stability  result  needed  for  the  discrete 
Gronwall’s  Lemma  of  the  form: 


At  (C  +  Ch~ 2  +  SC  +  5Ch~2  +  S2Ch~ 2)  <  1. 
This  leads  to  the  computationally  restrictive  stability  relationship 

At  <  Ch 2 . 


(5.258) 


(5.259) 


At  this  junction  in  order  to  complete  the  proof  of  the  a  priori  error  estimate  we  need  to  consider 
the  bounds  on  each  of  the  Hp,  Hg ,  Hupa ,  and  Hupb  terms  in  expression  (5.254)  using  interpolation 
properties,  and  some  of  the  lemmas  from  the  appendix. 
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5.1.8  Bounding  Hp 


Using  interpolation  properties  of  the  approximation  space,  and  elements  of  order  m  for  the  stress 
we  approximate  the  terms  in  At  Yln=o  \  Hp-  With  C  and  C)  denoting  constants  independent  of 
the  discretization  and  upwinding  parameters  h,  At,  and  5,  exemplary  terms  are  handled  as, 


Thus, 


Et  Ur(”+I) 


n=0 


E 


At 


n=0 

l 

£  E 

n=  0 
l 


i  ++1  dr  , 

At  Jt  ~dtdt 


\  [■  ftn+l 

2A t 


l2  dt 


tn+l  /gp' 

dt  dA 
at 


E; 

n= 0 


tn  +  l  /dY\2 

—  dtdA 


<  Ch2m+2  11  —  1,2 


m\ 


0,m+l  ’ 


(5.260) 


(5.261) 


^At(Ci)  T(n+0)  <  (Ch2m+2)  ]T  At 


n= 0 


M+e) 


<  Ch2m+ 2 


\n= 0 


l0,m+l  > 


m+1 


(5.262) 


J2 +  <52C3)  Vr^  <  (Ch2m  +  5Ch2m)  E  At 


n= 0 


\n=0 


rM 


m+1 


<  Ch2m  |  Ill'll lo,m+i  +  52Ch2m  ||Mllo,m+i  i 


(5.263) 


At^2jHp  <  C(h2m+2  +  h2m  +  62h2m)\\\cr\ 


1 2 

lo,m+l 


n=0 


+  Cft2”+2  |k,||*ra+1 . 


(5.264) 


5.2  Bounding  h9 

For  each  of  the  terms  in  Hg  a  typical  term  is  bounded  here.  To  complete  the  bounds  we  use  Taylor 
series  expansions.  The  first  bound  is  shown  in  detail  here,  and  then  the  reader  is  referred  to  a  series 
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of  useful  lemmas  given  in  the  appendix  that  aid  in  bounding  the  remaining  terms.  Examine 


<  C(At)4  11^11^ . 

(5.265) 

Similarly, 

y “At  Oa^+9*^  a(n+l) 
^  2A 

n— 0 

<  yC(At)4  [tn+1  \\att(;t)\\2  dt 

n= 0  Jtn+8 

(5.266) 

<  CiAt^Wauiy. 

(5.267) 
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Next  consider  the  convective  term,  and  use  lemma  C.2 


E  y  Ou (n)  •  Vff(n)  +  0«(n+^)  •  Vo-(n+0')  -  u(n+^)  •  Vcr(n+5) 


n=0 


n=0  Jtn 


l 


E 

n=0 


(At) 

48 


4  /•tn+ 1-9 


|«  •  V<rtt  +  2ut  •  Verj  +  utt  •  V<t||2  dt 


< 


< 


< 


ftn+i-e  /  \ 

y  /  (^4  ||w  •  V cr 1 1 2  +  4  ||2wt  •  Vcrt||2  +  4  ||wtt  •  V<r||2J  dt 

n= 0  ^tn 

E  y  JY"  (4  IHI~ d 11  v<r“112  + 16  l|u|l“  ** Ir Va,]l2 + 4  l|u|l“ rf  l|v<r| 


E 

n= 0 


(At) 

48 


4  /•tn+i_0 


4M2d  || Vo-ttf  +  16M2d  ||V<Ttf  +  4M2d  ||Vtrf  )  dt 


r‘2,1 


2J  I 


<  (At)  C  y|o-tt||0)1  +  ||crt||0  !  +  ||er||0 j y  . 
Again  using  lemma  (C.2)  the  ga(-,  •)  terms  are  handled  as: 


At 


YY-  0ga(<rW,VuW)  +  6ga(t 


~  9a{<T 


n= 0 


l  t  a  ,\4  /*tn+l 


< 


E 


(AO4 


n=0  48  y+0 


ll(5a(o-,  Vu))tt||2  dt 


l  ftn+l 

<  E(At)4Cr  /  ||(crVri)tt||2  dt 


n= 0 
Z 


'tn+e 

fN+1 


< 


< 


< 


/‘^n+i 

E(At)4Ci  /  ||cr(Vri)tt  +  2crt(\7u)t  +  <xtiVti||2  dt 

n=  0  -tin+e 

E(At)4C2  f  ^||cr(Vri)tt||2  +  ||crt(Vrt)t||2  +  ||<T«VM||2j  dt 

„ n  Jtn- \-G 


1  rtn+ 1 

E^*^  / 

n=0 

(ll(Vti), 

l 

r  (»■ 

E  ( At) 4  C2  M2  d2 

n=0 

An  +  fl  ^ 

(At)4C*  (||o-it||2j0  +  ||o-f||20 

tt 


IL  ^  ll^ll2  +  ll(V«)*||2o  rf2  Ikzf  +  II Vnll^  d2  \\att\ 


cr\\2  +  Ikzll2  +  ll^ttll2)  dt 


0,0  1  ii"  110,0  Z  • 


2)  dt 


(5.268) 


dt 


(5.269) 
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Lastly  for  the  Hg  terms  we  consider 


i 


E 

71= 0 


aAt 

~T~ 


d(Jn+5)) 


0d(u^)  -0d(w(n+0')) 


< 


< 


< 

< 


(A  t)*CCT, 


||d(w)tt||2  dt 
||(Vu)tt||2  dt 
||(Vw)it||^0  d2  ||1||2  dt 


where  Ct  is  a  constant  dependent  on  the  final  time  T. 


In  summary  we  have 


AtE 


n= 0 


a** 


<  C(At)4 


|  (Tttt\ 


0,0 


+  llcrii  llo,l  +  ll°"t 


10,1 


+  k 


1 2 

1 0,1 


+ 


(5.270) 


(5.271) 


5.3  Bounding  Hupa  and  H 


upB 


For  each  of  the  Hupa  terms  a  characteristic  term  is  treated  below.  All  parts  of  the  Hjnt  terms  will 
match  one  of  the  following  cases,  and  is  bounded  similarly.  Thus, 


S2M2d 


EAt 


n= 0 


>+§) 


<  C62 \\\(Tt\ 


1 0,0 


(5.272) 


Using  lemma  C.7  we  bound 


^  S2u9M2dAt 


71=0 


2A 


M+°)  -  o-k+l) 


<  ^2  PuOM2d{At)2 


ii=0 


2A 


f'<n+2  II  || 2  T, 

Ikdl  dt 


<  c(At)2k|k< 


1 0,0  • 


(5.273) 
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Using  lemma  C.8 


and 


E 

n= 0 


S29XM2dAt 


2A 


u 


(n+^)  •  Vcr(n+0')  -  u(n+0  •  Vcr(n+5) 


< 


E 

n=0 

l 

-  E 

n=0 

l 

£  E 

n=0 


S29XM2d(At)2 


2A 


(w  •  V<r)t||2  dt 


829XM2d{At) 

2A 


2  rtn+1_g 


|| Ut  ■  Ver  +  u  ■  VctjH2  dt 


n+  4 


S26XM2d(At)2  ftn+ i-e 


A 


r tn+1-e 

/  ||iti  •  Vcr |p  +  ||u  •  Vcrt |p  d 

Jtu+h 


<  C52(At)2  (||er||2 a  +  ||<xt 


io,i  /  > 


829M2d(At) 

2-J  9 


n=0 


5a(>+0,V>+*))-9a(>+0, 


< 


E 


S26M2d(At) 


2  /■tn  +  l-S 


ll(5a(o-,  Vw))t||2  dt 


n=0 

l 

<  ^Ci(52(At)2 

n=0 

l 

=  ^C1(52(At)2 

n=0 

l 

<  ^C282(Atf 

n=0 


It  ,  i 

n+2 


In+l-O 


||(<rVu)t||2  dt 


't  m 


zn+l-6 


\\crt\7u  +  a\/ut ||2  dt 


ln+l-6 


HctjVmII2  +  ||<rViif ||2  dt 


<  CS2(At)2  (  ||<Tt||o0  +  IMIon 


Lastly  using  lemma  C.14 


E 

71=0 


S2iaM2d(At ) 


S  E 

71=0 

l 


A 


d(w(n+0-iO)-d(«(n-0)) 

||dO)i||2  cit 


529aM2d(At)2  f^n+e-t 


itn-0 


^  829aM2d(At)2  fK+e-t  2  , 

<  E - E — -  /  iKv^di2  dt 


n=0 
s:2/  a  +\2/ 


<  82{AtyCCT. 


(5.274) 


(5.275) 


(5.276) 
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This  gives 


(HiipA  +  Hupb )  <  CS2(At)2  ^||ct||q j  +  ||<Tt ||o,i  +  II II 0,0  +  II  ^ 


n= 0 


1 2 

1 0,0 


+C52\\\at\\\l0  +  62(At)2CCT. 


(5.277) 


5.4  Apply  the  discrete  Gronwall’s  lemma 


Using  the  defined  bounds  on  the  H  terms,  the  lemma  A.l  and  lemma  A. 2  for  bounds  on 


p(9) 


and 


p(9) 


,  and  assuming  At  and  h  satisfy  the  relationship  given  in  (5.259)  the  discrete  Gronwall’s 


lemma  is  applied  to  (5.254)  giving 


pl~\~  1 


+ 


?i+ l-e 


+ 


?i+e 


+  AtY^  \  ^  -  26)  (SX  -  e\52  (2  +  3A  +  2a)) 


n= 0 


u(n+9)  .  v iKn+^) 

<  Go- {At,  h,  5),  (5.278) 


where 


Ga(At,h,5 )  C{At)4  ^  ||cr Ho,0  +  ll<Tttllo,1  +  Iloilo, i  +  Iloilo, i 


+  ll^ttllo.o  +  Iktllo.o  +  Iloilo, 0  +  Ct 


+  C{At,)~52  (llo-llo.i  +  II II o,i  + 
+  C  ( h2m+ 2  +  h2m  +  52h2m) 


,  1 2  i  ||  1 1 2  i  s~i 

<T||n  n  +  ||Crt||0,o  +  O T 


10,0 


l°"lllo,m+l 

+  Ch2m'+2\\at\\lm+1+  C52\\\*t\"2 


io,o  • 


Using  (5.278)  it  can  be  seen  that 


\F\\\2^0<Ga(At,h,d), 


1 0,0 


<TGa{At,h,S), 


and 

Then 


\u-VF\\\l<CGcr{At,h,6). 


\a  O'aIIIoo.O  —  Ills'll  loo, 0  +  lll^llloo.o 


<  Ga{At,h,5)  +  Ch2m+2\ 


loo,0  ’ 


(5.279) 

(5.280) 

(5.281) 


(5.282) 


\<r  ~  <Th\\\0,o  S 


< 


l^lllo,o  +  lllr"'2 


0,0 


<  TGa{At,h,S)  +  Ch2m+2 


l0,m+l  ’ 


(5.283) 
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and  (see  (4.11)) 


m-v(o--o-a)  mi  = 


< 

< 


|||n-V(F  +  r)|||| 

|||w  VF  +  u- VT)|||| 
|||wVF||||+|||wVr)|||| 
CGa(At,h,S)  +  Chm\\\a\\\ 


establishing  the  estimates  in  Theorem  4.1. 


(5.284) 


6  Proof  of  Theorem  4.2 


In  this  section  the  proof  of  Theorem  4.2  is  established.  The  proof  is  done  in  a  manner  similar  to  that 
used  for  Theorem  4.1.  For  notational  simplicity  the  ‘tilde’  is  dropped  from  Uh  and  p,  the  discrete 
approximations  of  velocity  and  pressure  using  the  true  stress.  We  start  by  examining  appropriate 
linear  combinations. 


Linear  combination  1:  Consider  the  linear  combination  of 


6  x  (4.3) 


+  (1-2 6)  x  (4.4) 


t - tv 7 


+  e  x  (4.7) 


t - tv 7 


(6.1) 


t=t„ 


When  we  multiply  the  expression  resulting  from  (6.1)  be  At  (done  in  a  later  step)  we  have  a  unit 
stride  of  size  At  between  tn  and  tn+i  with 


Re 


u 


(n+l) 


—  U 


(n) 


At 


+  6  ^n+e\  d(v))  +  6  (*(n+1\  d(v))  +  2 a6  (d{u^+1)),  d(v)) 

V  v  <E  Zh. 


+2 a~e  (d(4n+e)),d(v))  =  e  (f(n+0),  v)  +  e  (Vn+1\v)  , 

=►  B2  {uh ,  v)  =  6  (f (n+0) ,  v)  +  e  (f  (n+1) ,  v) 


(6.2) 

(6.3) 


Terms  in  (6.3)  are  defined  by  those  in  (6.2).  The  true  solution  of  the  momentum  equation  evaluated 
at  the  midpoint  between  tn  and  tn+ \  satishes: 


0  =  ^f(n+2);  +  (p(n+2\  V  •  v j  —  Re  v^j 


—  2d  (u(n+^  ,  d(v)^  —  ^cr(n+2),  d(v)^  ,  Vv  6  Zh.  (6.4) 
Adding  B2{u,v)  to  both  sides  (6.4)  yields 

B2(u,v)  =  (f(n+iO,v)  +  (p(n+^\  V  •  v)  +2^(v),  Vve4  (6.5) 


where  the  term 
2j(v)  :=  Re 


(dtu(n+1)  -  4n+^,v)  +  (tf<x(n+0)  +  0<r(n+1)  -  o-(n+§) ,  d(v)) 

+  2d^d(«(n+1))  +  0d(w(n+0))-d(w(n+5)),d(v)),  Vv£  Zh.  (6.6) 
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Subtraction  of  (6.2)  from  (6.5)  gives  the  error  equation: 

B2{eu,v)  =  (f(n+5)  -0f(n+0)  -0f(n+1\v)  +  (p(n+i),V-v)  +  X4  (v),  Vve4  (6.7) 
This  procedure  is  now  repeated,  setting  up  two  more  error  expressions. 


Linear  combination  2:  Using  the  linear  combination  of 


6  x  (4.3) 


+  (1  -  26)  x  (4.4) 


t=t„ 


+  6x  (4.7) 


t=tn 


t — tn—1 


(6.8) 


will  provide  a  unit  stride  between  tn-g  and  fn+i-0  when  we  multiply  the  resulting  expression  by 
At  in  a  later  step. 


+  02d(dC4n)),d(v))  =^(f(n+0),v)  +0(f(?l\v)  ,  Vve4  (6.9) 

=►  m K,v)  =  0  (f(n+V)  +e  (f("\ v) .  (6.10) 

Again  the  terms  in  (6.10)  are  defined  by  those  in  (6.9).  The  true  solution  of  the  momentum  equation 
evaluated  at  the  midpoint  between  tn_g  and  tn+i_g  satisfies: 


0  =  (f  ,  v)  +  (p(n+12-<>) ,  V  •  v)  -  Re  (4n+l  d) ,  v) 

—  2d  ,d(v))  —  d(v))  ,  Vv£  Z^.  (6.11) 

B2(u,v)  is  now  added  to  both  sides  (6.11)  giving 

Bg(u,v)  =  (f(n+M,v)  +  [p(n+^~e),V-  v)  +1 £(v).  (6.12) 

Here  we  define 


2£(v)  :=  (4u(n+0~)  -  4"+*  4  +  (0<r(n+(?)  +  -  cr(n+|-e) ,  d(v)) 

+  2d4d(M(n+e))  +  0d(-u(n))-d(w(n+i-e)4d(v)),  Vve4  (6.13) 
Lastly  the  error  expression  is  obtained  by  subtracting  (6.9)  from  (6.12)  giving 

H£(ew,v)  =  (f(n+12~e)  -0f(n+e)  -0fW,v)  +  (p(n+i-e),V-v)  +  2£(v),  Vve4  (6.14) 
Linear  combination  3:  Multiplying  the  linear  combination  of 


0  x  (4.3) 


+  (1-20)  x  (4.4) 


t — 


+  6x  (4.7) 


t — tn  —  1 


t — tn  —  1 


(6.15) 
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in  a  later  step  by  At  provides  a  unit  stride  between  tn+g_i  and  tn+g. 


Re 


+  6 


/  (■ n+0 )  _  (n+0- 1)  \ 

[- - A T - ’Vj  +2“0  (d(4n+0))>d(v))  +  2d0  (d(4n)),  d(v)) 

+  2d  (1-20)  (d(u(£+0~1)),d(v)')  +  0(<x(r!+e\d(v))  +(1-20)  ^o-(n+0_1),  d(v)^ 
(Vn),d(v))  =  0  (f  ("+*),  v)  +  (1  -  20)  +0  (f(n),v)  ,  Vv£  Zh. 

=>  ^K,v)  =  0(f(n+0),v)  +(1-20)  (f('!+M))v)+0(fW)v)  . 


(6.16) 

(6.17) 


The  terms  in  (6.17)  are  defined  by  those  in  (6.16).  The  true  solution  of  the  momentum  equation 
evaluated  at  the  midpoint  between  tn+e- 1  and  tn+g  satisfies: 


0=  (f(n+0~i),v)  +  (p(n+0-5),V-v)  -Re(u\n+(>  *\v) 

—  2d  ^d  (u(n+e~2^  ,  d(v)j  —  ^(j(n+0_2))  d(v)^  ,  Vv£  Z^.  (6.18) 
Adding  £>g(w,v)  to  both  sides  of  (6.18)  we  obtain: 

06V,v)  =  (f(n+0“^,v)  +  (p(n+0-i),V-  v)  +  2£(v),  Vv  G  Zh,  (6.19) 


with 


2£(v)  :=  Re  (dtu{n+9)  -  u[U+ °  ^ +  ( 6ain+ +  (1  -  20)  ^n+e~^  +  0o»  -  tr(n+°- s) ,  d(v)) 
+  2d  (0d(w(n+0))  +  0d(M(n))  +  (1  -  20)  d(w(n+0"1))  -  d(-u(n+f^)),d(v))  ,  V  v  g  Zh.  (6.20) 
Subtracting  (6.16)  from  (6.19)  gives  the  final  error  expression  for  the  conservation  equations. 

B$(eu,  v)  =  (f(n+e-^)  _  0f(«+0)  -  (1  -  20)  fb^-1)  -  0fW,  v) 

+  (pKI-^V-v)  +  2£(v),  Vv6  4  (6.21) 

In  order  to  complete  the  analysis  bounds  are  now  found  on  (6. 7), (6. 14),  and  (6.21),  with  the  value 
of  the  test  function  v  set  as  E^n+l\E^n+9\  and  E^n+e^  respectively. 


6.1  Bounding  B2{eu,  E{n+1)),  B%(eu,  E^6)),  and  B%{eu,  E(n+e)) 


In  order  to  complete  the  analysis  for  the  conservation  equation  we  need  to  consider  bounds 
on  (6. 7), (6. 14),  and  (6.21).  In  order  to  consider  the  bounds  on  the  terms  in  B\ (eu,  E^n+1^), 
BZ(eu,E(n+°)),  and  B^(eu,E^  n+^)  we  use  the  following  estimates: 

Re  (dte{u+1\  £(n+1))  =  Re  (dtE{n+1\  £(n+1))  +  Re  (dtK{n+l\  E[n+l)^j  ,  (6.22) 


with  bounds 


Re 


(dtE{n+1\  £(n+1))  >^f  (y  E(n+l) 


(6.23) 
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and 


Similarly, 


Re 


Re  (dtA{n+1\  E^n+1^  <  ^  ||(itA(n+1)||Z  +  i?e||^(n+1) 


and  is  bounded  as 


Re  (dtE(n+e\  E(n+e)^)  > 


and 


We  bound 


as 


Re 


(dtA(n+0),  E(n+§^ 


< 


Re 
2A t 

Re 


E{n+e) 

dtA(n+°) 


E(n-e) 


+  Re 


E 


(n+e) 


Re  (dtE^n+e\  £(n+0))  > 


and 


Re 
2A t 

Re 


j^in+O)  _  f](n-\-Q—  1) 


+  Re 


E  (n+0) 


and 


2a6  (d{\(n+V),d{E(n+V)^  VA(,,+1)  ' +  2 dde2  d(£(n+1)) 


4e2 


2d6  (d(E^n+e'>),d{E^n+1'>)^  <  ^  V£(n+0) 


+  2d0e2 


d  (£;(n+1)) 


and 


2d0(d(A(n+0)),d(£;(n+1)))  <  ^||vA(n+0)|p  +  2d0e2||d (E^n+1^) 


4e2 


(6.24) 


dtJ£+(>\  E(n+d)'\  =  Re  (dtE(n+e) ,  E(n+e)^j  +  Re  (dt A(n+f?) ,  #(n+(?))  ,  (6.25) 


(6.26) 

(6.27) 


Re  (dte{u+<>\  £(n+0))  =  Re  (dtE^n+0\  E^n+0^  +  Re  (dtA^n+0\  E(n+0)^j  ,  (6.28) 


(6.29) 

(6.30) 


Re  (dtA(n+0\  E(n+0^  <  ^  d*A(n+0) 

For  the  deformation  tensor  terms  we  have, 

2d0  (d(e^+1)),d(F;(n+1)))  =  2 dd  (d[E{n+V),d{E{n+1^  +  2d6  (d(A(n+1)), d(^n+1)))  ,  (6.31) 
with  bounds 

2a0  (d(£(n+1)),d(F;(n+1)))  =2dd  d (_E(n+1)) 


(6.32) 

(6.33) 


The  term 

2  dd  (d(e&+<?)),  d(£(n+1)))  =  2d6  (d^™^),  d(£(n+1)))  +  2  dO  (d(A^+0^) ,  d(E^n+1^  ,  (6.34) 
and  is  bounded  such  that 


(6.35) 

(6.36) 


We  have 

(1-0)25  (d(e^+0)),d(£;("+e'))) 

=  (1  -  9)2a  (d{E(n+V),d{E(n+0))^  +  (1  -  6)2d  (d(A(n+0>),d (tf(n+*)))  ,  (6.37) 
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with  bounds 

(1  -  0)2 a  (d(E(n+^),d(E(n+^)^J  <  ^~^2a 

1 

(1  -  0) 2d  (d(A(n+e)),d (EM))  <  (1  ~^2' 

.,-4- 


Next 


VE(n+°)  2  +  (1  —  #)2de3  d(£(n+0))  \  (6. 


9)2d(d(A(n+^),d(£;(Tl+0')))  <  (1  2^  ||vA(w+e)||2+  (1  ^2a  ||v£;(n+e')||2.  (6.39) 

(d(e^}),d(£(n+e~)))  =  02d  (d(E^),d(E(n+d))^  +  02d  (d(A^), d(£;(n+0')))  ,  (6.40) 

92d(d{E^n)),d{E(n+§))Sj  <0d  V.E(n)  2  +  0d  VE^n+'d^  2,  (6.41) 


Splitting 
2 a.0  [d{e^+e^ 
with  bounds 


02d  (d(A(n)),d(£;(n+^))  <  #d  VA(n) 
),d(£(n+0)))  =  2d0  (d(E{n+d)),d{Ein+e)) 


V£(n)  2  +  0d  V^(n+e)  2, 

2  ~  2 
VA(n)  +0d  VE(n+e)  . 


Then 


±  iien 

2d0  (d{eiu)),d{E(n+e^  =  2d0  (d{E^),d(E^n+e^  +  2d6  (d(AW),d(£(n+s)))  ,  (6.46) 

with  each  terms  handled  as: 

2d0  (d{E^),d{E^n+6^  <  ^-  V£(n)  2  +  2d#e3  d{E(n+d))  2,  (6.47) 

2d6  (d(A(n)),d(£(n+0)))  <  ^  VAW  2  +  2d0e3  d(£(n+f?))  2.  (6.48) 


(6.38) 


(6.40) 


=  2d0  (d{E{n+d)),d{Ein+9))^j  +  2d6  (d(A(n+0)),d(.E(n+0)))  , 
(d{E^n+^),d{Ein+^)^j  =2d0  d{E^n+e))  2, 


(6.41) 

(6.42) 

(6.43) 

(6.44) 


>  (d(A(n+0)),d(.E(n+e)))  <  VA(n+f?)  2  +  2d0e3  d(£(n+0))  2.  (6.45) 


Lastly 


2d  (1  -  20)  (d(e£+0_1)),  d(£(n+0))) 

=  2d  (1  -  20)  d(£(r!+e)))  +  2d  (1  -  20)  (d(A(n+0_1>),  d(£(n+0)))  ,  (6.49) 

where 

2d  (1  -  20)  (d(.E(n+*_1)),d(.E(n+fl))) 

<  a(yl  ~  2 ^  V£l(n+e_1)  2  +  2d  (1  -  20)  e3  d (E{n+e))  2,  (6.50) 
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and 

2d  (1  -  20)  (d(A(n+fl-1)),  d(£;(n+0))) 

<  a{l~2^  VA (n+e~V  +  2d  (1  —  20)  e3  d (E{n+e))  2.  (6.51) 
2e3 
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6.2  Bounding  T'l (_E(n+1')) 

Bounds  on  the  interpolation  error  are  now  considered  for  I2(E^n+1'>)  (see  (6.6)).  Thus 


I£(£(n+1))  :=  Re  ^dtn(n+1)  -  u^\ E^n+1^  +  (Wn+^  +  6^n+^  -  cr(n+5) ,  d(^(n+1^)) 

+  2d  ^d(w(n+1))  +  6d {u^n+e^)  -  d  ( u(n+ s))  ,  d(£(n+1)))  .  (6.52) 


We  bound  each  of  these  terms  as 


^cWn+1)  -  u{n+^\  E(n+1A  <  ^  dt«(n+1)-4n+")  +  -Re  ||£i^n+1^  || 2  ,  (6.53) 


i+0)  +  6a^n+1)  -  o-(n+i),d(^(n+1))) 


<  —  6c An+d^  +  0cr(n+1)  -  cr(n+§)  +  e2  d (E^ 

4e2 


,  (6.54) 


2d  (6»d(u(n+1))  +0d(u(n+e))  -  d  (u(n+^))  ,d(£(n+1>)) 

<—  6d{u^n+1^  +  6d{u^n+^)  -d  fu(n+2))  2  +  2de2  d{E^n+1^)  \  (6.55) 

6.3  Bounding  l£(E(n+°)) 

Similar  to  the  previous  section  we  now  bound  the  term  (see  (6.13)).  Here 

1 2{E^n+~e))  :=  Re  (dtu(n+°)  -  u^~d\  E(n+^  +  (Wn+(?)  +  6 -  a^~e\  d(^(n+0'))) 

+  2d  4d(n(n+0))  +  6d(u^)  -  d  (u(n+^-0))  ,d(^(n+0')))  .  (6.56) 


Then 


Re  ( d, u(;+i)  -  u^'" 


4”+’  +  ife  £(’*+")  2,  (6.57) 


>+0)  +  gcrin)  _  a{n+\-e)  )d(_E(n+e)^ 


<-  6a{n+^  +  6a(n^  -  a(n+12~e)  2+  VE^ 


,  (6.58) 


2d  (dd{u(n+e^)  +  6d{u^)  -d  (w(n+^-0))  ,d(£l(n+0'))) 

(l-0)d(u(n+0))  +  #d(u(n))-d(u(n+i-0))  2  +  2d  V_e(’ 


(6.59) 
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6.4  Bounding  Zg  (£Kn+0)) 

Here  we  exarnin  the  bounds  given  for  the  Zg  term  (see  (6.20))  with  the  choice  of  v  =  E^n+0\ 
jn^E(n+0))  ■-  Re  ^dtu(n+V  -  u[n+(>~^ ,  E(n+eA 

+  (Wn+0)  +  (1  -  20)  <T(n+0_1)  +  0a(n)  _  cr(n+0-|) ,  d(£(n+(?))) 

+  2d  (0d(u(n+e))  +  0d(i*(n>)  +  (1  -  20)  d(u(n+0-1))  -  d(u(n+e-^)),  d(Zj(n+0)))  .  (6.60) 

Each  piece  is  bounded  as 

Re  (dtu(n+6)  -u\n+6^\E^n+e^\  <  ^  (Wn+(?)  -4”+0~^  +-Re  £(n+0)  2  ,  (6.61) 

(0<r(n+e)  +  (1  -  20)  a(n+0-V  +  0O-W  _  cr^-l),  d(£tn+0))) 

<  0Cr(n+0)  +  (l-20)cr(n+0-1)  +  0Cr(n)-Cr(n+0-i)  2  +  e3  d(£(n+0))  2,  (6.62) 

4e3 

and 

2d  (0d(«(n+0))  +  0d(«W)  +  (1  -  20)  d(n(n+0-1))  -  d(n("+e-2)),  d(E^n+0^ 

<—  9d{u(n+V)  +  6d{u^)  +  (1  -  29)  d{u(n+0~V)  -d(u(n+e-^))  2  +  2de3  d(£;(n+*>)  2. 
2e3 

(6.63) 


6.5  Bounding  the  f  and  p  terms 

Bounds  for  the  terms  containing  the  forcing  function  f  are 

(f (n+§)  _  §f(n+d)  _  0{(n+l)  E{n+l)^  <  1  f(n+i)  _  £f( 
^f(n+i_0)  _  Qf(n+0)  _  Qf (n),E(n+0)^  <  I  f(n+l-0)  _  ^ 


<  I  f(n+i 
“  4 


i)  -  0f(n+(?)  -  0f(n+1)  2+  E(n+1)  2,  (6.64) 


<1  f("+5-0)  -  0f(n+O)  _  0f(")  2+  E^n+0)  2,  (6.65) 


f(n+0-§)  _  gf(n+9)  _  (1  _  20)  f(n+0~1)  _  0f(n);  E^n+0^ 

<  I  f(n+e-|)  _  0f(n+e)  -  (1  -  20)  f(n+6l_1)  -  0f(n)  2  +  _E(n+0) 

;  pressure  terms  are  bounded  using 

(p(n+^),  V  •  _E(n+1^  =  (p(n_4)  -  p(n+i);  V  •  £(n+1)^ 

<  p(n+l)  --p(n+ 1)  S  vs(n+1^  2 

^  ^  IL(«  +  4)  -7-)(™+l)!|2  I  Ilv7  B’(n+1)  l|2 


(6.66) 
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>  +  5-0) 


y -e 


<—  p(n+h~e)  -p(n+  h~9)  2  + 
~  4 


(6.68) 


(p(n+e~^\  V  •  E(n+p  <  p  p{n+B-\)  -'p(n+e~ I)  2+  v^(n+e>)  2. 


(6.69) 


6.6  Bring  all  the  parts  together  for  the  conservation  equations 


Using  the  estimates  given  in  (6.22)  -  (6.69)  a  single  error  expression  is  given  by 


f  E^n+1)  2  -  E 2  +  e( 

2A t  V 


n+e)  2  _  2  _J_  ^(n+0)  2  _ 


+  [2a0  -e2  (50  +  25  +  450  +  lJJ  d(^n+1J) 

+  (2d0  —  e3  (25(1  —  0)  +  6d0  +  4a  (1  —  29)  +  2a  +  1))  d (E^n+e^) 

<  (2Re  +  1)  ||£;(n+1)||2  +  (2Re  +  1)  ||.E(n+e)||2  +  (2Re+  1)  ll^7^  if 


+  (da  +  VE W  2  +  ( 1  +  ^  j  V£l(r!+e)  2  + 


2  (  5(1-20) 
1  A, 


Y^ln+e-i) 


^2  +  (1  +  9)2 a  +  (1  +  (1  || X7E(n+°)  ||2  +  ||v^(n+1) 


+  //A  +  H  T"  +  tiff  +  Hp 


where 


if  A  :=  ^  N|dtA(”+1)  ||2  +  ||dtA(n+0~)||2  +  |dtA(n+0)||2 


2d0 

+  —  VA 
4e2 


(n+i)  2  +  fe&  +  VA(n)  2  +  VA(n+e- 

V  4e3  y  2e3 


/  2a6  (1  -  0) 2a  d9 

+  [  4e^+  2  +  2es 


+  —  VA(»+0)  , 

2e3  / 


(6.70) 


(6.71) 


dtu^  -  u\n+^  + 


in+^-e)  2+  dtuM  -  tS1**-® 


+  —  9a(n+V  +9a(n+V  -  cr(n+l)  +  -  0o>+0)  +  9a^  -  a(n+^~0 

4e2  4 

+  —  0o-(n+0)  +  (l-20)cr^+e”h  +  0crW_(T(n+e-i)  2 

4e3 

+P  ||0d(w(n+1))  +  9d{u^n+0))  -  d  (Vn+^))  ||2 

+  “  || (1  -  0)d(n("+e))  +  0d(i*W)  -  d  (™(n+^))  ||2 

+  —  ||0d(w(n+^)  +  9d{u^)  +  (1  -  29)  d{u(n+e~V)  -  d(u(n+0-5))||2  . 


(6.72) 
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Hf  :  = 


1 

+  4 


f(n+i)  _  §{in+d)  _  0{(n+l)  2  _|_  1  f{n+h~e)  -  §{(n+s)  _  (9f(n) 

4 

f(n+6»-|)  _  Qf(n+e)  _  (1  _  29)  f(n+e,_1)  _  0f(n) 


(6.73) 


and 


H  -  — 
hp  —  4 


^  p("+5)  _£>(™+0  V  p(n+h~d)  -J>(n+k-e)  2+  p{n+9~h)  _p(n+0-§)  ^ 


Multiplying  expression  (6.70)  by  and  summing  from  time  n  =  0  to  l  gives 

_g;(z+i)||2  _  ||_e;(0)||2  -g  II  £0+0)  II 2  _  ll^Wlh  4_  |£0+0)||2  _  II £(e) 


(6.74) 


+  (4a9  -  e2  (2a0  +  4d  +  8a9  +  2))  ^  ^  ||d(£;(n+1)) 

77.— 0 

At  1 

+  (4d0  -  e3  (4d(l  -  0)  +  12d0  +  8d  (1  -  20)  +  4d  +  2))  —  V 

ire  ' 

n= 0 

<  AtJ]||.E(n+1)  2  —  ((2Re  +  1)  +  Ch~2)  +  AtJ^  ||^(n) 

n=0  n=0 

+A II  E(n+§) 


d(E^n+e)) 


n= 0 

2  2Ch~2 


9a 


a9 

2es 


77. =0 

l 


+A  t^  £1 


l(n+0) 


n= 0 
l 


2  2 

Re 

^{2Re  +  1)  +  Ch~2 

2  —  | 
Re 

( (2 Re  +  1)  +  Ch~2  | 

2d9 

4e2 


+AtVllR(n+0-1)  2  2(~^  ^O"2^ 

^11  2i?ee3 


+At  E  +  At  E  i~eHln  +  At  E  +  At  E  ^eHv  (6-75) 

77.— 0  72=0  72=0  72=0 

In  order  to  apply  the  discrete  Gronwall’s  lemma  we  need  the  inequalities  stated  in  (6.76),  (6.77), 
and  (6.78)  to  hold. 

((2 Re  +  1)  +  Ch~2)  +  Ch~2  ( ' 9d  +  Y, 

[{2 Re  +  1)  +  Ch~2  [2  +  —  + 


1  > 

1  > 


1  > 


2At 
Re 
2  At 
Re 

2  At 
Re 


4e2  J 

Summarizing  (6.76),  (6.77),  and  (6.78)  yields  that 

At  (C  +  Ch~2)  <  1. 

This  gives  the  stability  condition 

At  <  Ch 2. 


(6.76) 

(6.77) 

(6.78) 


(6.79) 
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6.7  The  H  terms  in  (6.75) 


Using  interpolation  properties,  and  elements  of  order  k,  and  q  for  velocity  and  pressure  respectively, 
bounds  are  found  for  the  and  Hp  expressions  in  (6.75).  Here 


<  Ch2t+2  KIIJU+J  +  Ch2k 


u 


lo.fe+l  • 


n= 0 


Bounding, 


(6.80) 


o  At 

-  C(Af)4  ||«(ft||o,o  +  C{At)1  ||(Trt||o,0  +  C(Ai)4  ||(Vm) 


2 

tt  ||  0,0 


n= 0 


<  C(At)4  \\uttt\\l0  +  C(At)4  ||(Vu)tt||gi0  +  C(At)4CT 

<  C(At)4  \\um\\lfi  +  C(At)4  WuuW2^  +  C(At)4CT. 

For  the  forcing  function  and  the  pressure  we  have: 

i 


T.  Ah<  S  C(At)4  ||f, 


n= 0 


2 

til  1 0,0  > 


and 


]T  <  Ch 2q+2 


n= 0 


Re 


1 2 

lo,qr+l  ' 


Applying  the  discrete  Gronwall’s  lemma  to  (6.75)  we  obtain 


E 


rZ  +  1 


E(0) 


+ 


£{i+e) 


E& 


+ 


£jO+0) 


E(d) 


+ 


At 


(4 ad  —  2e2  (6a  —  3 ad  —  1))  —  d(.E 

Re 


'(«+!)' 


At 


n= 0 
l 


+  (4a£-2e3(8d-4d0  +  l))  —  V  d (E^^) 

Re  z — J 


n= 0 


<  Gu(At,h,8), 


where 


Gu(At,h,8)  :=  Ch2k+2\\ut\\2+Ch 


2k 


0,fc+l 

+  C(At)4  ||wttt|lo,o  +  C(At)4  \\uuWl,  +  IMo, 

+  C{At)4CT. 


u 


l0,fc+l 


+  Ch2q+2 


0,g+l 
,0 


(6.81) 


(6.82) 


(6.83) 


(6.84) 


(6.85) 


Applying  lemmas  5  and  A. 3  for  bounds  on 


E (0) 


and 


E^  respectivly,  and  assuming  E^0) 


II  II  I  I  I  I  II  II 

0,  then  from  (6.84)  we  have  (6.86).  Lemma  5  and  lemma  A. 3  are  stated  and  proved  in  Appendix 
A. 

(6.86) 

(6.87) 


which  gives 


\E\\2  +  C\\\d(E)\\\l0<Gu(At,h,5), 
\\\E\\\2<Gu(At,h,5), 
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and 


|  -E  lo,0  <  TGu{At,h,5). 

(6.88) 

Using  Korn’s  lemma  with  (6.86)  gives 

\\\WE\\\l0<CGu(At,h,8). 

(6.89) 

Hence 

|||£||lo,i  <CGu(At,h,6). 

(6.90) 

To  establish  the  error  estimate  in  Theorem  4.2  note  that 


u  U^IHoo,0  —  I  Ill'll  loo, 0  ""t  I  ll^ll  loo, 0 

<  Gu(At,  h,  6)  +  Ch2k+ 2  |  ||u||  &ifc+1 ,  (6.91) 

ll«-  uh\\\l,i  <  |||-E||lo,i  +  !llAlllo,l 

<  CGu(At,h,6)  +  Ch2k\\\u\\\lM1.  (6.92) 


Justification  of  (IH1)  for  U)x. 

Assume  that  (IH1)  holds  for  n  =  1,  2, . . . ,  l  —  1.  Then  using  inverse  estimates,  interpolation  prop¬ 
erties,  and  (4.15) 


uh 


< 


< 


ulh  -  ul 


+ 


u 


El 


Ch 


+ 


E1 


A1 


+  M 


+  Ch~£ 


A1 


+  M 


<  C  (hk~i  +  hq~ 5+1  +  (A t)2h^  +  hk~i+l^j  + 


M. 


(6.93) 


Setting  k  >  |,  and  q  >  i  —  1  and  choosing  h,  and  At  such  that 


hk~i,hq~i+1  <  and  At2  < 

(_y  O 


we  have 


Similarly  it  follows  that 


uh 


<  M  +  4. 


u 


(n+0) 


<M  +  4. 


(6.94) 


7  Numerical  Results 

In  this  section  we  present  numerical  results  for  the  0-nrethod  on  two  test  problems.  The  first 
example  uses  a  known  analytical  solution  to  verify  numerical  convergence  rates  (Cvge.  Rate)  for 
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the  0-method.  In  the  second  example,  we  considered  a  prototypal  problem  of  viscoelastic  flow,  flow 
through  a  4:1  planar  contraction.  In  both  examples  finite  element  computations  were  carried  out 
using  the  FreeFem++  integrated  development  environment  [13].  Continuous  piecewise  quadratic 
elements  were  used  for  modeling  the  velocity,  and  continuous  piecewise  linear  elements  were  used 
for  the  pressure  and  stress.  The  constitutive  equation  was  stabilized  using  a  SUPG  discretization 
with  a  parameter  <5. 

For  the  (optimal)  value  of  0  =  1  —  V2/2  ~  0.29289  the  local  temporal  discretization  errors  are 
0((At)2).  The  influence  of  the  value  for  0  on  the  numerical  approximations  is  illustrated  in  Figure 
7.1.  for  computations  performed  on  Example  1  (described  below).  All  the  other  computations 
reported  in  this  paper  were  obtained  using  0  =  1  —  \/2/2.  The  constitutive  equation  splitting 
parameter  oj  was  set  to  1/2  in  all  computations. 


(a)  Error  ||it  —  uLHo.o  as  a  function  of  6.  (b)  Experimental  convergence  rates. 

Figure  7.1:  Optimal  0  value 


7.1  Example  1 

In  order  to  investigate  the  predicted  convergence  rates  we  consider  fluid  flow  across  a  unit  square 
with  a  known  solution. 

Let  12  =  (0, 1)  x  (0, 1),  Re  =  1,  a  =  1/2,  A  =  2,  and  a  =  1.  For  the  true  solution  we  use 

=  (  e^-i‘)(x^-xW-y)  \  ,  . 

y  —  e(x+y  x 2  —  x)(y2  —  y)  J  ’ 

p  =  cos(27 xx){y2-y),  (7.2) 

a  =  2ad('u).  (7.3) 

Remark:  A  right-hand-side  function  is  added  to  (2.1)  and  f  in  (2.2)  is  calculated  using  (7.1)-(7.3). 
For  this  example  three  sequences  of  computations  were  performed: 

(i)  approximation  of  ah,  assuming  u  and  p  (Theorem  4.1), 

(ii)  approximation  of  Uh  and  ph,  assuming  cr  (Theorem  4.2), 

(iii)  approximation  of  Uh,Ph  and  ah- 
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7.1.1  Approximating  a ^  with  u  and  p  known 

We  first  consider  the  approximation  of  the  stress  &h  assuming  the  velocity  and  pressure  functions 
are  known.  This  is  analogous  to  implementing  Step  la,  Step  2b,  and  Step  3a  of  the  0-nrethod. 
From  Corollary  1  we  have  predicted  asymptotic  convergence  rate 

1 1| tx  •  V  (<r  —  &h)  || \g  <  C  ((At)2  +  A t8  +  hS  +  h  +  8) , 

which  is  consistent  with  the  numerical  convergence  rates  presented  in  Table  7.1.  Note  Table  7.1 
shows  the  effect  of  the  upwinding  parameter  5  on  |||<r  —  &h |||0q- 

For  the  proof  of  Theorem  4.1  we  required  the  restriction  At  <  Ch2.  The  numerical  computations 
were  performed  with  At  ~  h.  It  is  an  open  question  if  the  restriction  At  <  Ch 2  is  a  necessary 
condition  for  (4.12). 


Table  7.1:  Approximation  errors  and  experimental  convergence  rates  at  T  =  2 


Si 


(At,  h) 


1  \/2 
2 


1  \/2 
2’  4 


1  j/2 

4  *  -8 


i  ■jM 
8’  16 


J_  \/2 
16 !  32 


~JlT 

C2 


ll0-  O'hlllo.O 
Cvge.  Rate 


2.1235e-l 


6.6773e-2 

1.7 


1.9191e-2 

1.8 


5.0437e-3 

1.9 


1.2830e-3 

2.0 


\u-V{cr  -  ah)\\\§ 
Cvge.  Rate 


1.1544e-l 


6.1759e-2 

0.9 


2.9979e-2 

1.0 


1.5185e-2 

1.0 


7.6944e-3 

1.0 


Im  tThlWofi 
Cvge.  Rate 


2.0070e-l  8.4563e-2  3.7449e-2  1.6980e-2  8.0428e-3 

1.2  1.2  1.1  1.1 


\u-V{a  -  ah)\\\g 
Cvge.  Rate 


1.0263e-l 


6.4336e-2 

0.7 


3.2570e-2 

1.0 


1.6325e-2 

1.0 


8.2071e-3 

1.0 


_h_ \  2 
C2 


A.)2 

C2 


ll0-  o-hlllo,o 
Cvge.  Rate 


2.0174e-l 


7.3678e-2 

1.5 


2.3575e-2 

1.6 


6.9629e-3 

1.8 


2.0645e-3 

1.8 


\u-  V(er  —  &h)\\\g 
Cvge.  Rate 


1.0447e-l 


6.2239e-2 

0.7 


3.0497e-2 

1.0 


1.5304e-2 

1.0 


7.7199e-3 

1.0 


ll0-  o’hlllo.o 
Cvge.  Rate 


2.0346e-l 


6.9501e-2 

1.5 


2.0245e-2 

1.8 


5.3281e-3 

1.9 


1.3546e-3 

2.0 


\u-V(cr  -  arh)\\\g 
Cvge.  Rate 


1.0664e-l 


6.1737e-2 

0.8 


3.0104e-2 

1.0 


1.5203e-2 

1.0 


7.6968e-3 

1.0 


7.1.2  Approximating  Uh  and  ph  with  known  a 

Numerical  results  for  the  approximation  of  velocity,  Uh ,  and  pressure,  ph  for  a  known  stress,  er,  are 
presented  in  Table  7.2.  These  results  correspond  to  the  analysis  of  step  lb,  step  2a,  and  step  3b  as 
stated  in  Theorem  4.2. 

The  numerical  convergence  rates  observed  are  consistent  with  those  predicted  in  Corollary  2  where 

|||«  -  ithlllo,!  <  o  ((At)2  +  /i2)  . 
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Similar  to  the  case  in  the  previous  subsection  we  make  note  that  the  proof  of  Theorem  4.2  required 
the  restriction  At  <  Ch2 .  Here  the  numerical  computations  were  performed  with  At  ~  h. 

Table  7.2:  Approximation  errors  and  experimental  convergence  rates  at  T  =  2 


(At,  h)  -> 

Ilk  —  ^h\\ 1 0,1 

Cvge.  Rate 


7 1  Vl\  (\  (l  n~V2\  /  J  y/2\  7  J  V2\ 

\  2  ’  4  )  (4’  8  j  ^8’  16  )  /  16  ’  32  J  ^32’  64  J  64  ’  128  y 


4.4196e-2 


1.1485e-2 

1.9 


2.9707e-3 

2.0 


7.5759e-4 

2.0 


1.9142e-4 

2.0 


4.8129e-5 

2.0 


llU-^llloo,0 
Cvge.  Rate 


1.4734e-3 


1.7996e-4 

3.0 


2.5014e-5 

2.8 


4.1759e-6 

2.6 


8.5692e-7 

2.3 


2.1636e-7 

2.0 


\\\P  ~  Ph\\\o,o 
Cvge.  Rate 


1.0859e-l 


6.6842e-3 

4.0 


1.5033e-3 

2.2 


3.9097e-4 

1.9 


1.2086e-4 

1.7 


4.7884e-5 

1.3 


lb-Phllloo.0 
Cvge.  Rate 


8.4003e-2 


4.9703e-3 

4.1 


1.1343e-3 

2.1 


3.2878e-4 

1.8 


1.2797e-4 

1.4 


6.0659e-5 

1.1 


7.1.3  Full  0-method  approximation  for  viscoelasticity 

Table  7.3  contains  the  results  for  the  approximation  of  u  and  a  using  the  0-method  described  in 
Step  la  -  Step  3b.  The  numerical  convergence  rates  are  consistent  with  our  expectations  based  on 
Theorems  4.1  and  4.2,  i.e. 

Ilk  -  Mftlllo,!  +  Ilk  -  O'hlllo.o  +  Ilk  •  V(CT  -  ah) Hie  ~  o  ((At)2  +  (At)  8  +  5  +  hJ  . 


7.2  Example  2 

For  a  second  example  the  numerical  approximation  of  viscoelastic  flow  through  a  planar  4:1  con¬ 
traction  channel  is  presented.  This  has  been  a  long  standing  benchmark  problem  for  viscoelastic 
flow  [16,  19,  20].  A  diagram  of  the  flow  geometry  is  given  in  Figure  7.2.  It  is  assumed  that  the  chan¬ 
nel  lengths  are  sufficiently  long  for  fully  developed  Poiseuille  flow  at  both  the  inflow  and  outflow 
boundaries.  In  the  computations  the  value  of  L  in  Figure  7.2  is  set  at  1/4. 


16L 


Figure  7.2:  Plot  of  4:1  contraction  domain  geometry. 


The  flow  at  t  =  0  is  assumed  to  be  stationary  and  then  slowly  increased  for  t  >  0  using  A(t)  =  1— e  t. 
The  boundary  conditions  at  the  inflow  of  the  channel  are  defined  by 

u  =  .4(0  (  55  k  )  ,  (7.4) 
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Table  7.3:  Approximation  errors  and  experimental  convergence  rates  at  T  =  2 


Si 

nr 


“7T“ 

s/2 


(At,  h)  — ■> 

i  m  _  uk\\  1 0,1 

Cvge.  Rate 


1  s/2 

2’  4 

4.7608e-2 


1  s/2 
4’  8 

1.2323e-2 

1.9 


1  s/2 

8’  16 

3.2034e-3 

1.9 


_L  vd 

16’  32 

8.3187e-4 

1.9 


_L  vd 

32’  64 

2.1793e-4 

1.9 


Cvge.  Rate 


6.5569e-2 


1.9248e-2 

1.8 


5.0845e-3 

1.9 


1.3000e-3 

2.0 


3.2981e-4 

2.0 


\u-V(a-ah)\\\d 
Cvge.  Rate 


6.1217e-2 


2.9982e-2 

1.0 


1.5191e-2 

1.0 


7.6984e-3 

1.0 


3.8800e-3 

1.0 


|||W  Wft||lo,l 
Cvge.  Rate 


5.0150e-2 


1.6193e-2 

1.6 


6.4636e-3 

1.3 


2.9147e-3 

1.1 


1.3989e-3 

1.1 


m  °7iiiio,o 

Cvge.  Rate 


1.1922e-2 


3.5905e-2 

1.2 


1.5791e-2 

1.2 


7.3743e-3 

1.1 


3.5789e-3 

1.0 


\u-V(a  -ah)\\\§ 
Cvge.  Rate 


6.4087e-2 


3.2589e-2 

1.0 


1.6309e-2 

1.0 


8.2007e-3 

1.0 


4.1180e-4 

1.0 


_h_ \  2 

V2 


If" 

s/2j 


|||«  uh\\ lo,l 
Cvge.  Rate 


4.8620e-2 


1.3135e-2 

1.9 


3.6334e-3 

1.9 


1.0192e-3 

1.8 


2.9513e-4 

1.8 


m  o-ftiiio.0 

Cvge.  Rate 


7.1941e-2 


2.3392e-2 

1.6 


6.8309e-3 

1.8 


1.9976e-3 

1.8 


6.0598e-4 

1.7 


\u-  V(cr  -  crh)|||0- 
Cvge.  Rate 


6.1860e-2 


3.0510e-2 

1.0 


1.5305e-2 

1.0 


7.7221e-3 

1.0 


3.8851e-3 

1.0 


|||«  uh\\ lo,l 
Cvge.  Rate 


4.8022e-2 


1.2507e-2 

1.9 


3.2644e-3 

1.9 


8.4818e-4 

1.9 


2.2193e-4 

1.9 


m  o'hiiio.o 

Cvge.  Rate 


6.8104e-2 


2.0308e-2 

1.7 


5.3614e-3 

1.9 


1.3680e-3 

2.0 


3.4645e-4 

2.0 


\u-  V(o’-cr/l)|||0- 
Cvge.  Rate 


6.1279e-2 


3.0111e-2 

1.0 


1.5208e-2 

1.0 


7.7004e-3 

1.0 


3.8803e-3 

1.0 
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and 

XA(t)2y2a(l  +  a)  —16  aA(t)y  XA(t)2y2a(a  —  1) 

1711  = - W) - ’  =  W)  ’  a  1722  = - W) - ’  (  1 

where 

D(t)  =  256  +  A(t)2y2 A2(l  +  o)(l  -  a). 

The  outflow  boundary  condition  is 

=  2(“0_*2)  )•  (7-6) 

No  slip  boundary  conditions  were  imposed  for  the  velocity  on  the  solid  walls  of  the  contraction, 
and  a  symmetry  condition  was  imposed  along  the  bottom  of  the  computational  domain.  The 
computations  were  performed  on  a  uniformly  refined  version  of  the  mesh  shown  in  Figure  7.3  with 
Ax,nin  =  0.0625  and  A ym;n  =  0.015625. 


Figure  7.3:  Sample  Contraction  Mesh 


The  computations  were  done  using  the  full  0-method  approximation  given  by  (4.2)  -  (4.7)  for  an 
Oldroyd  B  fluid  (a  =  1),  with  A  =  2  and  Re  =  1.  The  value  of  a  was  set  to  8/9,  which  is  commonly 
used  in  the  literature  [19].  The  discrete  time  step  and  upwinding  parameter  were  set  to  At  =  1/32 
and  5  =  (2/A ymin)2.  Figures  7.4  and  7.5  show  streamlines  for  the  fluid  at  times  t  =  1/8, 1/2, 1,  and 
4  superimposed  on  a  contour  plot  showing  the  magnitude  of  velocity.  Note  that,  consistent  with 
expectations,  as  the  velocity  is  increased,  a  vortex  appears  in  the  upper  corner  of  the  domain  and 
grows  with  the  magnitude  of  the  velocity. 


0  0.002  0.004  0.006  0.008  0.01  0.012  0.014  0  0.005  0.01  0.015  0.02  0.025  0.03  0.035  0.04  0.045 


(a)  t  =  0.125  (b)  t  =  0.5 


Figure  7.4:  Streamlines  and  magnitude  contours  for  u 
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(a)  t  =  1 


(b)  t  m  4 


Figure  7.5:  Streamline  and  magnitude  contours  for  u 
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A  Initial  Bounds 


In  this  section  we  establish  error  estimates  for 


F (0)  ,  F 0)  , 


and 


used  in  the  analysis  of  the  0-method. 


A.l 


p(9) 


JT0) 


and 


u 


0)  .  yjrW 


Lemma  A.l 


_p(0) 


<  Ch2m+2  +  CAt2h2m+2  +  C52  At2  h2m 

+  CAt2h2m  +  CS2At2h2m~2  +  CAt4  +  CS2At2h~ 2  +  CAt4S2h~2. 


(A.l) 


Proof:  Subtracting  ((4.2)  with  known  u  and  n  =  0)  from  ((3.2)  at  t  =  0): 


A 

JA t 


(0) 


e(j  ,T)  +w  ey.r^o)  = 


A 


0At 
-  A 


d°) 

(T  ) 


,(°) 

(T  ) 


r)  -  (1  -  w)  (e§},  r5(o))  -  A  (u(n)  •  Ve^,rj(o)) 
V«(°)),t5(o))+A 


rW  _  n-(O) 


(0)  . 
t  J 


+  cj  ^er(0)  -  cr(0),  r^o)^  +  A  •  Vrj 

Recall  that  +  F^\  Let  r  =  F^  and  (A. 2)  yields 


F0)  Z  +  |  F(e)  *  +  ^ 


u>0At 


'pW 


+  - 


uOAt 


(f0\8uW  •  VF(0))  +  (r(0\ <5n(0)  •  VF(0)) 


A 


uOAt,  /  uOAt 


A 

uOAt 

A^ 


p(°)^  _ 


A 

uOAt 

A 


F(0),<5m(0)  •  VF(0) 
r(0),5u(0)  •  VF(0) 


-0Ai  (u(0)  •  VF(O),F(0))  -  0At  (V0)  •  VF(0),<5w(0)  •  VF(0) 
-0At  (u(0)  •  vr(0), -  6 At  (u{0)  •  Vr(0), du(0)  •  VF(0) 


-0Ai 

-0Ai 


^a(F(0),Vw(0)),F(e))  -0At  ^a(F(0\V«(0)),  W0)  •  VP(0)) 
r(°),  Vm(0)),  F(0))  -  0At  f5a(r(°),  Vm(0)),  <5u(0)  •  VF(0)) 


/ 


-0Atw  (V0)  -  ffC),#)  -  0Ahj  (V0)  -  <r^°\  5u^\/ F^^j 


+0A  t. 


-  ^°) 
6At 


(A.2) 


-crj0),F(eA  +  0Af  (-o-J0),W0)  ■  VF((?))  .  (A.3) 
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Appropriate  bounds  can  now  be  placed  on  each  of  the  terms  in  (A. 3).  Thus, 


+  ei 


F(e) 


co8At 

A 


< _  r(0) 

rw  F(0))  <  uj202(At)2  T(0) 

’  )-  4e2A2 

y),k(°).v^)=o, 


+  e2 


pW 


A 


^  (r^,  fe(°)  •  vfW)  <  f)'2  y2(/(  7'  yg2(At)2  ||rw 


4e3A2 

y(o),F(0))  =  0 

<  — —  r(0) 

V  ’  )  ~  4e3 

u8At 


+  C3 


p(9) 


A 


+  e4 

F(°),F(0))  =  0, 


p(d) 


CjQAt 

A 


F(0),<ki(0)  •  =  0, 


fr(0)  F(0)\  ^  u282{At)2 

p(0) 

2 

-\-  €5 

F(d) 

V  ’  /  4A265 

-2"r(°)|12 


^  (r(°),W°)  •  vf^)  <  PM2fo2P{M)2Ch 

6 A t(V0)  •  VF(O),F(0))  =  0, 

0At  (w(0)  •  VF(0),  Ju(0)  •  VF(0))  =  0, 

0A t  (u(°)  •  VT(°),fW)  <  M2d02^'2  |  |vr(0) 


eg 


_p(0) 


0At  (V0)  •  vr(0),  <5-u(0)  •  VF(0)) 


< 


4e7 

52d2M482{At)2Ch~2 

4e8 


+  €7 


vr(0) 


p(d) 


es 


jtA) 


0At  (ffa(F(°),Vu(°y^)  =0, 

0A£  (#a(F(0),V«(0)),<W0)  •  VF(fl))  =0, 

0At  <  LUW:A/i 


eg 


p(°) 


+  eg 


p(9) 


8At 


^a(r(°),  Vw(0)),  5u^  .  VF^) 


< 


624M2di62{At)2Ch 


-2 


0At 


'<")  -  <r(0)  (0)  „(»)  \  „  92(Af): 

^  ’F 


11 


eio 

r(«)  -  n-(0) 


p(0) 


eio 


jr(0) 


0At 


8 At  (  —a 


Ao) 


,<W°)  •  VF(S)) 


< 


52M2d92(At)2Ch~ 2 

4ei2 


—  (T 


AO) 


(0) 


+  en 


jri0) 


+  ei2 


_p(0) 


(A.4) 

(A.5) 

(A.6) 

(A.7) 
(A.8) 
(A.9) 
(A.10) 
(A.ll) 
(A.  12) 

(A.13) 

(A.14) 

(A.15) 

(A. 16) 

(A.17) 

(A.18) 

(A.19) 

(A. 20) 

(A.21) 

(A. 22) 

(A. 23) 
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9  A tuj  (V0)  -  a(0\F(e < 


62(At)W  ||w(0)  _  w(o)"2 


4ei3 


+  ei3 


0Afcu  -  cr(°\6u(o)VF(0)j 


< 


5292(At)2uj2dM2Ch~2 


4ei4 


-  <r(°) 


pi9) 


+  £14 


pi9) 


(A.24) 


(A. 25) 


Using  (A. 4)  -  (A. 25)  in  (A. 3)  we  obtain: 


pi9) 


1  + 


ujOAt 


A 


ei  —  €2  —  €3  —  €4  —  —  6q  —  €7  —  es  —  €9  —  €10  —  6ll  —  642  —  £13  “  £14 


< 


]+) 


2  /  1  cu202At2  52M2dCh-2L0262At2\ 

^4ei  4A2£2  4e3A2  J 


+ 


+ 


+ 


p(o) 

vr(0) 


2  /  1  Cb292At2  52M2dCh~2uj262 At2  AM2d202At2  4M2d3Ch~26292At2\ 

- 1 : — To 1 - ; — To - 1 - 1 - 


4ei  4e,5A2 


4e6A2 


£9 


£io 


2  f  M2d92At2  d2M452Ch~292At2\ 

l  4e7  4es  J 


<+)  _  CT(0) 


2  ( 92{At)2io2  5292At2Lo2dM2Ch~2' 


\  4ei3 


4e 


14 


+ 


92At2 

4£11 


a(0)  _  0.(0)  (0) 


9At 


+ 


52M2d92At2Ch~2  ||^(0)|'2 
4fl2 


(A. 26) 


Choose  fj  =  1/24,  for  i  E  {1, . . . ,  12}.  Then  using  interpolation  properties  and  lemmas  C.17  and 
C.20  with  (A. 26)  gives: 


pi9) 


<  Ch2m+ 2  o-W 
+  Ch2m+2  H"-(0) 


m+1 
2 


+  C  At2  h2m+2 


A9) 


C52  At2  h2m 


+  C  At2  h2m 


+  CAt4S2h~2 , 


A9) 


m+1 

2 


+  C  At2  h2m+2 


r(0) 


m+1 
2 


rW 


C52  At2  h2m 


m+1 


r(0) 


m+1 
2 


m+1 


»n+l 


+  CAt252h 


2  jr2  l  2m— 2 


r  m 


m+1 


+  CAt4  +  C62At2h~2 


,(°) 


(A. 27) 


_p(0) 


<  Ch2m+2  +  CAt2h2m+2  +  C52At2h2m 

+  CAt2h2m  +  C52  At2h2m~2  +  CAt4 
+  C52At2h~ 2  +  CAt452h~2. 


(A. 28) 


Lemma  A. 2  For  sufficiently  smooth  solutions  u,  a,  and  p, 


pi9) 


<S(h,At,5),  and 


u 


(9)  .  \/pi9) 


<  S(h,At,6), 


(A. 29) 
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where 


S(h,  At,  5)  = 


C[h2m+ 2  +  A  t2h2m+2  +  52At2h2m  +  A  t2h2m 

+52At2h2m~2  +  A t4  +  52At2h~2  +  CAt4S2h~2]  (l  +  At2  +  AtS ) 
+  C  ( h2m+ 2  +  A t2h2m+2  +  A tdh2m+2) 

+  C  (A t2h2m  +  A tdh2m)  +  C(At)4  +  C52  At2. 


Proof:  Subtracting  ((4.5)  with  known  u  and  n  =  0)  from  ((3.2)  at  t  =  9): 


(1  -  26)At 


ea\T)  +^(ea\TS(o))  +A(w0-  Ve^T^)  +A 


,  V«(^),- 


A  ((e)  \  ((e)  \  .  ,  ( v{6)  ~  cr(0)  0) 

(1  -  2B)At  ’ • T)  -  U  (e°-  ■  V>)  +  M  (1  -  20)  At  ~ 


+  A  (-<3-Je),5«(e)  ■  Vt)  -  W  (it"  -  Te,Tl(h')  .  ( 

Observing  that  e ^  =  F1'^  +  and  choose  r  =  F^'1 .  Then  from  (A. 30)  we  have 

p0)  2+ 

+A(1-2«)A t  l^ip  +  a(l-2fl)At  (p{i)  Su{i) ,  VF(^ 

(r<»).  F(«))  +  ALl  AA  (r<»),  W»>  •  vf(«)) 

+  (1 -20)  At  (V  •  VF(f?),F(0))  +<5(1 -20)  At  ue  •  VFW  2 
+  (1  -  20)  At  (V~  •  Vr(e'},  +  (1  -  20)  At  (y?  •  VT ^,Su^  •  VF(e)) 

+  (1  -  20)  At  ^(fW.Vw^),^)  +  (1  -  20)  At  ^(ff^Vu^),  Su®  ■  VF(^ 
+  (1  -  20)  At  ^ga(r^,Vu^),F^j  +  (1  -  20)  At  ^a(r^,  Vw(^),  <J«W  •  VT^) 
=  irW'j  +  fr(0),  F^) 


(A.30) 


•  vf(0) 


u  (1  ~  20)  At  ^(0)^(0)^  _  ^(1  ~  20)  At  (F(0)tSu(0)  .  vf-(«)J 
^(1  ~  2g) Aj;  ^r(0;_p(0^  _  ^(!  -26)  At  fy(e)^§u(0) .  VF(0)^J 

+  (1  -  20)  At  +  (1  -  20)  At  (-<rf\8u®  •  VF^) 

-w  (V  -  a9,  F ^  +  5u{§)  ■  VFi§))  . 


Bounding  all  the  terms  in  (A. 31), 


<  1  r(0~)  2  +  ei  f(0)  2, 
J  4ei 


a)  (1  -20)  At 


\Su^-VF(0)j  =  0, 


(A.31) 


(A. 32) 


(A. 33) 
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“(1  ^)At  (r®,F(fl))  <  “2{1  x2^At2  |r(0~)||2  +  e2|.F(g)l|2,  (A. 34) 

fi^-A2fl)At  (r(g),Wg)  •  vf(6>))  <  ^V4fAtl|r^l|2 

+  e3-5  (1-20)  At  wW-VF(0)  2,  (A. 35) 

+  (1  -  20)  At  (V  •  VF{e),  F(d =  0,  (A. 36) 

(1  -  20)  At  (V  •  vrW,  F(^)  <  ^  ~  At2M2d  vr(«)  2  +  £4  f(«)  2;  (A. 37) 

V  /  4e4 

(1  -  20)  At  (V  •  vr^,  <5uW  •  VF(^  <  vrw  2 

+e5-5  (1  -  20)  At  u ^  •  VF(®)  2  ,  (A. 38) 

(1  -  20)  At  ^a(F(0),Vw(0)),F(e))  <  4(1 -20)  AtMct  F(0)  2,  (A.39) 


(1  -  20)  At  [ga{F^\Vu^),8u^  ■  VF&'j 


<  e6  F^  + 


2  4-52  (1  —  20) 2  At2M2d2 


t(fl)  •  VF(e)  ,  (A.40) 


(l-20)At(ffa(r^,V-a(0')),F^)  <  ^ — 20)^  ||r(e)||2  +  e7  |f^|2  ,  (A.41) 


(1  -  20)  At  (ga(T^\  Vm(^),  6u&  • 


2r/'2  II  ,^l|2 


4-5(1-20)  A tM2d 


+  <5(1-  20)  Ate8  t*W  •  VFW  ,  (A.42) 


<_L  pW  2  +  eg  F(0)  2;  (A.43) 

(,rW)jr(»A  <  J_  r(0)  2  +  eio  F {§)  2,  (A.44) 

V  /  4eio 

"  (1  At  ^  (1  ~f)2  At'2  fW  2  +  €U  ",  (A.45) 

A  V  /  4A2en 


u;  (1  -20)  At  (F{f>),5u{j))  .  VJp(g)^  <  ^  (1o7  2g)Af  l|F(e)||2 


+  (5(1-20)  Atei2  u^-VF^ 


(A.46) 


w(l  -20)  At 
A 


>,  F^")j 


w2(l  -  20)2  At2 


(A.47) 
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to  (1  -  26)  At 
A 


r0\8u®  •  vf(^)  < 


8to2  (1  -  26)  At 

A24ei4 

+ei4<y(l  -  2 6)  At 


p(0) 


u 


0)  .  \7F0) 


(1  -  26)  At 


-  a 0 ) 


-  cri^  F0) 
(1-26)  At  4  ’ 


< 


(1  -  26)2  At2 
4ei5 


+ei5  F® 


C T0 )  -  <T0) 


(1  -  26)  At 


—  a 


0) 


Using  u |  ,)(1  =  0,  and  V  •  u  =  0, 

(l-26)At(crf\5u0)-VF^  =  -  (l  -  26)  8 At  (y  ■  u0\  <t{9)  :  F^ 

-  (1  -  26)  8 At  [u^  •  Va(9\  F^ 

=  -(1  -26)8At(u0)  •  Vcrf^F^) 

(l-  26)282At2 


Lastly 


and 


tO  \ 


(1  -  26)  At  (ae  -  ae,F (e))  < 


4ei6 


w2(l  -26)2  At2 

4ei7 


u W  ■  Verf } 


„9  rr® 

or  —  o 


+  ei7 


+  ei6 


JP0) 


F0) 


to  l 


(l-26)At(a§-ae,8u^-VF^)  < 


82to2  (1  -  2 6)2  At2dM2Ch~2 


+  ei«  F ^ 


4ei8 

2 


0 

a 

—  cr 

(A.48) 


(A.49) 


.  (A. 50) 


(A.51) 


(A. 52) 
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Using  the  bounds  stated  in  (A. 32)  -  (A. 52)  with  (A. 31)  yields 


p0) 


1  + 


u>(  1  -  20)  At 

A 


e\  —  62  —  €4  —  —  ej  —  eg 


—€io  —  en  —  ^13  —  ei5  —  ei6  —  epsl7  —  epsl8  —  (1  —  28)  A  t4Md\ 


+ 


< 


u(0)  .  <5(1-  28)  At  j^l  -  e3  -  e5  -  e8  -  e12  -  e14 

46  (1  —  20)  AtM2d2\ 

«  ) 

2  ( J_  u2(l  ~  28)2  At2  6uo2  (1  —  28)  A  A 

1  4e9  +  4A2en  +  A24ei2  / 


+ 


+ 


pi9) 

p(®) 


2  /  J_  cu2(l  -2fl)2At2  5lo2  (1  -  28)  At 
y4eio  A24ei3  A24ei4 

2/  1  A2  {1  -  20)2  At2  6A2  (1  —  28)  At 


46!  + 


A24e9 


+ 


A24e3 


+ 


+ 


vr^ 


4(1-  28 )2  At2M2d2  46  (1  -  20)  AtM2d2  \ 

£7  es  J 

2  /  (1  -  2 8)2  A t2M2d  6  (1  -  28)  AtM2d\ 
y  4e4  +  4e5  J 


rr9  rr9 

(J  —  CT 


2  (u2  (1  -  2  8)2  At2  <52a;2  (1  -  2  8)2  At2dM2Ch~2' 


4ei7 


4eig 


(1  -2d)2  At2 

4ei5 


+ 


(rW  -  o-W 


(1  -  20)  At 


—  CT 


0) 


+ 


(1  -  28)2  62At 


4eie 


2  -*>2  A  +2  ,,  -  ||2 

UU)  .  Vcrje)l 


(A. 53) 


Note  that  choosing  et(i  £  {1, . . . ,  18})  and  At  sufficiently  small  the  L.H.S.  of  (A. 53)  will  remain 
positive.  Thus,  using  interpolation  properties,  and  lemmas  C.19,  A.l,  and  C.20  we  obtain 


where 


<  S(h,  At,  6) 


and 


u(9)  .  yp0) 


<  S{h,At,6), 


S(h ,  At,  6)  =  C  [h2m+2  +  A  t2h2m+2  +  62At2h2m  +  A  t2h2m 

+62At2h2m~2  +  At4  +  62At2h~2  +  A t462h~2]  (l  +  At2  +  At6) 

+  C  (ti2m+2  +  A t2h2m+2  +  At6h2m+2) 

+  C  (A t2h2m  +  A t6h2m)  +  C(At)4  +  C62At2  +  C62At4h~2.  (A.54) 
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A. 2 


E(e) 


and 


E(§) 


Lemma  5 


E(e) 


2  |  02(1  -  a) At 
+  Re 


d(L;(0))  <  Ch2k+ 2  +  C(At)4  +  CAth2k. 


Proof:  Subtracting  ((4.3)  with  known  er  and  n  =  0)  from  ((3.3)  at  t  =  6 At)  yields: 

Re  /  (o)  \  _  /  u ^  —  u 


^  (<4?, v)  +  2(!  -  «)  (d(en),  d(v)) 


eu  +Re 


(0) 


OAt 


—  u.  ,  v 


t  >  V  > 


Recall  that  =  A^  +  E^.  Let  v  =  E^9'1  and  from  (A. 56)  we  obtain 


E (0) 


2  |  2(1  —  a)9At 
+  Re 


d(£(e)) 


=  („<o>  -  «J»,  J5«)  +  *At  (  -  u?\E^ 


OAt 


2(1  -  a)6At 


Re 


d(A{0)),d(£;(e))) 


3 

<  - 
“  2 


«(°>  -  ni0) 


2  1 
+  6 


2  3d2  At2 

+  — ~ — 


U^d)  —  -ui0) 


+  2 


A(0) 


2  1 
+  6 


E(e) 


2  (1  —  a)6At 


Re 


OAt 

d(A(e)) 


—  u 


(0) 


1  II  ,„n  1 1  2 

+  6 


2  |  (1  -a)0At 
+  Re 


E(d) 

d(L;(0)) 


Simplifying  expression  (A. 57)  we  obtain 


E(d) 


2  |  2(1  —  a)6At 
+  Re 

(o)  _„(o) 


d(£(0)) 


<  3 


uK  —  u 


+  3d2  At2 


id6*)  —  1£,(°) 


—  u 


(0) 


Make  note  that 


+3 


t*(°)  -  u® 


A^ 


2  2(1  -  a) 6 At 


Re 


OAt  t 

1 2 


d(A(0)) 


A(°) 


<Ch2k+2\\u(-,0)\\k+1<Ch 


2k+2 


Applying  lemma  C.18  and  interpolation  properties  completes  the  proof. 


Lemma  A. 3 


E ^  <  Ch2k+2  +  (1  +  h~2At)  \ch2k+ 2  +  C(At)4  +  CAth2k]  +  C(At)2h 


■2k-2 


(A. 55) 

V  v  e  Zh. 
(A. 56) 

|2  .  (A. 57) 

(A. 58) 

(A. 59) 

■ 

(A. 60) 
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Proof:  Subtracting  ((4.4)  with  n  =  0)  from  ((3.3)  at  t  =  6 At): 


Re 


(1  -  29)  At 


eu^)  = 


Re 


JO)  v 

eu  )V 


(1  -  29)  At 
2(1  -  a)  (d(efP),  d(v))  +  Re 


u(0)  -  u(8)  m  \ 

(A.61) 


Let  v  =  then  (A.61)  gives 

2  =  +  (aw,E*J  -  (a(°\e^  +  (1-29)  At 

2{l-a)(l-29)_A t(d(Jsm)id(£r 


{1-26)  At 

h\  2(1  —  a)  (1  —  29)  At 


U{0)-U(9)  -u?\E9 


Re 


4 

<  — 
“  12 


E9 

2 

+  3 

E{9) 

2 

+  3 

A(«) 

2 

+  3 

AW 

Re 


+  3  {1-29)  At 


(d(A(S)),d  {E°j) 
u<«>  -  „(*) 


(1  -  29)  At 


—  u 


(o) 


2(1  —  a)  (1  —  29)  At 
Re 

2(1  —  a)  (1  —  29)  At 
Re 


d{Ew) 

d{E9) 

d{A(9)) 

d{E9) 

(A. 62) 


Note  that 


d(i;H|  <  VE® 


<  Ch 


-l 


Considering  the  terms  on  the  R.H.S.  of  expression  (A. 62)  we  have 


2(1  -  a)  (1  -  29)  At 


Re 


d{E W)  d  {E9 


< 


2(1  -  a)  (1  -  29)  At 


Re 


d(£(e)) 


*  3C  h 


2 1,-2 


2{1  -  a)  {1  -  29f  At  \  (2{l-a)9At\ 

Re  ) 


9Re 


d{E 


W' 


2  1 
+  12 


Ch-1 
2 


Eb 


E 


<  3C2/r2  ( 2<1  -  °)  t1  -  M)~ At 

\  ORc 


Ch2k+2  +  C{At)4  +  CAth 


2k 


1 

+  12 


Eb 


(A. 63) 


and 


2(1  -  a)  (1  -  29)  At 
Re 


d(A^)  d{E9) 

2 


\  Re 2 


Ch2k  +  — 


Eb 


(A. 64) 


Applying  lemma  C.19,  with  the  bounds  established  in  (A. 63),  and  (A. 64)  to  (A. 62)  and  we  obtain: 


E (0) 


<  Ch2k+2  +  (1  +  h~2At)  \ch2k+ 2  +  C{At)4  +  CAth2k 1  +  C{At)2h 


2k— 2 


(A. 65) 
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B  Taylor  Expansions 


The  following  are  expressions  for  u  at  different  times.  Specifically  the  expressions  listed  are  com¬ 
putations  of  Taylor  series  about  the  expansion  points  uin+  2  \u(n+^~9\  and  u^n+e~^>  that  give  an 
exact  value  for  the  truncation  of  the  series.  Integration  by  parts  proves  each  of  these  statements. 

Expanding  about  Jn+5): 


i(n+1)  =  n(n+^)  + 


r*(n+ 1) 


dt 


("+2) 


=  u(n+  2)  +  tu[n+^  +  f  utt(-,t )  (t(n+ 1)  -  t)  dt 

=  u(n+^  +  tu[n+2^  +  +  ]-  f  um(-,t )  (t(n+i)  -  tf  d1{BM) 

“  '  *(»+*) 


u(n+0) 


Jn+i)  -  f^n+hKt(;t)dt 

n+e ) 

u(n+h)  -  -q  \  A tu[n+2^  +  f  utt{-,t)  (t-t{n+e))  dt 

\z  /  •/*(„+«) 


(B.67) 


dt  (B.68) 


u 


H  =  „("+£)  -  A"+*)  mb 


dt 


J(n) 


=  u 


W  -  (1)  Ai  4"+4)  +/'("+’)  ««(•.*)  (<  -  *<»)) 


dt 


2  4 

_1  f  ("+2)  Uut{-,t)  (t  —  t(n))2  dt 
z  dt(n) 

Expanding  about 

ti(")  =  n(n+i-0)  -  [t(n+h~eKt(;t)dt 

A") 


(B.69) 


=  ti("+2  e)  -  Q  -  At ?4 ”+  2  ^  (',1+2  ^  ua(-,t)  (t  -  t(n))  dt  (B.70) 
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u 


(n+0)  _  u(n+h~6 

=  u(H+ 3“6 


*(«+») 


ttt(-,t)  dt 


(n+h~d) 

+  (  26-  ^ 


r*^(n+0) 


(n+^-6») 


Utt(’,t)  (t(„+0)  -  t) 


u(n-e)  =  u{n+\~e)  _  /*("+*  »)  u^.^dt 

dt(n-9) 


=  U(n+5-®)  -  —  Jn+3- 


^  +  f  ^+2  ^  ««(•,*)  (i  -  i(n-e))  dt 

dt(n-6) 


U 


(n+~e)  =  u(n+^~e)  + 


t(n+S ) 


ut(-,t)  dt 


(n+^-6>) 


=  «(n+2  e)  +  ^4"+a  ^  +  [{+)  utt{-,t)  (f(n+e)  - 1)  dt 


(n+^-e) 


Expanding  about  «("+*“ 5): 


«(")  =  «("+M)  +  [Hn) 


ut(-,t)  dt 


{n+e~h) 


=  u 


(n+e  2)  +  Q  -  At4n+0  +  Jt  )  (t(n)  -  t)  dt 

{n+e~h) 


uy 


(■ n-0 )  =  (n+0-1)  _  /’i("+«-2) 


«t(-,t)  ^ 


l(n-0) 


=  u(n+d  5)  -  ("20  -  ^  At4n+0  ^  ^ n+6  ^  utt(-,t )  ( t  -  t(n_e)) 

V  AJ  Jt(n-e) 


u 


(n+0—1)  _  u{n+e-\) 
=  u(n+e-^ 


{n+e~h) 


’t(n+e-i) 
U 


Ut(-,t)  dt 


[n+(>  ^  +  /  ^+S  ^  utt(-,t)  (t  -  t(n+e_ i}) 

Jt(n+B- 1) 


u 


<•**>  =  u(^-i)+  [Hn+$)  Ut(;t)dt 

Jt / 


(n+0—  ^ ) 

=  „("+»-§)  + (A)  4”+"4) 


/■*("+») 

/  uu{-,t)  (t(n+g) -t)  dt 

Jt,  1  \ 


dt  (B.71) 


(B.72) 


(B.73) 


(B.74) 


dt  (B.75) 


dt  (B.76) 


(B.77) 
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Expanding  about  u(n  2); 

(n)  A  t  (n—  i)  1  (At)2  (n—  i) 

u[n)  =  uyn  *>  +  —u\  2>  +  -K-^u\t 

Cn_i'i  (n—  I'l  At  (n—  4)  1  (At)2  (n—  A 

w(n  1)  =  [n  2) - V  2)  +  _  V - )_  \  2) 

2  t  2  4  44 


uttt(-,t)(tn  -  t)2  dt,  (B.78) 

-tn— i)2  dt  (B.79) 
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C  Lemmas 


Lemma  C.l 


0it(n+1)  +  (1  -  0)  u{n+e)  -  u(n+5) 


2  (At)3  />t("+1)  o 

<  /  IM+t)f  dt 

*(»+«) 


24 


(C.80) 


Proof  of  Lemma  C.l  using  (B.66),  (B.67)  and  choosing  the  value  of  0  =  1  —  : 

0u(n+1)  +  (1  -  0)  u(n+e)  -  u(n+5)  2  =  J  ^u(n+1)  +  (l-6»)u(n+e)-u(n+5)j  dA 


9  |  u(ra+2)  +  -AtUf  +  f  utt{-,t)  (t(n+ 1)  —  t)  dt 

Vi) 


+  (1-0)  (u("+4)-  Q-f)  Atu[n+h\  jy+hKtt(-,t)(t-t{n+e))  d?j  dA 


-6>2  +  20-  -  )  At 4  2> +9  utt{-,t)  (t(„+i)  -  t)  dt 

J  Vi) 

2 


t  \ 

+  (1-0)  f  +2"*  «tt(-,t)  (t-t(n+0))  dt  J  dA 

dt(n+e)  / 


^(ra+1)  /*£/  ]\  \ 

0  /  Utt{-,t)  (t(n+i)  -  t)  dt  +  (1  -  0)  j  +2’  utt(-,t)  (t  -  t(n+e))  dt  )  dA 
(n+s) 

2 

*("  +  !) 


'*(«+») 


9  />t("  +  l)  o 

(uu(-,t))  dt  (t(„+i)-t)  dt 

Jt,  x 


<  2  /  I  6»2  |  /  u«(-,t)  (t(„+i) -t)  dt  |  +(1-0)J(  f  *'"+2')  utt(-,i)  (t-t(n+e))  dt  ]  )  dA 

'*(n+9) 


(n+s) 


Vi) 


+  (1  -  0)2  |  f  +2')  (utt(-,t))2  dt  j  ^  +2-*  (t  -  t(„+e))2  dt  j  j  dA 
VW)  dt(n+s)  j  j 


t))  dt  dA 
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(««(•,*))  dt  |  +(i-0)L(i-2 ey 

*(»+!) 


/‘("+S)  (».«(•■  t)fdt 
'*("+»)  ) 


vn  1  1  1  (utt(-,t))z  dt  ]  dA 

I  Cl 


(1  -  6>)2 

12  v  ; 


f‘(n+l) 


||wtt(-, t)\\2  dt 


*(»+«) 


dA 


Lemma  C.2 


6ht(n)  +  (1  -  0)  u(n+§)  -  u(n+5) 


2  <  (At)- 


24 


£(n  +  0)  « 

\\utt(;t)\\  dt 


*(») 


(C.81) 


Proof  of  Lemma  C.2  using  (B.68),  (B.69),  and  the  value  of  9  =  1  — 


+  (1  -  6<)  u(n+e')  ^Jn+5)  =  J^0u{n)  +  (1  -  9)  u(n+§)  ~  u(n+^  dA 
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<  2( ^ 
-  3  I  2  , 


> 

max  |d2,  (1  —  9 )2  ,  (1  —  20)c 


'("  +  »)  9  \ 

(utt(-,t)Y  dt  dA 


n  \  Jt 


(n) 


12 


Lemma  C.3 


{At)  (1  -ef  [t{n+e)  \\utt(;t)\\2  dt 


*(»> 


6u{n)  +  (1  -  9)  u{n+e)  -  M(n+5-e) 


2  <  (At) 

-  24 


3  ft 


(n+e) 


IKtM)ll  dt 


(C.82) 


*(»> 


Proof  of  Lemma  C.3  using  (B.70),(B.71)  and  choosing  9  =  ^1  — 

9u{n)  +  (1  -  9)  u{n+6)  -  J"+w“e)  2  =  J  (Wn)  +  (1  -9)u{n+0)  -w(n+3- 


dA 


0  ^(n+*-fl)  -  Q  -  o)  At4n+i  ^  +  /t("+i_S)  ««M)  (*  -  *(»*)) 

+  (1-0)  ^u("+3-e)  +  ^20-^  Atr&n+*-e)  +  J‘ln+B>  utt(;t)(t(n+e)-t)  dij  -J"+M)j  dA 


-02  +  29  -  A tuin+i~e)  +  9  [  ("+^e) 


+  (!-0)^  ««(■,*)  (i(n+e)  -  0  dt  \  dA 


r ,  I - t  W  V  I  utt('i  t)  (t  ^(n))  dt 

Z  '  JHn) 

2 


(n+h~e) 


utt(-,t)  (t(n+0) -t)  dt  j  dA 


(r.+  1-e) 


««(-,*)  (t  -  t(n))  dt )  +(1-0)" 


f*(n+e) 


^^(*5^)  (^(n+0)  dt  j 


>+!-) 


(Utt(-,t))  dt 


(~M  (*-*,„>)+<’ 


Hn) 
ft(n+S) 


2  7  ("+0)  2  \  \ 

(utt(-,t))dt  (t{n+e)-t)  dt  dA 

Jt,  ,  X  /  / 


("•■2  ") 
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uws-'H'fr” . . 


+  (l-0)2  ((20- 0  At)  j*™  (utt(;t))2  d?JdA 

I  (t)3  /  (1  -  20)3  r("+M  t))2  dt  +  (1  -  0)2  (40  -  l)3 


(utt(-,t))  dt  dA 


*(n  +  |-9) 


max  |02,  (1  -  20)3  ,  (1  -  O)2  ,  (40  -  1)J 


{utt(-,i)f  dt  dA 


w: 

ygv-vr” 


(' Utt{-,t ))  dt  dA 


\\utt(;t)f  dt 


Lemma  C.4 


-  20)  u(n+§)  +  0u<n>  +  0M(n-e)  -  w(n+5- 


2  <  (At)3 
“  3 


(3-2V2 


f  (n+6)  \\utt(;t)\\2  dt  (C.83) 
*(»-«) 


Proof  of  Lemma  C.4  using  (B.70),  (B.72),(B.73)  and  choosing  0  =  (l  — 

(1  -  20)  u(n+e~)  +  0u(n)  +  0u{n~e)  -  u(n+5"e) 

=  f  ((1-2 0)  uf+§)  +  0u{n)  +  Ou^-8'1  ^  u(n+^-s)y  dA 

J  Q 

=  |^(l  -20)  +  M;t)(t{n+g)-t)  dij 


+0  |u(”+2  ^  -  Q  -  d'j  Atu[n+2  ^  +  J^n+2  ^  utt(-,t)  (t-t(nf  dtj 
+e  ^(”+5-*)  -  ^4n+i-e)  +  utt(;  t )  (t  -  t („_,))  ettj  -  u(n+l>-e)  j  dA 


((1  -  20)  -  0  (1  -  20)  -  0)  ^t4"+i“9)  +  (1  -  20)  f  ("+S) 

Z  1  \ 

(»+*-•) 


««(■,*)  ( 


,  s\  —  t )  dt 


)2 

dA 

J  ((1-20)  J(  }  utt(-,t)  (t(„+ei)-t)  dt  +  0  J  ^  +2  ^  utt(-,t)  (t-t(n))  dt 
n  \  tf+h~s)  t(n) 

t  \  2 

+0  f  ^  +2  ^  utt(-,t)  f  -  t(n_g))  dt.\  dA 

■'*(»-«)  J 
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<  4  /  (1-20) 

In  \  \  Jt 


n+S) 


(n+|-e) 


««(■,*)  (i(n+e)  -t)  dt  |  +0 


2  /  /**("+ 5-#) 


««(•,*)  (*  -  i(n))  dt 


*(n) 


+02  ^  ^"+2  ^  (t-t{n_e))  dtj  j  cL4 

<  4  /  (\l  —  20)2  f  (utt(;t))2  dt  I  in+"]  [t 


("+5-«) 


(*("+*)  *) 


—  £ )  dt 


+e2  /  ("»-)  („„(,() )2  dt  /‘f-1-)  (t-t(„,)2  dt 

Jt{rt)  dt(rl) 


+  0! 


2  /vm 

'i(n-e) 


(««(-,  *))2  dt  f  ^  •*  ( t-t(n_6) )2  dtj  rM 

dt(n-e)  ) 


'  (1-2 0)2  T("+e)  (u«M))2  di  +  02  f(l-20)^V  /t(n+i_e)  («*,(-, t))2  dt 


M  ») 


(») 


(««(■>*))  d*  d-4 


2  /  /  /  jv\2  j,  ,  n2  fi  om3  f  (n+2  8)  /..  /  a\\2 


(«t*(-,t))^  dt  +  02  (1  -  20)'1 


(utt(-,t))  dt 


(»+*-*) 


t(„) 


+02  f  ^  (itM(-,f))2  dt  J  gL4 
■'*(«-«)  / 

<  ^  2max  |(1  —  20)2 , 02,  (1  —  20)3|  J  \utt(-,t))2  dt\ dA 


(t>-8) 


^(i-20)2  /t("+e)  IKM)II2^ 

d  An-8) 


Lemma  C.5 


01t(n)  +  (1  -  0)  W(n_e)  -  u("+M) 


2  (At)3  /■*<«)  9 

I  \\uu(;t)\\2  dt 


24 


H»-9> 


Proof  of  Lemma  C.5  using  (B.74),  (B.75)  and  the  value  of  0  =  ^1  —  . 

0u(n)  +  (1  -  0)  u(n"e)  -u(n+0-5)  2=  t  (0u<n>  + (1-0)  «("“*> -u(n+0_*))2  dA 

J  £7 
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J  (Wn)  +  (i  -  e)  u{n~0)  -  u(n+e~ I ) ) 2  dA 

(n+O-l)  +  (1  a\  Atu("+<>-k)  +  /'<(")  dt) 

V  7  V-i)  / 

+  (1  ^  0)  (U"+e^)  -  ^2d  -  0  A tu[n+e~h)  +  M;t)  (t  -  t(n_fl))  dfj  -  dA 


02  _  20  +  ^  Atot("+^} 


'(»*) 


+  6  Uu{-,t)  (t(n)  -t)  dt 


2 


+  (1  -0)[(+  Utt{-,t)  (t  -  t(n_g))  dt\  dA 
JHn-o)  J 


f(  +  ^  uu(-,t)  (t  -  t(n-9))  dt  \  dA 


Utt( •,  i)  (t(n)  -t)  dt  +  (1  -  oy 


\^-i) 


un{-,t)  (t  -  i(n-e))  dt  ]  )  dA 


r(n) 


(t(n)  t^)~  dt 


{utt{-,t))2  dt  f  (  +  2-*  (t-t(n_g))2  dt  \  dA 

dt(„-9)  / 


r<») 


(««(•,*))  df+(l-d)2(4d-l) 


3  /■*(»+•- J) 


(utt(-,t))2  dt  dA 


V'-i) 


r("-8) 


:(n) 


(utt(-,t))  dt  dA 


/-.  m2 


rW3 


(i  -ey  \\utt{.,t)\\  dt 

^  Jt(n-e) 


Lemma  C.6 

1(1  -  2d)  U(n+e-1)  +  dM(n)  +  6ht(n+e)  -  m("+0-5) 


2  <  (A0 


(3-2V2)  [kn+6)  II utt( 


,  i)||2  dt 


Proof  of  Lemma  C.6  using  (B.74),(B.76),(B.77)  and  the  value  of  6  =  ^1  —  : 

(1  -  2d)  u('n+e~1)  +9u(n)  +0u{n+9)  -u(n+e-5)  2  =  J  ^(1-2 :d)w("+0-1)+dw(")+dw(”+e)-w(”+^?)j  dA 
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f  (( 1  -  20)  (u(n+9-^  -  U„(;t)  (t  t{n+e_ !))  di) 

+0  ^u(n+e_5)  +  Q  -  0^j  Atw^"+e  ^  +  Jt  i  utt{-,t)  (t(n)  -  t)  dt. j 

+0  u[n+e-l)  +  jHn+6)  uu(;t)  (t{n+e)-t)  dt'j  -u("+e-5)NJ  dA 

[  (  (-202  +40  -l)  4"+0'l)+(l-20)  ft{n+e~h)  utt(;t)(t-t{n+e_ d)  dt 

Jfl  y  V  /  J  t  (n-H-0  _  1) 


/  Utt(-,t)  (t,(n)  -t)  dt  + e  [  Uu{-,t)  (t(n+g)  ~  t)  dt\  dA 

Jt{™+e-h)  Jt(r>+o-h)  J 


(»+•-*) 


f  [  (1  -  20)  f  ^+<>  ^  utt(-,t)  (t-i(n+0_i))  di  +  0  f  (t(„)  ~t)  dt 

V  0t(  n  +  9-1)  ''t(n+»-4'\ 

rHn+9) 


t(n+9-  1) 

\2 

^ti(*>^)  (^(n+0)  0^  J  0^ 


(»+»-4) 


V-i) 


<  4  /  (1  —  20)2 


Utt{-,t)  (t-t(n+g_1))  dt  +0 


(^-4) 


««(■,*)  (*(n)  -  0  dt 


^ti(*)^)  (^(ro+0)  dt  dA 


(»+•-*) 


<  4  /  f  (1-20)2  f  ^  ^  {utt(-,t))  dt  (  ('"+8  2-*  (i-t(n+e_i))2  dt 

J  Q  \  J  t(n+0-l)  J  t(n  +  d-l) 

/t(n)  r^i71) 

(utt(-,t))2  dt  J  (t(n)  -  i)2  dt 


V 

*(n  +  0-l) 

+02 

f 

Ott(-,£))2  dt 

l+e-i) 

+02 

/" 

n  +  d)  9 

(««(',*))  dt 

‘+*-4) 

V’-l) 

/■*(»+«) 


/•*("+®)  .  ,2  \ 

/  (i(n+0)  -  i)  dt  J  dA 

Jt("+«-4)  ' 


(1  -  20)2 


/“(»+«- i) 

't(„  +  8-  1) 


“  3  1  2  , 


*(“+•-4)  ' 

<  max  {(1  —  20)2 , 02,  (1  —  20)3 


(«tt(-,*))  dt  +  02(1  —  20)2 


(««(•.*))  dt 


\^-h) 


(■ ««(■,* ))2  dt  dA 


(utt{-,t))  dt  dA 


(1  ^  20) 2 


||««(-,i)||  dt 


*("+«- 1) 
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Lemma  C.7 


,("+»)  _  ,,("+s) 


<  A t  f  ^  +2-*  ||wt(-,i)||2  dt 
-'*0+0) 


Proof  of  Lemma  C.7  using  (B.67): 


,(«+»)  _  u0*+f) 


/«+»)  _  U("+J)  )  dA 


f  (u(n+i)  -  [t^+iKt{-,t)dt-u(n+^\  dA 

J  &  \  -'*(n+0)  J 


< 


< 


*0+1) 


O) 


l2df  '  -2 


dtdA 


<  At  j  (  +2->  ||ut(-,f)||2  dt 
-'*("+«) 


Lemma  C.8 


U+e)  _  (n+ 1) 


<  A t  f  {  }  ||tit(-,f)||2  dt 

'(-I) 


Proof  of  Lemma  C.8  using  (B.68): 


i(n+®)  -  u(n+5) 


r(n+e~)  -  u(n+a)  1  dA 


u(-+l)  +  /  ‘("+9)  dt  -  «(n+»)  | 


:(-+«) 


+  /  ut(-,t)dt  —  u\ 

\*l) 

v 

1  ut(-,t)  dt  dA 


dA 


n  \  Jt 


(”+i) 


< 


< 


Af 


r 

>+») 

l2df 

«+») 

M 

'(■ 

*+3) 

•> 

+3) 

p 

J"+3- 

fl)  l2df  , 

r 

% 

*•3) 

J 

t( 

(  n+ 

s) 

r 

O+o) 

IKM)II 

2  dt 

(ut(-,t))  dtdA 


( uf(-,t ))2  dtdA 
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Lemma  C.9 


+  !«("+*)-«("+«)  <  ^  [Hn+e)  |M-,t)||2  dt 


Proof  of  Lemma  C.9  using  (B.67)  and  (B.68): 


n+e)  +  ^u(n+®)  -  Jn+s) 


I  -u(n+e)  +  ^u("+e')  -  n(”+5)  dA 
a  \ 2  2  / 


i  f  um)  -  /V*) 


ut(-,t)  dt 


+  2  “ 


-  (  U(n+5) 


+  /  ("+£,)  ut(-,t)dt  -Jn+5)  )  dA 


\  J  t  /  2  > 


*("  +  «)  +  \ 

Ut(-,t)dt —  /  v  2^  i(t(-,f)di  dA 


«t(-,t)  dt 


ut(-,t)dt  dA 


1  f  /W)  .  /■*("+») 


-  1  dt  /  dt 

JnJt(n+i)  \n+i) 

+  [t{n+k)l  dt  ft{U+i)  (Ut(;t)f  dtdA 

^^(n  +  O)  *'*(n  +  0) 

At  r  /■*(->+«) ,  ,  .  /•‘iwa , 


(ut(-,t))2  dt  +  [  +2^  (ut(-,t))2  dtdA 


At  /■>+»)„  .  ..,2  , 

<  -y  /  IKM)II  dt 

z  Jhn+<>) 


Lemma  C.10 


/  ,  \  2  /*^(n+0) 

i(n+fi)  -  u(”+5-e)  <  At  /  |M-,t)||  dt 

Jtt  ^  \ 


(~+M 


Proof  of  Lemma  C.10  using  (B.71): 


(n+0)  _  ?.(n+i-< 


t(n+fi)  -  u(”+5-e)  I  dA 
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[  Ln+h~*)+  [Hn+B) 

V  “'Vn+i-e'l  , 


dA 


r(»+i-0 

v 

1  ut(-,  t)dt  I  dA 


l2df 


IKM)||  dt 


(ut(-,t))  dtdA 


Lemma  C.ll 


i(n+§)  _  u(n+ 1 -A 


<A t  \\ut(-,t)\\  dt 

•) 


Proof  of  Lemma  C.ll  using  (B.73): 


,("+«)  _  u("+5-e) 


t("+®)  -  «("+*-*)  |  dA 


(n+3-®)  +  /  ("+e)  Ut(;t)dt-  u(n+i-e)  )  dA 

*(n+i-e)  / 


lut(-,t)dt  dA 


l2  dt  ^  (ut(-,t))2  dtdA 


IM’A)I|2  dt 


Lemma  C.12 


\(”+e')  +  -M(n+0)  -  U(n+5-fl 
2  2 


2  „t 

£At/ 


("+«)  „  ,  .x  ,,2 


IKM)II  dt 


("+3-s) 


Proof  of  Lemma  C.12  using  (B.71),  and  (B.73): 

I  2 


iM("+e)  +  £„("+«)  -  U("+M 


'  «("  • ®)  +  ^”+e)  -  u("+5-e)  )  dA 
/a  \  2  2 


1 


«("+=-0)  +  f  ("+8)  «t(-,t)dt 


(n+|-") 

+  i  /«("+*-*)  +  f("+e)  «*(•,*)  dtj  -  n("+5-e) 
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Iff  /•*("+»)  /•*("+»)  \ 

-  /  /  ut(-,t)  dt  +  /  ut(-,t)dt\  gL4 

4JnV‘H-9)  / 


1 

<  - 

2  ./n  \  \  Jt 


'(«+») 


dt 


KM 


*("+«) 


ut(-,t)dt  dA 


KH) 


(■*(»+») 


l2  d£ 


r*("+«) 


<  if  [(n+B)  1 2dt[in+e)  dt  + 

JaJt(n+k-*)  ^KM 


(ut(-,f))2  dtdA 


At 


A  (n+9) 


<  —  /  /  (  )  (uf(-,t))2  dt+  I 

JaJt(n+i-s)  !,) 


( Ut(-,t ))  dtdA 


<  At 


*(»+•')  „  ,  ,„|2 


IKM)II  dt 


'(•••*  M 


Lemma  C.13 


(n+e-l)  _  (n+0-i) 


<  At  f  (  +  ^  ||ttt(-,f)||2  dt 

Jtin+e- 1) 


Proof  of  Lemma  C.13  using  (B.76): 


(n+0-1)  _  (ra+0-i) 


< 


<  At  f  ^  +  ^  \\ut(-,t)\\2  dt 

J  t(n  +  6-l) 


Lemma  C.14 


(n-0)  _  (n+e-|) 


<  At  /  2'  |M-,f)|r  dt 


Proof  of  Lemma  C.14  using  (B.75): 


■»-0)  _  1) 


1  C  _  ]  ,/,l 


[  ( u(n+e- *)-  ^  ut(;t)dt-u(n+0-^)\  dA 


ut(-,t)dt  dA 


n  \  Jt 


<»-») 
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< 


< 


("+8-i)  ,2 


l  (it 


^  (ut(-,t))2  dtdA 


(n-0) 

t 


At 


^+S  ^  IM-,i)l|2  dt 


*("-») 


Lemma  C.15 


^(n+e-i)  +  lu(n-fl)  _  u(n+9-i) 


2 

<  At 


IM-,i)||2  (ft 


Proof  of  Lemma  C.15  using  (B.75),  and  (B.76): 

1 2 


£  (n+fl-i)  +  iu("-fl)  _  u("+*-5) 
2  2 


dt 


< 


< 


Ut(',t )  (ft  + 


«t(-,f)  dt 


'  t(n  +  0-l) 


ut(-,t)dt  dA 


U»-e) 


ut(-,t)dt  dA 


r("-8) 


(ut(-,f))2  dt 
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At 
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Lemma  C.16 
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Af 
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Proof  of  Lemma  C.16  using  (B.78)  and  (B.79): 
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tn 
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"At"  2  +2U*  _2~Utt  2Af  Jt  uttt{-A)ii-tn-i)  dt-ut  J 
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1 


1 


— —  /  Uttt{- -  t)z  dt  +  — — /  um(-,t){t-tn-{)  2dt)  dA 


2A  t 


2A  t 


< 


< 


'tn-l 

2 


Uttt{- -  t)2  dt  +2  /  2  uttt{-A)(t-tn-i)2  dt  dA 


dt  ( tn-t)4dt. 


„  i 


rt_  i 


( uttt{-,t ))  dt+  (utttiA))  dt  dA 


dt  dA 


\\uttt{- ,t)\\  dt 


Lemma  C.17  Assuming  ||utt(-,i)||  <  M  then 
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(0)  _  u(0) 

9  At 


—  u 


(0) 


<  C{At)2. 


Proof  of  Lemma  C.17:  Make  note  that 
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(tf\  an  / 

u 
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I  f  \utt(-,t)\2  dtd A 

J  J  to 
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<  ^{9At)2M2 

<  C(At)2. 
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Lemma  C.18  Assuming  \\uu(-,t)\\  <  M  then 


I  —  u 


-  u\U)  <  C(Atf 


Proof  of  Lemma  C.18:  Make  note  that 
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Lemma  C.19  Assuming  ||u#(-,t)||  <  M  then 
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Proof  of  Lemma  C.19:  Make  note  that 
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Thus 
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Lemma  C.20  Assuming  ||ut||  <  M  then 
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Proof  of  Lemma  C.20:  Make  note  that 
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